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MODERN INTRODUCTORY MECHANICS 


MATHEMATICAL REVIEW 


1 MATHEMATICAL REVIEW 


TRIGONOMETRY 


Mathematics is the language of physics, so we must all have a certain fluency. The first order 
of business is to remind ourselves of some basic relations from trigonometry. 


opposite 

sin 0 = -, 

hypotenuse 

„ adjacent 

cos 0 = -, 

hypotenuse 

sin 0 opposite 

tan 0 - -- -. 

cos 0 adjacent 


For right now just think of a vector as something with both a magnitude and a direction: 



Vector notation: A = (a : , A 2 ) . Vector addition: 
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Unit Vectors: 


3 



These are handy guys which point along the 1,2 or 3 directions with unit magnitude. We 
always choose a right-handed coordinate system in this class. The right-hand rule identifies 
such a system: curl the fingers of your right hand from e x to e 2 in the above figure; your 
thumb will point in the e 3 direction. 

Matrices 

A matrix is a collection of entries which can be arranged in row/column form: 

-► columns 

A 12 A 13 

I A 2 l A 2 2 ^23 

\A 3 i A 32 A 33 / 


I 

rows 


A single generalized matrix element is denoted: 


row 



column 


Addition of matrices (number of rows and columns the same for both A and B matrices): 


A X1 

A i2\ + 

(Bn 

B 12\ _ 

(An 

+ 

Bn 

A 12 

+ 

B 12 

A 2 i 

A 22 / 

\b 21 

B J 

\A 21 

+ 

b 21 

A 22 

+ 

B 22 


( 1 . 1 ) 
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In more abstract language (“index notation”) this is just 


(A + B)ij = A ±j + B ijf (1.2) 

where i and j are taking on all possible values independently. In the above equation i and 
j are said to be “free” indices. The free indices on one side of an equality must always be 
the same on the other side. 


Multiplication of matrices. (Here we only require that the number of columns of A equal 
the number of rows of B): 



A n B i + A 12 B 2 ^ 
V A 21 B 1 + A 22 B 2^ 


(1.3) 


Another example: 


11 ^12^ 

( B 1X 


f A ll B ll 

+ 

A 12 B 21 

A 11 B 12 

+ 

A J 
21 ■ r ^22 y 

vb 21 

B 22^ 

^ A 21 B 11 

+ 

A 22 B 21 

A 21 B 12 

+ 


A 12 B 22 

A R 

22 22 


(1.4) 


Notice that the result has the same number of rows as A and the same number of columns 
as B. In index language, these two examples can be written much more compactly as: 


(AB)i = 2 A ik B k 

U_ 


k = l 


"dummy" index, k 


(AB)ij — ^ Ai k B k j 


k = 1 


"dummy" index, k 


(1.5) 


( 1 . 6 ) 


Note that dummy indices are ones which are summed over all of their values. Unlike free 
indices, which must be the same on both sides of an equation, dummy indices can appear 
on either side. Also notice that dummy indices always appear twice on a given side of an 
equation. These rules trip up many beginning students of mechanics. 


For reference, here is a summary of the understood “index jockey” rules for index manipulations: 


1. “Dummy” indices are those which are summed. Each such index always appears 
exactly twice. One can interpret this sum as matrix multiplication only if the 
indices can be placed directly next to each other. Separate summation symbols 
must be used for independent summations. 
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2. In general, one can not change the order of indices on an object, such as A^j 
(Occasionally one knows effect of interchanging indices; see later comments on 
symmetric and antisymmetric matrices.) 

3. Free indices are those that are unsummed. In general, each free index appears 
once on both sides of a given equation. 

Identity matrix (3x3 context): 


(1 0 0 \ 

0 10 
VO 0 1 ) 

The identity matrix is often simply written as the number “1”, or is absent altogether in 
contexts where its presence is unambiguous. (Physicists must learn to read behind the lines 
for the meaning!) 

We will need three additional matrix operations. 

1. Inverse: AA' 1 = 1 (1.7) 


The “1” on the right hand side here means the identity matrix. A Theorem from linear 
algebra establishes that AA 1 = 1 implies A _1 A = 1 (Can you prove this?) Finding 
A -1 in general is fairly complicated. For most of the matrices we will encounter, finding 
A -1 will be easy. (I’m thinking of rotation matrices, which will follow shortly.) Notice: 

(AB) -1 (AB) = 1, 

(AB ) -1 A = B -1 , 

=> (AB)" 1 = B^A- 1 . ( L8 ) 

2. Transpose: A- = Ajj_ (1.9) 


Examples of the transpose operation: 


B 


B i) 

B 2 ) ' 


"column matrix" 


A = 




B T = (B, B 2 ) 

"row matrix" 
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Also 


(AB)^ = (AB) j± , 

- 

- Sbix,. 


T 


= (BA ) ij , 


T T T 

=> (AB) = BA. 


( 1 . 10 ) 


3. Determinant: 


2x2 case: det 


A 11 A 1 2 ) 

A 2 iA 22 y 


AnA 22 


Ai 2 A 2 i 


also written as 


A ll A l 2 

A 2 iA 22 


/Xl 

■^12 

A 13 ) 



x3 case: det 1A 21 

A 22 

A 23 

- A 11 (A 22 A 33 - A 23 A 32 ) 


U31 

^32 

A 33^ 



- A 12 (A 

2 1 A 33 

- a 23 a 31 ) + a 13 ( a 21 a 32 - a 22 a 31 ) . 

(1.12) 


Note that det(AB) = (det A) (det B) and that det A T = det A and det( A 1 ) = 
l/det(A). We’ll find a more elegant definition of the determinant later. An important 
point to realize is that the inverse of a matrix, A, exists only if det (A) is not zero. 

An important point about linear algebra will also be called upon in later chapters. A system 
of linear (only x x 2 3 appear, never ( x i, 2 , 3) 2 or higher powers) homogeneous (the right 
hand side of the following equations are zero) equations, 



A 11 X 1 A 12 X 2 

+ a 13 

X 3 = 

0 , 


- 

A 21 x 1 + A 22 x 2 

+ a 23 

X 3 = 

0 , 



A 31 X 1 A 32 X 2 

+ a 33 

X 3 = 

0 , 

(1.13) 

has a nontrivial ( x : f 2 f 3 

are not all zero!) solution if and only if 



( A n A 12 

A 13 ) 





det I A 21 A 22 

A 23 

= 0 

• 

(1.14) 


\A 3 i A 32 

A 33' 
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One often encounters matrices which are said to be symmetric or antisymmetric. A symmetric 
matrix is one for which 

A t = A or Aj i = A ± j 

An antisymmetric matrix has 

A T = -A or A^ = -A^ 

Matrices are not generally symmetric or antisymmetric, but such combinations can always 
be constructed. For example, 


C = A + A t 


is a matrix combination which is symmetric (f-i j — C ji) even though A is not itself. Likewise 


D = A - A t 


is antisymmetric (Dij = - Dji) 
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Orthogonal Transformations 

Let us study a little about the above mentioned rotation matrices. Special case: 

2 


P 



1 


P at (x 1 , 0) 
(old coordinates) 


Let’s do a “passive rotation,” where we rotate the axes and not the point P. (Define 0 >0 
for counterclockwise rotations.) 



(0 >0 here) 


P at (x 1 ,x 2 ) = (Xj^cosB, -x^inB) 


= (Xj^cosO, Xjcos (0 + — ) ) (new coordinates) 


Another special case: 
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P ' at ( 0 , x 2 ) (old coordinates) 


Tt 


P' at (x 1 ,x 2 ) = (x 2 sin0, x 2 cos0) =(x 2 cos( j-0),x 2 cos0) (new) 

Let’s now put these two special cases together into the general case of rotation in a plane: 


P'» 


x 2 p" 


x 


P" at (X x ,X 2 ) 
P P' 

: (*!,0) + (0,x 2 ) 
(old coordinates] 


Now characterize the point P" with respect to rotated axes: 



p P' 

P" at (x lf x 2 ) = (x lf x 2 ) + (x 1 ,x 2 ) 


JT JT 

= (x-^cosQ, x 1 cos(0 + — ) ) + (x 2 cos( — - 0) , x 2 cos0) (new coordinates) 


Therefore, comparing the old and new descriptions, we have 


, Jt 

x x = Xj cos 0 + x 2 cos (-0) , 

2 , 
TC f 

x’ = X, COS(— + 0) + X cos 0. 

2 


(1.15) 


To simplify the notation, introduce the concept of direction cosines: 
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cos(x! L ,Xj) = cosine of the angle between the x ± and Xj axes. 


start measurement of angle 


The angle may be measured in a clockwise or counterclockwise sense but must be consistent. 
We will adopt counterclockwise as in the above examples. For the above rotation, 


fcos(x 1 ,x 1 ) 

| cos (x' lf x 2 ) 

COS(X 2 ,X 1 ) 

cos (x 2 , x 2 ) 


cos 0, 

TC 7t 

cos(2it - — + 0) = cos(— - 0) , 

. TC 

COS(— + 0) , 

2 

cos 0. 


So, our general rotation may be written 

x x = x : cos (x lf x x ) + x 2 cos (x lf x 2 ) , 

= 

x 2 = x : cos (x 2 , xj + x 2 cos(x 2 , x 2 ) . 


(1.16) 
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In three dimensions, we have 


x[ = x 1 cos (x{, x 3 ) + x 2 cos (xj, x 2 ) + x 3 cos (x' lr x 3 ) , 

• x 2 = x x cos (x 2 , x 3 ) + x 2 cos (x 2 , x 2 ) + x 3 cos (x 2 , x 3 ) , (1.17) 

[x 3 = x 3 cos (x 3 , x 3 ) + x 2 cos (x 3 , x 2 ) + x 3 cos (x 3 , x 3 ) . 

I wrote this out explicitly so that you would notice a pattern in the indices which indicates 

the above can be written in matrix language. Let us define 


X i:j = cos (x ^ ,Xj) . 


(1.18) 


If we arbitrarily choose to represent the position (x lf x 2 ,x 3 ) as a column matrix, the 
above relationship can be written, 


/ , \ 
X 1 

/ 

x 2 

= 

X 3 

\ y 

\ 


Xn X 12 X 13 
X 2 i X 22 X 23 
X 3i X 32 X 33 


\ 

/ \ 
x i 


X 2 

/ 

X 3 

V / 


or more abstractly 


(1-19) 


x ' = Xx, 


( 1 . 20 ) 


where it is understood x ' and x are column matrices and X is a 3x3 matrix. The easiest 
way to see this is just to do the matrix multiplication on the right-hand side to see that it 
reproduces the rotation. 


Written out in explicit index notation, the above relation may also be written as 

3 

x ’i=2^ X i (1.21) 

j=i 

where it is understood that the free index "i" takes on values 1,2,3. 

Question: What if we knew the x^ instead of the Xj ? In other words, what is the inverse 
relationship between these quantities? Using matrix notation we have 

x ' = X x 

=> X" 1 x' = X _1 Xx = x. 


so 


1 


x = X -1 x ' 
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or 


X 


i = 2 * 

j-1 


( 1 . 22 ) 


As stated above, there is a simple way of getting X (given X) for rotations. Consider: 


2*2 - 


x/ + xd + xd. 


By definition 

2 x '2 = 2 

= 2 X ij X ik x j x k . 

i,j,k 

Require: 2 x i 2 = 2 x i. 2 (kngth preserved) 

i i 

2 X ij x ik x j x k = 2 x i 2 * 


(1.23) 


=> 


(1.24) 


i,j,k 


The coefficient of each XjX k term (j ,k independent) must be the same on both sides: 
2 MCk = 0, j * k, 

i 

2 X ij X ik = = k - 


=> 


(1.25) 


More simply 

2 ^ij^ik = ^jk r 

i 

t 


(1.26) 


"Kronecker delta 1 


where 


i] 


[0, i ^ j , 
If 1 = j • 


(1.27) 


Notice that the above is not a matrix statement (why?), but that the right-hand side, 6j k ,, 
is just index notation for the identity matrix. However, it may be cast into matrix language 
with the help of the transpose: 
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i = hj r 

=> S X ji X ik = 8 jk or ^ = 1* (1.28) 

i 

This last statement establishes that 

A, T = X' 1 , (1.29) 

for rotation matrices. In fact, the above equation is sometimes taken to define such 
transformations. Actually, the types of coordinate transformations allowed by this equation 
are a bit more general than simply rotations, as we will see shortly. The name orthogonal 
transformations are usually used for real matrices which satisfy X, T = A. -1 , and therefore 
preserve the lengths of vectors. Using the above, the relationship between the X| and the 
x ± may now be cast as 

x ' = Ax <=> x = A t x ' , 
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in matrix language, or 


x 


i= 2 


j X D 


X i = 2 X ji x i ' 


in component language. It is important to note that in place of 

2 ^ij^ik = ^jk ' 
i 

we may also write 

2 ^ji^ki = ^jk * 
i 

Deriving one from the other will be a homework problem. 


(1-30) 


(1-31) 


How many independent elements are there? Notice that the equation ^ ^ij^ik = ^j k is 

symmetric under j x k, so that this system of equations actually represents 6, not 9 equations. 
(The number of independent elements of a real symmetric 3x3 matrix are 6.) This means 
the number of independent is 9 — 6 = 3. This makes sense from the physical point of 
view of rotations in three dimensions, which require three independent angles, in general. 


Let me make three additional points about the ^ j: 

1. x" = (X 2 ki)x = k 3 x. 


(1.32) 


A -3 is an orthogonal transformation if Ax and A 2 also are. (We can view this as a 
transformation Xx followed by a second transformation X 2 .) That is, the product of 
orthogonal matrices is also an orthogonal matrix. The proof of this statement is left as a 
problem. 


x ' = X x x, 

x " = X 2 x ' = k 2 k x x. 


(1.33) 


Do the rotations in the opposite order: 


x = k 2 x r 
x = k x x = k x k 2 x 


(1.34) 


In general k{k 2 * k 2 Xx, (“noncommutative”) and the order in which rotations are 
performed is important. 
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3. As stated above, we started out describing rotation, but the A, i j can also 
represent operations which are not rotations. For example, 


X, 


0 -1 

0 0 


0 

-1 

/ 


This satisfies X:[ nv X inv = 1 , but describes a spatial inversion. This is not an operation 
which can be physically carried out on objects in general. We will consider such orthogonal 
transformations only sparingly. 


All orthogonal transformations satisfy 
f+1, rotations 

detX = \ . . (1.35) 

[-1, inversions 

A. i j which have det X = - 1 are called proper orthogonal transformations, and those that have det 
X = 1 are called improper orthogonal transformations. In general, improper orthogonal 
transformations can be thought of as an inversion plus an arbitrary rotation since 

det (X^ nvers ^ on X rotation ) — det (X^ nvers ^ on ) (det X rotat ^ on ) , 


which means the combination is still improper. Proving that det ( X) = 11 for orthogonal 
transformations is left as a homework problem at the end of the chapter. 

Scalar and Vector Fields 

You probably have an intuitive feeling as to the meaning of a “scalar”. It is something which 
is unchanged under a change in coordinates. Under rotations, a scalar field behaves as 

Scalar: cj>' ( 2 C' ) = 4>(x), (1-36) 

where x i = ^ ^-±j x j • It is understood that the x and x 1 represent the same point 

j 

described in two coordinate systems. See the below: 
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There are also pseudoscalars, which behave like scalars under rotations, but acquire an extra 
minus sign under inversions. 


Pseudoscalar: cj>' ( x ' ) = ( detX ) (j) (x) . 


(1-37) 
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Although we’ve informally defined a vector as a quantity with both magnitude and direction, 
we need a quantitative definition. Here it is: 

Vector: A i (x')= ^ XijAj(x). (1-38) 

j 

The A ±j are elements of an orthogonal transformation. Obviously, the transformation law 
for the coordinates x i , Xi is the template for this definition. 

There are also pseudovectors, which transform as 

Pseudovector: A^x' )= (detA)^ AijAj(x). (1-39) 

j 

These quantities do not change direction under inversions, as vectors do. We will see an 
example of a pseudovector shortly. 

Vector Algebra and Scalar Differentiation 

Given two vectors, A and B, we may form either a scalar or a vector. Let’s study these 
two possibilities. 

A* B= ^ AiBi. (1-40) 



denoted defined 

It is also called the “dot product.” Here is the proof it is a scalar: 



Notice I changed the order of the sums in some intermediate expressions, which is always 
allowed for a finite number of terms. Another, more geometrical, definition of the scalar 
product follows from the diagrams below: 
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Same situation with a new coordinate system: 



Therefore 



= cos 0 


A ' • B ' , 

=> A • B = | A | | B | cos0. 

The vector product (or cross product) of A and B is defined as: 
( A X B ) ^ Eij^AjBk. 
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denoted defined 


Define the “permutation symbols”, £ijk : 


U 23 

U 31 

U 12 

£ 321 

= £ 213 = 

£ 132 

all 

other 

£ ijk 


(1-43) 


Important property, Ujk = £jki = Ekijf £ ijk = -£ ikj = -e jik • This is known as the 
cyclic property of the permutation symbols. Using the Eijkf a more formal definition of 
the determinant of a 3x3 matrix, A, can be written 

detA = 2 £ nCmAm A 2E A 3m . (1.44) 

n,£, m 

We will use this definition shortly. In order to see how A X B transforms, we will need 
the following identity for orthogonal transformations, 
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2 jk ^jC^-km ( det A) 2 ^nCm^-inf 

j,k n 

Tliis identity will be proven later (end of the chapter.) Let’s begin: 


(A'XB 1 ). 


= 2 

j/k 


£ijkAjB k 2 


jk 7 -j C A[ A km B n 


■ 2 


^ijk ^-jC ^-km A(B m- 


(1.45) 


(1.46) 


Now use (1.45) on the right-hand side of (1.46): 

—> (A X B — ( detA) 2 EnCm^-in A(B m — 

n,(,m 

(j)(k) (j)(k) 

\ 


(detA) 2 Ai n ( £ njkAjB k ) 

n, j,k 


notice I am 
indicating a change 
in the dummy indices! 


= (detA) J ^in ( AX B ) n . 

n 


Therefore we have that 


( A ' X B ' ) ± = ( det A) ^ A in ( AX B ) n . 

n 


(1.47) 


The extra factor of (det A) indicates that A X B actually transforms as a pseudovector (assuming 
A and B are themselves vectors.) 


As an aside, there are many useful identities we can form with the £ijk by multiplying and 

summing on indices. For example, let’s evaluate: 2 E ijk £ Cmk* Following the index jockey 

k 

rules, the only objects we can build out of Si j and the £ijk which have 4 free indices are: 
hijSfmf SimScj r hi£6j m . 


Therefore, we must have 

21 ^ijk^tmk = CiSijShn + f-kdimScj + CsditSjm, 
k 
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where the Ci f 2,3 are unknown constants. One immediately can see that Ci = ( 
Now multiply both sides of the above by 6ij and sum over i and j : 

=> 0 = C 2 Scm + C 3&Cm => C 2 = -C 3 . 

Now multiply by 6j m and sum over j and m:: 

— > ^ jk Ef.jk — — 2C 2 Si[ . 

jrk 

Then, taking a special case, say i = E = 1, we then see that 

2 £ ljk £ ljk “ 2, 
j/k 

so therefore C 2 =-l and 

^ Ojk £[ jk — 26 iC , 
j/k 

— ^ ^ijk ffimk — SiE^jm — Simd£j . 
k 

From this last expression, we easily see that 

^ jk jk — ^ ( Sufi j j — ) — 3X3 — 3 — 6. 

hj/k if j 

There is an alternate definition of the vector product as well. It says: 

AX B = | A | | B | | sin0 | n , 

|a| = magnitude of A 
where |b| = " " B 

0 = angle between A and B 

n = unit vector using the right hand rule 

It doesn’t matter in which direction 0 is measured. This angle is measured in 
defined by A and B if they aren’t colinear. This definition is illustrated below. 


Download free eBooks at bookboon.com 


. (Why?) 


(1-48) 

(1-49) 

(1-50) 
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the plane 
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A,B in the 1,2 plane; 
C is directed along 
the +z axis. 


0 measured in the 1,2 plane 


Let’s now reconnect to the useful concept of unit vectors in this context. The scalar products 
are summarized as 


e i * 


hij. 


(L52) 
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Now, using the second definition of the cross product, we find 

eiXe 2 = |sin90°|e 3 = e 3 . 

Similarly, 

e 3 x ei = e 2 , e 2 x e 3 = e 1 . 


Notice that order matters here since 
e 2 X e i = - e 3 , 

and similarly for the other non zero products. Also 

eiXei = e 2 xe 2 = e 3 xe 3 = 0 . 

We can summarize all properties of unit vectors under the cross product as 

e i x e j = ^ £ijk e kf 

k 

and with a cyclic order wheel for permutation symbol indices: 

1 

=> ^23 = 1 , 

^213 6"tc« 



( 1 - 53 ) 


I can now show the first and second definitions of the cross product are equivalent by 
resolving the components of A and B into unit vectors. From the second definition, 

AX B = (Aiei+A 2 e 2 +A 3 e 3 )X(Biei + B 2 e 2 + B 3 e 3 ) 

= (A 2 B 3 - A 3 B 2 ) e i+ ( A 3 Bi - A 3 B 3 ) e 2 + (A 3 B 2 - A 2 B 3 ) e 3 . 


On the other hand, from the first definition, 

( AxB ) 1 = ^ E ljk A j B k = ei 23 A 2 B 3 + ei 32 A 3 B 2 = (A 2 B 3 - A 3 B 2 ) , 

j,k 
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( A X B ) 2 - £231 A 3 B l + £213 A l B 3f “ (A 3 B 1 - A 1 B 3 ) , 
( AX B )3 = £312AiB2 + £ 321 ^ 263 , = (A 1 B 2 - A 2 B 1 ) . 


Comparing, we see the two definitions agree. By forming the determinant and comparing 
with the above definitions for the cross product, one may also show that 


AX B = 


©20203 

A 1 A 2 A 3 , 
B i B P 3 


where the right hand side is understood as a symbolic determinant. 


From either definition of the cross product we have that 
A x B = -B x A , 


(1.54) 


(1.55) 


as well as 

A x (b x c) = B(A • c) - C(A • §) . (1.56) 

The last statement follows from Eq.(1.49) above. In addition, there is the cyclic rule, 

A • (§ x 6) = B ■ (c x a) = C • (a x b) , (1.57) 

which is a direct result of the cylclic nature of the permutation symbols. Given that A , B , 
and C are vectors, the quantity A • ( B X C ) is an example of a pseudoscalar. (This is part 
of a problem.) 

So, we know how to add and multiply vectors together. Question: does it make sense to 
define division? Consider 


AXB= C. 


Given that A and C and are known quantities, this equation clearly does not uniquely 
define B since any component of B along A will not contribute to the cross product. So 
the answer is in general, no. However, if we know both A • B and A X B , we can solve 
for B , say. (We will have a homework problem along these lines.) 

Differentiation of vectors with respect to scalars leads to new vectors. 


Download free eBooks at bookboon.com 


27 




MODERN INTRODUCTORY MECHANICS 


MATHEMATICAL REVIEW 


Given: 


Then 


Al = 

2 S 

j 

(a 

scalar). 


dA'. 

„ dA, 


dA, 


1 

ds' 

= y x ±j -4 

f 3 ds 

ll 

ds 



=> dAi = y 

ds' f 

Xij 

{ dA A 

\ ds j 

(1-58) 


A fundamental postulate of mechanics is that the time parameter transforms as a scalar quantity. 
Given this, these derivative transformations imply that velocity and acceleration are vectors. 


Alternate Coordinate Systems 

The radius vector, r , is easy to characterize in the three most common coordinate 
j. dr 

systems. We have ( r =- ; tis considered a scalar ): 

dt 
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Rectangular: r = Xi<ii + x 2 e 2 + x 3 e 3 (all ei = 0) 

Cylindrical: f = pe p + ze z (only e z = 0) 

Additional unit vector: 

Spherical: J = r e r 

Additional unit vectors: e^, eg 
v = r, a = v= r. 

As an exercize, let’s work out v and a in cylindrical coordinates. Follow a particle’s trajectory 
at two closely spaced moments in time: 



Displace the unit vectors to the origin: 


d<J) > 0 here 

2 



Can now see that (de p = e p - e^) 
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I dip | = | d(t> | , 

de p = direction X length, 
— e (j) d(|) • 


Likewise (de^ - e| - e| ) 

I di(j) | = | d<j) | , 
de^ = -e p d(j). 


Therefore 


dep 

dt 


d4> - a •- 
—e^ or e p = 4>e^ 


and 


de^ 

d4> 


dcj) ^ x 

TT e P or e * 



Thus in cylindrical coordinates, 


v = pe p + pe p + ze z , 

=> v = pe p + p 4>e ^ + ze z . 

a = pe p + pe p + pfjte^ + p(j)e 0 + pfjte^ + ze z , 
=> a = (p - p(j) 2 ) e p + (2p<j> + p<|)) e^ + ze z 


Since the unit vectors in rectangular coordinates do not change, it is 
the decomposition of the unit vectors in other coordinate systems in 
coordinates. 


Cylindrical: 


e P 

e 4> 

©z 


(cos 4>, sin 4), 0), 

(- sin <j), cos (]), 0), 

( 0 , 0 , 1 ). 


Spherical: 


Also 


e r = (sin 0 cos 4h sin 0 sin 4>, cos 0), 

- e e = (cos 0 cos 4h cos 0 sin 4>, - sin 0), 
e^ = ( - sin 4>, cos 4v 0). 
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P 



cyclically in spherical coordinates. 
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Angular Velocity 

Another important concept for particle motion is angular velocity. Here we will rely mainly 
on intuition to understand the concept. First, identify an instantaneous circular path: 


particle 

trajectory 


circle tangent Instantaneous 



Flip this picture on it’s side. Choose a cylindrical coordinate system with the origin along 
the axis of rotation. Use this origin to define the radius vector, r. 



Let 0 ) be directed co-linear to the axis of rotation, which is the third axis above. Angular 
velocity and velocity are related by the equation 

v = to x r . (1.65) 

Notice, using a cylindrical coordinate system, 

w = co z e z , v = p<j>e + , r = pe p + ze z , 

=>v = co x r => pcjce^ = co z e z x (pe p + ze z ) = co z pe r 
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This means 


co z = (j> , (1.66) 

so that 0 ) z has the meaning of instantaneous time change in the angular position. Notice 
the direction of w is given by a right hand rule. (Curl the fingers of you right hand in the 
direction of and your thumb points in the direction of c5). A common case is uniform 
circular motion: 



co is a constant 
vector pointed along 
the 3-axis. 


v = 

a = -P<i> 2 e p , 
q = <j>e z . 


Differential Operators and Leibnitz Rule 

A frequently occurring mathematical operation is the gradient: 


V - 



d 


dx i 


(1-67) 


The x,y,z components of this are just the usual partial derivative operators. When 
operating on a scalar, it gives a vector: 


v 41 = ? e * ax : ° r ( v ' 4 

>■) - 

^ dx ± 

(1.68) 

Proof: The chain rule says, 



df v dx J 's 

atj) 9xj 


dxj A? dx^ dxj A 

r 1 dx^ dxj 



Now remember 



Ski 


(1-69) 
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Therefore 


d<t>' = V ^ dfo 

dx^ y 1 -’ ^Xj 




(1.70) 


which states that V(|) transforms as a vector. A physical interpretation of the gradient operator 
acting on a scalar is given by the following sketch. Think of lines of constant density in 2 
dimensions again. Construct the 1’ axis tangent to the constant 4> line, at x ' : 



Clearly, 
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ad) - ad> ~, 

—- = 0 => V<t> = —7 e 2 only. 

dx-L dx 2 

We conclude that the gradient at each point is a vector pointing perpendicular to lines of 
constant (|) at that point and points in the direction that (|) is increasing. 

In three dimensions we can orient a new coordinate system so that the 1’2’ plane is now 
locally tangent to the {[> = constant surface at x' : 



Again 

d(b - , 

V <t> = —— e 3 only. 

dx 3 

So the gradient is perpendicular to lines (2 dimensions) or surfaces (3 dimensions) of 
constant 4>, and points in the direction (j) is increasing. When dotted into a vector field, 
the divergence gives a scalar. (You should be able to prove this statement.) Thus, V ■ A 
transforms as a scalar if is a vector. 


A further mathematical operation, the curl of a vector, 

Curl: (Vx a)- ^ £ ijkVjA k , (1.71) 

j/k 



denoted defined 

can also be defined using the gradient V if but will not be used extensively in this course. 
If A is a vector, V x A can be shown to transform as a pseudovector. 


Another important operator is the Laplacian: 

v ■ v = y -^2 - v2 • 

4* ax? 


(1.72) 
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One can show it produces a scalar or vector when acting on a scalar or vector, respectively. 
Here is the proof it transforms as a scalar when acting on a scalar: 

v-7' jx' V 1 V 5 d 

v * ~ ~ Z- — *• 


i ^ x i i ^ X i dx i 


Since (|> = 4>, we have 


v 2 f - 2 


3Xj\ d ( dx } 


i^ k V ax i/ dx j \ dx 'i / dx f 


4 > • 


As shown above 


<9x-i 


dx i 


T - ^ij. 


j* V i j 


d d , 
ik f .. 4* r 




The proof that V 2 operating on a vector yields another vector is similar. 


Finally, there is a useful rule for taking the derivative with respect to a variable which is 
contained in the upper and lower limits of an integral. This result, called the Leibnitz rule 
for differentiation, states that 


d 

dt 


b(t) _ b(t) 

J dx f(x, t) = f(b(t) , t)b(t) - f(a(t) , t)a(t) + J dx 


a(t) 


a(t) 


df(x, t) 

at 


(1.73) 


When the differentiating variable appears in the upper and lower limits of an integral as 
well as in the integrand, one obtains three terms. The first term is the value of the integrand 
evaluated at the upper limit, f (b(t ) , t ), times the derivative of the upper limit, b(t) . 
The second term is similar but represents minus the integrand at the lower limit, a(t) . 
times the derivative of the lower limit. The third term represents the contribution from 
the t-dependence in the integral itself. Eq.(1.73) is the simplest case of the dependence 
on a single variable, t; we will use a more general version of the Leibnitz rule in a three 
dimensional context in the next Chapter. 


Complex Variables 

We will have occasion to use complex numbers and variables in our study of mechanics. 
The imaginary number i is given by 

i= -J— 1 r 
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and when raised to various powers gives, 


= -1, i 3 = -i, i 4 = 1. 


These numbers can be arranged on a unit circle in a plane: 
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This is called the complex number plane. As the figure suggests, any complex number can 
be written as a combination of a real and an imaginary number, 


z = x + iy. 


(1-74) 


where x and y are both real numbers. The complex conjugate of z is given by x - iy 
and is denoted as z . These two numbers specify a location in the above plane on the 
real (x) and imaginary (y) axes. The real and imaginary parts of z are separated off with 

Re(z) = x, Im(z) = y, 
and the distance from the origin to the point z is given as 

| z | = -Jz z = *Jx 2 + y 2 . (1-75) 

From plane trigonometry, any complex number on the above unit circle may be written as 
z = cos4> + i sin4>, (1.76) 

where the angle 4* is measured from the real axis, x. It was Euler who showed this may 
be written as 


z = e 


(1-77) 


for . Given this result we can now represent any number in the complex plane as 


z = z e 


(1-78) 


We often encounter ratios of complex numbers, 


z = 


We may always factorize such a ratio into real and imaginary parts by multiplying top and 
bottom by the complex conjugate of z 2 : 


z = 


Z 1 z 2 


12 ' 


which may also be written as 
, i ( < t ) l - <t>2) 


Z = 
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Complex numbers are extremely useful in solving linear differential equations, which is 
where we will see them next in this course. 

CHAPTER 1 PROBLEMS 

1. Let 


/ 



\ 


/ 



\ 


0 

1 

3 



6 

2 

-2 


-1 

2 

0 

, B = 


1 

0 

4 


2 

4 

0 



4 

1 

-1 

\ 



) 


V 



/ 


a) Find the matrix product AB. (A symbolic manipulator will be handy) 

b) Find det (AB), and show that det (AB) = (det A) (det B) 

c) Form the matrix product B T A T and show that B T A T = ( AB ) T , as it 
should. 

2. Try proving that 

A" 1 A=1, 
implies that 

AA _1 =1. 

3. Prove that the inverse of the transpose is the transpose of the inverse of a square 
matrix (assuming it’s inverse exists!). 

4. Given 2/ ^ij^ik = ^jk ' show that 

i 

2/ ^ji^ki = ^jk r 
i 

follows immmediately. [Hint: Use matrix notation and remember that X _i = X T for 
these matrices.] 

5. Starting with (X is a matrix; T denotes transpose) 

X- 1 = X T , 


prove that 
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det X = ±1 . 

[Hint: X = det X T .] 

6. Given that the 3x3 matrices Xi and X 2 represent orthogonal transformations 
on vectors x and x ' 



(x'= XiX, x" = X 2 x') 


show that the matrix product, X 2 X 3 ,, is also an orthogonal transformation. 

7. Using the definition given in the text, 

det A = ^ £ ijk A iiA 2 jA 3k , 

i,j,k 

show that the determinant of a 3 x 3 matrix may also be written in the forms: 
a) det A = ^ £ ijk A iiAj 2 A k3 , 



> Apply now 
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b) det A ^ E ijkGmn^a^mj^nk * 

8. Argue that the determinant of an antisymmetric 3x3 matrix, 

A T = - A , 

is zero. Is the determinant of an antisymmetric 4x4 matrix necessarily zero? 

9. We need to complete the proof of the relation we found necessary when proving 
that A X B actually transformed as a (pseudo) vector. This proof was based upon 
the supposed identity, 

^ £ ijk ^j^km = £ nta^in * 

j,k n 

By multiplying both sides by X ip (and summing on i) and choosing a special case, show 
that this identity may be reduced to Eq.(1.44) of the text (almost; see prob. 7), implying 

det A = C. 


10. (a) Show Eq.(1.56) follows from Eq.(1.49). 
(b) Show the cyclic property Eq.(1.57). 

11. Show that 


B x (a x b) = (a x b) . 

12. Given that A , B and C are vectors, argue that A • (b x c) transforms as a 
pseudoscalar. 

13. Does the angular velocity of a particle, u), transform as a vector or 
pseudovector? Explain. 


14. Evaluate the sum: 

LjkGntfTnjm * 




k, j, m 


[Hint: Use results already proven. 


Download free eBooks at bookboon.com 


41 




MODERN INTRODUCTORY MECHANICS 


MATHEMATICAL REVIEW 


15. (Adapted from Marion and Thorton.) Let’s investigate a point raised in the text. 
Consider 

A x B = C , 

where A and C are considered known quantities. Clearly, this equation does not uniquely 
define B since any component of B along A will not contribute to the cross product. 
However, consider 

A • B = S , 

where S is a known scalar. Show that we can solve for B and that ( A considered nonzero) 

B = -J- (SA + C x A) . 

A 2 v ’ 

16. Express the spherical unit vector, e e , in terms of the cylindrical unit vectors 
e p, e ^ and e z • The x, y, z components of these vectors are 

e 9 = (cos0cos4>, cos0 sin((), -sin0) , 

and 

e p = (cos(j), simj), 0 ) , 

= ( -sin4>, cos4>, 0 ) , 
e z = (0,0,1). 

17. Show that the velocity vector in spherical coordinates is given by: 
v = re r + r0e e + r sin0 (j) e^. 

[Hint: Use the decomposition of e r found in the text.] 

18. In spherical coordinates, show that 

a) e e = -0e e + (|>cos0e r , 

b) = -(|)cos0 e 0 - 4>sin0 e r . 

[Hint: See decomposition of ee, e^ in the text.] 
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2 

19. A particle follows the trajectory y = Kx 
S, in the xy plane as shown. 


(K is a constant) with constant speed, 



o 


a) Show that x • x =0 anywhere on the trajectory. 

b) Find the acceleration, x, when the particle is at the origin, 0. 

c) Find the instantaneous angular velocity, (O, at 0. 

20. Given the vector function v = 2 0xyi+ 2 5yzj+ 15xyk , find: 
(a) V • v (b) V x v 
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21. Given that 



and 


A : = 2 ' 


k 


prove that: 


a ) V 1 • A 1 transforms as a scalar, 


b) V 1 x A 1 transforms as a pseudovector. 

22. Show that (V 2 = ^ V 2 ) 

j 

V 2 (f g) = gV 2 f + fV 2 g + 2Vf*Vg. 

23. Show that 

z = cos([» + i sin4>, 
may be written as 
z = . 

[Hint: Consider the power series expansion of the exponential.] 

24.Show that the angle (j) in z = e 1 ^ .is given by 0 < (j> < 2 jt 
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2 NEWTONIAN MECHANICS 


REVIEW OF NEWTON'S LAWS 

Newton’s 3 laws: 

1) A body remains at rest or in uniform motion unless acted upon by a force. 

2) A body acted upon by a force moves in such a manner that the time rate of 
change of momentum equals force. 

3) If two bodies exert forces on each other, their forces are equal in magnitude and 
opposite in direction. 

We will see that these dynamical laws lead to certain conservation theorems. In fact, it is 

now known that the conservation laws have a greater range of validity that Newton’s laws, 

so we can regard these as a stairstep toward a more fundamental outlook. More sophisticated 
formulations of mechanics do not even use the concept of “force”. So Newton’s laws move 
us toward a more fundamental outlook in both physical content as well as formalism. 

Essentially, Newton’s laws are a framework applicable to macroscopic motions due to 
any kind of physical interaction as long as the characteristic velocities involved are small 
compared to the speed of light and characteristic length scales are large compare to atomic 
dimensions. By a “physical interaction” I mean the underlying force of nature responsible 
for producing the motion of objects. For example, electrodynamics often manifests itself 
in our macroscopic world in a sort of disguised form in the form of elastic and inelastic 
(frictional) forces. By elastic, in this context, I mean that no energy is lost in the form of 
heat during the interaction. 

Let me paraphrase Newton’s laws, as applied to point objects, to help you remember them 
better and to be a bit more quantitative. 

I. Let x(t) be the position of a particle at time t and let F(t) be the force on it. 
Then Newton’s first law says 

x(t) = vt + x 0 if F(t) = 0, (2.1) 

Of course x(t) = vt + x 0 is only going to be true if particle position is measured in a 
non-accelerating or so-called inertial reference frame. So, we need an initial inertial reference 
frame to define another one. The inertial frame concept is clearly an idealization, only true 
to a certain extent in applications here on Earth. (We already are idealizing to point particles 
as well!) However, this will not prevent us from an intuitive and fruitful application of 
Newton’s laws in real-world situations. 
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II. Let the momentum of a point particle be defined as 

p s mv, (2.2) 

where m is the particle’s mass and v is it’s velocity. Imagine summing all the individual 
forces to get the total force on the particle, 


- 2 L 


Then Newton’s second law says 

dp d , _. 

F = — = — (mv) . 
dt dt 

This defines force (undefined in Newton’s first law) as long as mass is also defined. 


(2.3) 
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III. I’ll use the following picture to help explain Newton’s third law: 






F 21 - 


F 21 and F 12 act on different objects! Note that the forces between particles do not necessarily 
point along the line connecting them, although it usually does. 

I invoked the concept of mass above. How is this defined? It can be regarded simply as the 
proportionality constant between force and acceleration in some system of units: 

F = ma. (2.4) 

However, this viewpoint is a bit illogical, because we are now using this one equation to 
define both force and mass. If we are going to simply regard mass as a proportionality factor, 
then we should have a means of defining both force and acceleration. We can always use 
a gravitational field to independently measure such quantities in the form of the equation, 

W = m g g, 

where W is the weight (measured on a spring scale for example) and 9 is the acceleration 
due to gravity. "m g " is then the mass as measured using gravity. It is now possible to imagine 
accelerating the particle using the same spring. If we displace the spring the same amount, 
we are guaranteed to exert the same force. If we allow the particle to be accelerated, we 
will now have that 

W = m a a, 

where W is the exact same force as before. Although the same force is being used in both 
situations, it is possible that the accelerations a (inertial) and g (gravitational) will not be 
the same and therefore the proportionality constants, mg and ma will not be equal. Let’s 
call ma the inertial mass and mg the gravitational mass. The question 


Download free eBooks at bookboon.com 


47 




MODERN INTRODUCTORY MECHANICS 


NEWTONIAN MECHANICS 


m g = m a 

is an experimental one and has been tested extremely accurately; one such test will be 
described in Ch.5. As far as is known, however, the equality m a = m g is exact. We will 
use "m" to denote the mass in all future situations. 

Let us investigate the solution of Newton’s force law for a point particle in two fairly 
common circumstances. The general form of the force law in one space dimension is just 

mx = F (x, x, t) . (2.5) 


x 


dx 

dt 


x 



What I am suggesting in the notation F (x, x, t) is that the force on the particle can be a 
function of it’s position, x ( t ) , it’s instantaneous velocity, x(t) , as well as having explicit 
time dependence. Without knowing the specific form of F (x, x, t) it is not possible to 
solve Eq.(2.5) in general. However, let’s say F (x, x, t) F(t) only. Then the first 
integral of Eq.(2.5) gives 

x(t) = k 0 + - f dt 1 F(t' ) , (2.6) 

m i 


from which the velocity at time t can be found. In the above, x 0 is just the constant of 
integration, which in this case just represents the initial velocity. (Just set t=t 0 above to 
see this.) If our goal is to find x (t) , we need only integrate again: 

ft t ■ 

x(t) = x 0 + x 0 (t - 1 0 ) + — f dt f dt ” F(t ” ) . (2.7) 

m . . 

Again, a constant of integration, x 0 , has emerged. Clearly, since Eq.(2.5) represents a 
second-order differential equation, we will always have two such constants to supply. In 
this case, we are simply requiring 


x(t 0 ) = x 0 , (2.8) 

x (t 0 ) = x 0 . (2.9) 

Eqs. (2.8) and (2.9) are often called initial conditions. Usually, one chooses t 0 = 0 . 
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A common application of Eq.(2.5) is when F(t) =ma, where "a" is a constant in time. 
(This is the case of gravitational attraction of a point mass near the Earth’s surface, for 
example.) Then, we obtain 

t t ' t -i 

/ dt' / dt" = / dt' (t t 0 ) = - (t - t 0 ) 2 

to t 0 to 2 

which gives 

x(t) = x 0 + x 0 (t - t 0 ) + | (t - 1 0 ) 2 , (2.10) 

which should be familiar to you as the equation for instantaneous position under constant 
acceleration. 

Another common case is where F (x, x, t j —> F( x) only. Then, we can solve for the 
motion again, as follows. Integrate both sides of Eq. (2.5) over x. Since 

x(t) t 

/ dxx = f dt ' x (t ') x (t ') ( 2 . 11 ) 

x n t 0 



Brain power 


% 


By 2020, wind could provide one-tenth of our planet's 
electricity needs. Already today, SKF's innovative know¬ 
how is crucial to running a large proportion of the 
world's wind turbines. 

Up to 25 t of the generating costs relate to mainte¬ 
nance. These can be reduced dramatically thanks to our 
^sterns for on-line condition monitoring and automatic 
lunhtation. We help make it more economical to create 
clearfflfcdneaper energy out of thin air. 

By snaring our experience, expertise, and creativity, 
industries can boost performance beyond expectations. 

. Therefore we need the best employees who can 
meet this challenge! 

The Power of Knowledge Engineering 


Plug into The Power of Knowledge Engineering. 
Visit us at www.skf.com/knowledge 


Download free eBooks at bookboon.com 




Click on the ad to read more 


49 



MODERN INTRODUCTORY MECHANICS 


NEWTONIAN MECHANICS 


1 

2 




( 2 . 12 ) 


we have 


1 x<t> 

- m (x 2 (t) - x„) = f dx'F(x'). (2.13) 

One may confirm that Eq. (2.13) is the first integral of (2.5) by taking the derivative of 
(2.13) with respect to t. 2 In this case x 0 is just a constant. Also, Eq. (2.13) is actually just 
a statement about energy conservation. One recognizes the difference in kinetic energies 
on the left hand side and the right hand side is just the work done. This will be discussed 
in more depth later. 


It is possible to integrate (2.13) again, as follows. First, solve for x (t) , 


x 


2 x(t) 

(*) - *J- I dx ' F(*') + *;. 


(2.15) 


Then using the general identity (which also applies in more than one dimension if the 
integration is along the particle’s path), 


x(t) 


/- = Jdt, 

J x i 


yields 


/ ax 

-j— - - (t - to). (2.16) 

- f dx‘ F(x") + xj 

The appropriate sign in Eqs.(2.15) and (2.16) must be chosen for a given problem. Another 
constant of integration has been supplied on the right hand side of Eq.(2.16), which in 
this case is just the initial time, Notice that while Eq.(2.16) is a formal solution to 
the problem, it gives t (x) rather than the more usual x (t) . However, if the relation of 
x(t) and t is one-to-one (we can imagine breaking the overall motion of the particle up 
into submotions for which this is true), then one can in principle invert the relation t(x) 
to get x (t) • 


Let us again take an example to illustrate the use of Eq. (2.16). Consider a stretched spring 
attached to a vertical support, as shown. (Imagine the spring is nearly massless.) 
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/ / / / / 


m 


////// 


at rest 

(a) 


F(x) 


m 


origin 


+x 


stetched 

(b) 


Let us pull the mass m down a distance x 0 and then release it from rest. Let us assume that 
(“k" is a constant called the “spring constant”) 

F(x) = -kx , k>0, 

for small extensions of the spring, a force law for springs known as Hooke’s law. (We’ll 
examine this force law in the next Chapter.) The goal will be to use Eq.(2.16) to find x(t) 
for an appropriate portion of the motion. 


Eq.(2.16) reads in our case (choosing x 0 = 0, t 0 = 0) 

dx' 


x(t) 

7 


, 2k r 

I -f x' ' dx' 


= t. 


(2.17) 


m 


where we have chosen the negative sign in Eq.(2.16) for our coordinate system. We now 
do the integral: 

x(t) 


r~ x ^, 

-Mi 


dx' 


V 


x 0 - X 


■ -M (“--if) - “'-*4 


■JfMfH- 


(2.18) 


Note that the inverse sine is not defined until it’s range is specified. I have chosen the 
principal branch above, symbolized by writing "Sin -1 ". (Any other branch would do.) 
Putting (2.17) and (2.18) together now allows us to solve for x(t) as 


x(t) = x 0 sin 


( 


-A— t + — 
m 2 




(2.19) 
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This type of motion is termed “simple harmonic”. We have now solved for the motion of 
the point mass along the portion of the motion where the connection between t and x is 
one to one. Obviously, since the motion is simple harmonic we can get x (t) at any time 
t by just writing the usual “sin” function above, as we have already done. Thus 


x(t) = x 0 cos 



( 2 . 20 ) 


The one to one submotions of the system are connected by “turning points”, where the 
motion of the particle reverses, occuring at x = ±x 0 in this example. We will discuss 
turning points more extensively in the material to follow. The period of motion of this 
system is given by the amount of time necessary to complete one spatial cycle of motion 
and is given by the amount of time necessary for the arument of the cosine in (2.20) to 


increase by 2it, namely 1^p r i n g 
amplitude of motion! 



Note that this period is independent of the 


Simple Examples Using Newton’s Laws 

We will look at four simple examples as further applications of the ideas in Newton’s laws. 
These should be considered as purely review. 
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Example 1 

We illustrate one-dimensional harmonic oscillation in the form of an oscillating spring, 
but now let us imagine a periodic force is applied to the system, initially at rest. We take 

mx = F (x, t) , (E.l) 

F (x, t) = - kx + F 0 sin oat. (E.2) 


The equation of motion can be put into the form 
k F n 

x + — x = — sin cot. 

m m 

Let us solve (E.3) subject to the initial conditions, 
x(0) = 0, x(0) = 0. 


(E.3) 


(E.4) 


Solution: 

This problem was chosen because it illustrates some techniques of solving linear differential 

k . 

equations. Let us call oo 2 = — . Then the complimentary equation is 

J m 


x c + « 0 x c = 0, 


(E.5) 


the general solution of which is 


x = A sin oo„t + B cos oa n t. 

c 0 0 


(E.6) 


“ O ) 0 " is now identified as the angular frequency of motion, co 0 = 2jtf 0 , where fO is 
frequency. The particular solution can be found by the method of undetermined coefficients. 
The trial solution is (assume to ^ oa 0 initially) 


x = C sin cot + D cos cot. (E.7) 

Substituting in (E.3) now gives 


sin cot 

• c — 

F 

■ L o 

(E.8) 


m(oo 2 - co 2 ) 

cos cot 

: D = 

0. 

(E.9) 

Thus, the general 

solution 

is 

F n 


x = x c 

+ x r = 

-^- — sin cot + A sin oo 0 t + B cos oo 0 t. 

(E.10) 


m(o) 0 - 0) ) 
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It only remains to satisfy the initial conditions. We have 


x(0) = 0 


B = 0, 


(E.ll) 


x(0) = 0 


A = - 


F 0 w 


mo) 0 ^o)g - ca 2 j ' 


x = 


F n [ . 0) i 

- — 11 -— sin cot-sin ca 0 t . 

mfoDg - to ) l to 0 / 


(E.12) 

(E.13) 


Let’s now examine the solution for to = to 0 . We have the same complimentary solution, 
(E.6), as before, but now the trial particular solution must be taken as 


x p = Ct sin to 0 t + Dt cos to 0 t. 


(E.14) 


Again, substituting in (E.3) gives (the coefficients of the t sin to 0 t and t cos to 0 t 
terms vanish) 


sin to 0 t : D = - --, 


2mto„ 


cos out : C = 0 . 


(E.15) 

(E.16) 


The general solution is thus 


x = A sin to 0 t + B cos <o 0 t - 


2mto 0 


t cos <o 0 t. 


(E.17) 


Again, supplying boundary conditions, 


II 

o 

w 

ii 

o. 

(E.18) 

x(0) = 0 => A 

F o 

2moOg 

(E.19) 

Tit is gives finally 



F 

x = — 5— sin oo n t 

2mff»o 

F 

--— t cos O) 0 t. 

2mo) 0 

(E.20) 


Notice that as t —> > the motion in (E.20) becomes unbounded. We will later study 
in detail a more realistic case with damping where the unbounded motion for 0 ) = O) 0 is 
removed and instead results in a phenomenon known as “resonance”. 


The following is a simple one-dimensional dynamics problem. 
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Example 2 

Two sprinters are running a 200m race. Both runners have the same top speed, lOm/s, at the 
end of their initial acceleration, but Michael reaches top speed after only 15m while Frankie 
takes 20m to reach the same top speed. Let’s idealize the acceleration of both sprinters as 
a constant. By how many seconds does Michael beat Frankie in running the race? How far 
ahead is Michael from Frankie when he crosses the finish line? 

Solution: 

The acceleration profile of the two runners is as shown. 

a 

- 1 ^-Michael 

I 

I 

. . ^-Frankie 


; ; a=0 for both j 

15m 20m 200m 

The velocity and position of the runners during the acceleration period are given by Eq. 
(2.10) and it’s first derivative (x 0 = x 0 = t 0 = 0) , 


v(t) = at. 

(E.21) 

at 2 

x(t). 2 . 

(E.22) 

Notice that we do not know either the acceleration or the time it takes to reach top speed, but 
Eq.(E.21) and (E.22) give us two equations in two unknowns. Plugging in v(t) = 10m/s 
in (E.21) and either x(t) = 15m or 2 0m in (E.22) then gives 

= _ m/ 

“michael ^ /S 2 ' 

(E.23) 

a ^ m/ 

“■Frankie ^ /S 2 " 

(E.24) 

Using these results, solve for the acceleration time 

from (E.21): 

Michael 3 s, 

(E.25) 

(t, ) . = 4s. 

\ a./F ran k ie 

(E.26) 
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We now simply solve for the time necessary to run the remaining distance at the constant 
speed of lOm/s from (different x(t) than (E.22)) 

x(t) = x 0 t. (E .2 7) 

For Michael x(t) = 185m, while x(t) = 180m, for Frankie; x 0 = 10 ■ Calling 

/ s 

this time t 0 , we find 

fioL.- .-i - 18 - 5s - < E - 2 «> 

(tolprank = 18 ' ° S ' < E - 29 > 

Adding t a and t 0 , we then obtain that Michael beats Frankie by only .5s. In addition, 
when Michael crosses the finish line, Frankie will be 

10 % x .5s = 5m 

behind him. (There are much simplier solutions of this problem!) 

Of course, in more than one dimension, Newton’s law of motion for point objects becomes 
a vector law, 

mx = F (x, x, t) . (2.21) 

If we assume two dimensional motion and use Cartesian coordinates, this becomes 


mx = F x (x, y, x, y, t) , 

(2.22) 

my = F y (x, y, x, y, t) . 

(2.23) 

An important special case occurs when 


F x (x, y, x, y, t) -»> F x (x, x, t) and F y (x, 

y, x, y, t) - F y (y, y, t) , 

because then the two equations can be solved independently, although the initial conditions 
still tie the two motions together. This is the case for the two dimensional harmonic oscillator 
(next Chapter) as well as for gravitational force near the Earth’s surface. In the latter case 
we have simply (x is horizontal, y is vertical) 

x : mx = 0 , 

(2.24) 
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y : my = -mg, (2.25) 

where g is the acceleration due to gravity. (We will use the approximate value 9 • 8m/s 2 for 
g .) The solution of these equations are of course quite easy, but orchestrating the solutions 
to reach a desired set of end conditions can still be challenging. 

We will consider two examples of two dimensional problems, one static and one dynamic. 


Example 3 


/ / / 



A rope is attached to the ceiling and a force F is applied to it, as shown. The rope will 
break when the tension in any portion is greater than a certain magnitude, T '. Find the 
downward force, F, which breaks the rope. Which side, in general, will break? 
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Solution: 

The force diagram is as below. 



We have 

i 

E. (30) 
E. (31) 


x : Tj cos 0j - T 2 cos 0 2 = 0, 
y : -F + Tj sin 0 : + T 2 sin 0 2 = 0. 


Since in the diagram 0!>0 2 , Eq.(E.30) means that TT>T 2 , so it is the left hand side of 
the rope which will break first. Given this, use (E.30) to eliminate T 2 from (E.31): 

F = T x (sin 0 : + tan 0 2 cos 0-J. (E.32) 

Setting Tr = T ' and using specific values for 0 12 in (E.32) gives the force which breaks 
the rope. 


Example 4 

The rim of a basketball goal at a carnival is designed so it will only allow basketballs to go 
through the net when the angle of entrance is greater than a , as shown. (We are clearly 
making idealizations in the geometry and the ball’s motion.) 
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The net is at a height, H, above the ground, and the player is located a distance D from the 
basket. The ball is released a distance h from the ground at an angle p, again as shown. 



With what angles, p, can the ball be released and still make it directly through the basket? 


Solution: 

We will locate the origin of the x,y coordinate system at the point of release of the ball. 
The instantaneous x and y positions at time t are then (v 0 > 0 ) 


x(t) 

= (V 0 cosp)t, 

1 2 

(E.33) 

y(t) 

- (v 0 sinp)t - — gt 

(E.34) 


If we require x - D and y = H-h, we will have two equations in two unknowns, b and t 
(the time to reach the hoop), 


D 

v 0 cos (3 

H-h . gD 

- = tan P-;-;— . 

D 2v o cos p 


(E.35) 

(E.36) 


The angle at which the ball intersects the hoop is now 
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y n { d u 

tan a. = -- - tan p + —5 - 5 — . 

v x v 0 cos (3 

The condition for making a hoop is clearly 
tan a > tan a 0 => a > a 0 . 


We can use the trigonometric identity, 

1 


cos (3 = 


1 + tan 2 (3 


to put this in the form 


2 2 

9 v n v n 

tan 2 (3 - — tan (3 + 1 - — tan a 0 > 0 


gD 


gD 


2 

(E.36) may now be used to eliminate v 0 from (E.40), with the resultw 




1 + tan 2 (3 


tan (3 - 


( H - h ) 


D 


tan 3 (3 - 2 a tan 2 (3 + tan (3 - 2a > 0, 


(H - h) 1 

a ~ --- + — tan a n 

D 2 


We get a cubic condition in tan P as our general result. 


Let’s put some numbers in (E.41) and (E.42) to get a feeling for the condition. Let’s 
(H - h) 1 

- = - , a 0 = 45 . 

D 2 

=> tan 3 p - 2 tan 2 p + tan p - 2 > 0 

There is one real root of the cubic equation tan 3 P 0 - 2 tan 2 P 0 + tan P 0 - 2 : 
is given by 


tan P 0 = 2 => p 0 - 63.4° 


Thus, P > P 0 is required to make a basket in this case. (One must also have the 
v 0 as well, given by (E.41) above.) 
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(E.45) 

> 0 It 

(E.46) 

correct 
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Single Particle Conservation Theorems 

As I said at the beginning of this Chapter, the conservation laws we will recover from 
Newton’s three laws constitute more general results than Newton’s laws themselves. We will 
now consider conservation laws for a single particle; multi-particle conservation laws will 
be considered later. 

I. The total linear momentum, p, of a particle is conserved when the total force, F, 
on it vanishes. 

Paraphrase of I: F = 0 => p = constant vector.. 

Proof: The proof follows immediately from Newton’s second law, 


dp 

r - ft' 

(2.26) 

F = 0 => P(t) = p 0 . 

(2.27) 

where p 0 is a constant vector. 


Angular momentum, L, and torque, N, for a single body are given by 


L = r x p. 

(2.28) 

N = r x F. 

(2.29) 


Clearly, both L and N are origin-dependent quantities, unlike linear momentum, P. For 
example, consider ( p , F in a plane, say) 



Here, r x p is “up” (along 3) and r x F is “down.” Relocate the origin, 0 : 
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Now r x p is “down” while r x F is “up.” Thus, it is important to maintain a consistent 
origin when using these concepts. 

II. The total angular momentum of a particle, L, subject to zero total torque, N, is 
conserved. 

Paraphrase: N = 0 => L = constant vector . 
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Proof: We have 

£ = A( f xp).fxp + rxg (2.30) 

Using Newton’s second law gives 

L = r x F = N. (2.31) 

Notice that angular momentum may be a constant of the motion relative to one origin but 
not to another for certain force laws. Consider uniform circular motion: 



Here, L is a constant of the motion measured from 02 (since r and P are co-linear) but 
not from Ol. This will be the case in general if (“central force”) 

F = e r f(r), (r = |f|) 
since then 

N = r x e r f(r) = 0. 


The definition of work done on a single particle in moving it from position 1 to 2 is given as 
W 12 = J' F • dr , (2.32) 


where F is the total force acting on the particle. It is understood the integration takes 
place along some specific path, as below. 
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(The motion is shown to lie in a plane for convenience.) Since F is the total force, by 
Newton’s second law 


F • dr = m 

dv 

dt 

dr , 

• — dt 
dt 

dv 

= m — • 
dt 

v dt 

m 

2 

d 

dt 

(v 2 ) dt. 

Eq.(2.32) 

=> 

W = 
vv 12 

- mv 2 | 2 

2 

= -m( 
2 v 

- 

Vr 2 ) 

1 = T 2 - T lf 


(2.33) 


where (“kinetic energy”) 


T = 


- 2 

— mv . 
2 


(2.34) 


Eq.(2.33) generalizes the one-dimensional case, (2.13), above. The left hand side of (2.33) 
is again the work done. 


Let us assume the special case, 

F (x, x, t) -»> F (x) . 

Let us also remember the connection between Newton’s force law and energy conservation 
we saw previously in one dimension: the energy conservation statement was the first (spatial) 
integral of Newton’s law. We have 


mx 


= F (k) , 


mf k • k dt dt ' • F (k) , 


dt ' 


\ m (x 2 - x 2 ) = jf dx '• F (x ') . 


Let us define, (“potential energy”) 


U (x) - -f dx'- F (x') , 

A 

: " reference point" . 

At this point we will make an assumption about U (x) :. We will assume, 


(2.35) 


(2.36) 


U(x) : A single-valued function of its 
argument, x. 
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Why do we do this? Consider: 



That is, assume a force law of the form, 


F = Fe, . 


Then 


J F • dx = 2 jtRNF. 

N turns 


N = number of turns around the circle 

Clearly, this force law does not result in a single-valued function of the final position. Thus, 
a single-valued force law does not necessarily imply a single-valued potential. With this 
assumption about U(x) , we have 


\ m (x 2 - Xo) = -U(x) + U(x 0 ) . 


Rearranging gives 


1 12 _ 1-2 

- mx o + U(x 0 ) = - mx + U(x) 

=> (t + u) L = (t + u) | t . 


(2.37) 


This leads to the third conservation law. 


III. The total energy a particle, 

E = T + U(x) , 

subject to a single-valued potential, U(x), is conserved. 

dE 

Paraphrase: U(x) single-valued => - = 0 . 

dt 
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Proof: By the above direct construction we have 


E(t 0 ) = Eft) . 


Since t is arbitrary, 



dt 


Note that the potential U(x) always needs a reference point for it’s definition. If we change 
it, we get 


X 

U new (x) = -J dx' • fix' ) 

r2 

x r 1 


= -Jdx' • F(x' ) 

ri 

U(xj 


J dx' • F(x' ) . 

r2 

constant in x 


(2.38) 


Thus, the old and new potentials just differ by a constant. This means the total energy of 
the particle, E, is arbitrary up to an overall constant. Clearly, given U(x), this arbitrary 
constant does not affect F(x) itself since by the Leibnitz rule, 
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- V U(x) = vf dx'- F(X ' ) = F(x) . (2.39) 

Since E is conserved and using Eq.(2.33), an alternate way of writing the work done in 
going from from x l3 to x 2 is 

W 12 = U(x x ) - U(x 2 ) . (2.40) 

Notice that the work done depends only on the initial and final positions and is therefore 
path independent. This is equivalent to the above assumption of single-valuedness for 
U(x). In addition, a force field which can be represented in the form of (2.39) is said to 
be conservative. This is again equivalent to the previous ideas of path independence and 
single-valuedness. This is summarized as 





Each implies (or assumes) the other; any of these concepts can be invoked in writing down 
conservation law III above. 3 

Potential Energy and Particle Motion 

We are familiar with the form of U(x) near the Earth’s surface: 


3 

( 



1 

, f = -mg 


/ / / / / / / 


We surmise that U(x) = mgx 3 . Is this correct? 
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F = -VU = -mg Vx 3 = -mge 3 . 

This is indeed correct. Notice that U = mgx 3 + const, would do just as well. 

With the above connection between F and U in a conservative force field, we can recast 
the solution for the one dimensional case when F(x, x, t) —* F(x) Eq. (2.16), in a different 
form. We have 


dU 


/ dx F(x) = ~f dx — = U(x 0 ) - U(x' ) . 


Since the total energy, E, is a constant of the motion, we can evaluate it at a single time, 
say t 0 , and it will continue to take on the same value at all later times. Therefore 

1 -2 

E = - mx° + U(x 0 ) , 


2 x r 2 2 

— I dx F(x) + x 2 0 = — (E - U(x ' ) ) . 
m „ m 


(2.41) 


and (2.16) becomes 


X(t) 


*/ 


dx 


(E - U(x ' ) ) 


= (t - t 0 ) 


(2.42) 


m 


Let’s find out how to use this equation by relating it to a picture of a general one-dimensional 
potential, U(x). 


U (x) 

E 4 "unbounded" 


"turning 


> t 3 X 

1 _ "turning 

"2 \ / \ s ^_ J 

_ 


"bound" 

"bound" 




x a x n x. x„ x 

a 0 b c 
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The potential U(x) is shown as the curvy line. On the same graph I have drawn horizontal 
lines which represent total energies Ex through E 4 . First of all, let’s examine the motion 
of a particle with total energy El. Clearly, the particle which has energy E = E 4 is bound. 
The two points labeled x a and x b are called “turning points”, mentioned previously in the 
harmonic oscillator discussion. They represent, for bound motion, the extreme values in 
position reached by the particle in it’s motion. x 0 is called a “stable equilibrium” point 
and will be discussed below. Note that a particle with total energy E 2 has both bound 
motion (in two different locations) and “turning” motion with a single turning point. The 
turning motion begins at position x = x c ; it represents where the velocity of the particle 
changes sign. In general, turning points are positions where E = U, which implying the 
kinetic energy, T, is zero. The motion of a particle with energy E 3 is separated into two 
regions with turning points as shown. It is also “turning” in the sense that the motion is 
not restricted to any finite region in x. The motion of a particle with energy E 4 . takes 
on all values -oo < x > oo along the x-axis. I will use the term “unbounded” to describe 
motion without any turning points. 

How would one apply (2.42) to the motion of a particle with energy E l5 say? First of all 
recall that in deriving (2.42) we assumed that the relation between t and x was one-to-one. 
This is clearly violated for periodic motion, as we saw previously, but remember we can 
always piece together various separate motions where the relation is one-to-one to describe 
the full motion. The rule here is that these periodic motions are separated by turning points. 
Let us say the initial conditions were x(0) = x 0 , x(0) > 0 • Thus, clearly one would 
want to choose the positive sign in (2.42) to describe this part of the motion since t > 0 
and x (t) > x 0 : 


x (t) 

E = E lf x(t) > x 0 : t = / . 

*0 J— (E - U(x' ) ) 

Vm 


(2.43) 


On the other hand, after the particle reached x b , it’s motion would reverse. For the same 
starting time as in (2.43) the time, t, associated with a point x(t) on this part of it’s 
motion would be given by: 


E 


E x , x(t) < x b : t 



(2.44) 


t b is the time necessary to reach the first turning point, x b : 

dx' 


^ = / 


i 


(E - U(x' ) ) 


(2.45) 
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Notice the minus sign in (2.44). It is required by the statement that t > t b for this part 
of the motion given that x (t) < x b . Clearly, by patching such equations together in a 
similar way, one can find the relation between x and t at any time. Similar observations 
can be made about the motions of particles with energies E 2 , E 3 a nd e 4 corresponding 
to their various bound or free motions. 


If one is interested in the period of bound motion in a one-dimensional potential, one has 
to add together the time intervals for the x > 0 motion and the x < 0 motion. These 
are given by (again referring to the figure) 


At, 


■ J 


dx' 


and 


f. 


(2.46) 


(E - U(x' ) ) 


At 


-f 


dx' 


i 


b , E - U(x' )) 


(2.47) 


respectively. Obviously, At ab = At ba ,, which is true even if a and b are not turning points, 
and one has for the total period T = At ab + At ba , 

dx' 


T = 2 


J 




(E - U(x' )) 


(2.48) 


where x a and x b represent the turning points. The turning points are to be found by 
setting the kinetic energy equal to zero, which results in 


E=U(x a , b ) , 


(2.49) 


where x ajb are the two possible solutions. 


One can always use (2.48) to write the period as (take x 0 = 0) 


T = 2 


J 


dx' 


J 


dx' 


l-(E - U(x' )) o /-(E - U(x' )) 

Vm llm 


For a potential even about the origin, U (x) = U ( -x) , 


X a = “X b , (X b > 0) 


and changing x' —■» -x' in the first term as the integration variable now shows 
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T = 



dx' 


I 


“(E - U(x' )) 
m 


As an example, consider Hooke’s law, where one may take U(x) = — kx 2 , k > 

I2E 2 

the period is given by (2.50) as ( x' = —y is the change of variables) 

V k 



This is just the result we found earlier in this Chapter for a Hooke’s spring. 
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Equilibrium and Stabilty in One Dimension 

In the above discussion, I invoked the idea of a stable equilibrium position. Let us examine 
the meaning of this concept. We will be asking the question: Is the position of a particle 
stable under an infinitesimal displacement? The following is an illustration of the concept. 



Assuming the potential U (x) may be expanded in a Taylor series about the point x 0 shown 
above, we have 


U(x) = U 0 + (x - x 0 ) ^ | X() 


(x - x Q ) 2 /d 2 U 


v dx z ; x o 


+ • • • 


( 2 . 52 ) 


The constant U 0 will not affect any aspect of the particle’s motion and can be eliminated 
by redefining U(x). For the above illustration, 


dU d 2 U 

xn H f , 2 x n 

dx dx 


> 0 . 


Motion with energy E slightly larger than U 0 will then result in bound motion, and we 
take this as meaning stable equilibrium. If 


d 2 U 

dx 2 


< 0 , 


the potential would have looked like: 


U(x) 


U„ 


x 
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Clearly, there is no stable equilibrium in this case. A position x 0 where 


du i 

dx l »' ' 


is called an extremum of U (x). 


Under the assumption on the expandability of U ( x ) in a Taylor series, we find that 


dU . 

- x = 0 r 

dx 1x0 


(C.l) 


is necessary (but not sufficient) for stable equilibrium. If we have in addition the condition, 

(C.2) 


d 2 U 


dx 5 


> 0 , 


then (C.l) and (C.2) together constitute sufficient conditions for stable equilibrium. (C.2), 
however, is clearly not necessary itself since we might encounter 


d 2 U 


dx 


2 Ixo 


= 0 . 


Under these conditions, one must then look at the first nonzero derivative of U ( x ) at x 0 
to decide the matter. 


These considerations allow us to solve for the small oscillations near an equilibrium position 
for a general one dimensional potential. The force at x 0 is of course 


F ■ - S u i- 


(2.53) 


and from (2.52) near equlibrium this gives 
F « - (x - x 0 ) 


d 2 U 


dx 


2 lx 0 * 


(2.54) 


Choosing the origin at x = x 0 now tells us that to lowest order 


m x = F « -x 


d 2 U 


dx 


2 lx 0 * 


(2.55) 


Thus for small oscillations we are back to the Hooke’s law equation m x = -k x where 


k - 


d 2 U 


dx 


2 lx = x 0 


as long as condition (C.2) holds. This means the formula (2.51) holds 


and that the angular frequency of small oscillations is determined by 


2 k 1 d U 
O ) 0 = — = — 


m m dx 2 1x0 ' 
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Although the usual meaning of the x variable is spatial position, it can also represent other 
things. In an end of Chapter problem it is seen that these ideas are aso useful when x 
represents an angle. 

If a Taylor series expansion can not be made for a point because of discontinuities, one 
must inquire more generally about the meaning of a stable equilibrium point. The more 
general meaning is associated with the turning point idea. Let us imagine increasing the 
energy, E, of the particle infinitesimally higher than Uq in the three following examples. 
(Increasing E is the same as giving an initial infinitesimal velocity to the particle.) 



2 turning points 
(a) 


1 turning point 

(b) 


no turning points 

(c) 


Clearly, for the motion to be stable, there must be two turning points, labeled (1) and (2) 
as in (a) above, located infinitesimally close to x 0 . If there is only a single such turning 
point, as in (b), the motion will not be bounded in one direction and we have an unstable 
equilibrium. If there are no turning points, as in (c), this also represents unstable equilibrium. 
Note that if there are 2 turning points, but one or the other is located a finite distance from 
Xq, resulting in bounded but finite motion, we also consider this as unstable equilibrium. 


Given that I = 
dx Xo 


0 , the more general meaning of stable equilibrium is: 


Under an infinitesimal increase in the total 
energy E = U(x 0 )+ 6 e, there are two turning 
points located infinitesmaly close to x 0 


dU 

A point, x 0 , which has - | = 0 but which does not satisfy the above criterion is 

dx 0 

considered an unstable equilibrium point. 


There is also the possibility of neutral equilibrium, by which I mean a point x 0 where 
all the derivatives of U (x) vanish. Then of course any particle motion is simply uniform 
in the vicinity of x 0 . 
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Equilibrium and Stability in D Dimensions 

What about equilibrium in multidimensional problems? If we again assume that U(x) can 
be expanded about the position x 0 in a Taylor series, we have (understood sums up to 
dimension D on repeated indices) 


U(x) = U(x 0 ) + (x - x 0 ) i V i UL + -(x - x 0 ) ± (x - xWVW U 




(2.57) 


Again, the constant, U(x 0 ), can be ignored in the discussion of particle motion. The analog 
of the extremum condition (C.l) is now 



(C.3) 


This is again a necessary, but not sufficient condition for stable equilibrium. Notice that 
(C.3) represents D conditions in a D -dimensional problem. 
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What is the analog of (C.2) in this case? We will just generalize the argument of the last 
section to more than one dimension. Assuming (C.3) holds and using (2.57), the force on 
the particle at x = x 0 is then given as (I’m using index notation) 


F, = -V,U 



U 


= - - V.V. U I- 

2 3 0 


d(x - x 0 ) ± d(x - x 0 ). 

- (x - + (x - x n ). 


dx, 


dx t 


= 5 . 


it / 


- " ^ ( V < V , » k ) < x " x .>, " “ o l«.) < x - x .)i 


= -A 0 (x - x 0 ). , 


(2.58) 


where 


A fj - U(x) 


(2.59) 


These set of constants, A fj = A j( , , are the analog of the single spring constant, k, in the 
one-dimensional case. Choosing our origin at x = x 0 , the equations of motion are now 


mx f = F, = -A /jXj , 


=> mx. + A t . x. = 0. (2.60) 

Tliis represents a set of coupled linear differential equations in the quantities Xj. Let us 
suppose a solution of the form 

x, = e imt . (2.61) 


We are now using complex numbers; the a^ are complex (and the A^ are real) but it is 
understood that x^ is represented by the real part of the right-hand side of Eq.(2.61). Then, 
the equations (2.60) become simply 

(A^.-mco^.) a. = 0, (2.62) 
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where the -to factor in the second term comes from taking the second time derivative 
in (2.60). Eqs.(2.62) represent D linear, homogeneous equations in the D unknowns, aj 
. As was quoted in the first Chapter, the condition for such systems of equations to have 
nontrivial solutions is that the determinant of the coefficients of the unknowns must vanish. 
Thus, defining 


C i3 " ^3 " 


mto 2 h 


13 


(2.63) 


we must have 


det (C) = 0. (2.64) 

Eq.(2.64) determines (the squares of) the characteristic angular frequencies, 0 ) , which 
enter in (2.61) above. In D dimensions, the determinant equation, (2.64), represents D 
equations, giving D solutions, to 2 , to 2 , . . . , to 2 . 

In a problem at the end of the Chapter we will see that the to are real quantities. Let 

us examine what would happen if one of these frequencies had a nonzero imaginary part. 

2 

Notice that since (2.64) is actually a condition on to , not to, the supposed complex roots 
occur in pairs, 

to = ± ito c , 

where to c is real. Thus, in this case the general solution would read 

/ -iv — COn t CO v 

Xj = a { (A e c + B e c ) , 

where A and B are constants. Requiring 

x(0) = 0 => B = -A, 

' x(0) = v 0 => a ± A = - , 

2to c 

v n 

=> Real(x i ) = —- sinh(to c t), 
to c 

=> lim Real x ± -> ± 00 . 

t -»oo 

Therefore, the motion is unbounded for to c ^ 0 . 
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A sufficient condition for stable equilibrium in higher dimensions is thus, 


Given (C.3) holds at x 0 and the squares 
of the characteristic angular frequencies, 
determined by det(C)=0 are positive, 
the potential possesses stable equilibrium 
at x 0 . 


(C.4) 


The characteristic angular frequencies, 0 ) k (k=l,...,D) are analogous to the single 
natural frequency, wO, in the Hooke’s law case. As we found, 

2 k 1 d 2 U 

w 0 - - , 2 lx 0 * 

m m dx 

'-T-I 2 

Thus the condition in (C.4), “the 0 ) k are positive”, is clearly a generalization of condition 
(C.2), which required 


d 2 U 

dx 2 

for stability. 


> 0 , 


Let us look at one example of the use of (C.4). 


Example 5 

In two dimensions, a particle near the origin experiences the potential, 

1 2 1 2 

U(x, y) = -k x x + — k y y + k xy xy, 

where k x , k y , k xy >0. Find the characteristic frequencies. Does the particle posess 
stable motions? 


Solution: 

Clearly the origin is an equlibrium position. We compute the nonzero Aj_ j: 


An — k x , A 2 2 — k y , A i2 — A 2 i — k xy . 


We have 


det 


k - moo 2 


xy 


k y - Itlft) 


= 0 . 


a) 2 = — (k x + k y ) ± — y(k x + k y ) 2 + 4(k; y - k x k y ) 


= V< k x - k y) 2 + 


4k 2 


xy 
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Notice we will have real solutions for co if and only if 

k 2 xy < k x k y , (E.52) 

which is the condition for stability. 

There are many other considerations related to stability in higher dimensions we do not 
have the time to cover here. One is often interested in stability under certain assumed 
types of small perturbations, asssociated with the various squared eigenfrequencies, oo k . 
One concept seen later is the idea of an “effective potential” for a moving particle in an 
external field. It is possible to examine such a potential for stability of a given orbit rather 
than particle position. 

CHAPTER 2 PROBLEMS 

1. Our sprinter friends Michael and Frankie run another 200m race. In this race 
both runners accelerate uniformly and then continue to run at different top 
speeds. If Michael’s top speed is 10 m/s achieved after 15m and if Frank’s 
acceleration is only 2/3 of Michael’s, how far must Frank accelerate to tie 
Michael in the race? What is Frank’s final top speed? 

[The methods, not the answer, are important, so show me your steps!] 

2. A steady wind is blowing on a pendulum hanging from a patio roof. The 
pendulum bob has mass "m" and the wire is approximately massless. The 
downward force of gravity is — mg. 


/// //////// 



The wind exerts a horizontal force, f , on the pendulum. Find the equilibrium angle, 0. 
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3. A baseball of mass m is thrown vertically during a game of catch. It experiences 
a downward force F = - mg - kv, where - mg is the force due to gravity 
and - kv, represents air friction (v is the baseball’s instantaneous velocity, k 

> 0). Given that the baseball starts at the origin of coordinates at t=0 with 
initial upward velocity Vo, solve for the position of the particle at all subsequent 
times, x(t). 

4. A particle is dropped in a medium for which the resistive force is given by 

D x 2 ( D> 0 ) where x is the particle’s velocity. The force of gravity is -mg, as 
usual. Investigate the time taken to fall a distance h from rest. Show that to first 
order in D, this is given by 



State the approximation that is being made here. 

5. Define Rh = horizontal distance between 2 pts. on the trajectory, a distance h 
above the horizontal. The particle is projected at an initial velocity of v 0 at an 
angle a above the horizontal. 
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a) Find the distance Rh in terms of h, v 0 and cc 

b) Find the angle of elevation, a, which maximizes the distance for a given 
velocity v Q , assuming zero air resistance. [Hint: take the derivative of R^ 
with respect to a.] 

6. A particle of mass m is subject to a force (one dimensional motion, k>0) 
F(x,t) = -kx + F 0 sin(o)t), 

where <x > 2 ^ ~ It also has initial conditions, 


x(0 ) = x 0 , 
x(0) = v. 

Find x ( t ) for t > 0 • 

7. Nolan Ryan (in his prime) throws a ball up an incline that makes an angle p 
with respect to the horizontal, as shown. The angle that the ball is thrown is a.. 
[Neglect air resistance.] 
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a) Given the initial magnitude of velocity, Vo, find the time it takes to hit an 
uphill point on the slope. 

b) Nolan wants to throw the ball as far as possible up the slope. Given the 
same initial velocity as in (a), at what angle a should it be thrown? [Ans.: 



[Hints: One thing that helps is the idea that maximizing the distance up the hill is the 
same as maximizing the distance along the x-axis above.] 

8. A magnetic monopole is a hypothetical particle with a magnetic rather than 
an electric charge. Assume a very massive magnetic monopole is sitting at the 
origin, as shown. A very light electric charge, q, is located at f relative to the 
origin and has a nonzero velocity. 



The magnetic field of the monopole is given by 
e r 

B = g r 2 , 

where e r is the unit vector point in the f direction and g is the magnetic charge. The 
force on the charge q is given by (the Lorentz force law) 
q J. 

F = — r X B . 
c 

Show that the equations of motion in cylindrical coordinates are 
(C*®* and r 2 ^p 2 +z 2 ) : 

C 

p: p-p4> 2 =—p4>z, 

r 

C • 

4): 2p 4) +p 4) = r3 (z p-p z), 

C_ 2 - 

z: z — r 3 P “4> . 
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[Hint: Recall the formulas for velocity and acceleration in cylindrical coordinates from 
the last Chapter.] 

9. Show that the motion of the monopole in problem 8 takes place on a cone, 
the end of which starts at the stationary magnetic charge g. Hint: A cone is 
characterized by 

p = z tan (a), 
where ct is the opening angle: 



10.A bead is released (at rest) from the top of a ramp of height "h" as shown. It 
moves without friction until it reaches the edge of the cliff, where it’s velocity is 
purely horizontal. The cliff is also of height h. 



Find the distance, x, that the bead moves horizontally before striking the ground in 
terms of h. 

11.A one dimensional potential, given by U(x) = — kx 2 , 
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U( X) 



is split in the middle at x=0 into two similar parts by inserting a flat segment of length 
L, as shown: 


new U(x) 



Find the period of motion of a particle with mass m and energy E in this new potential. 

12.A particle in one dimension experiences a potential, U(x) = K | x | n , where 
| x | is the absolute value of x, K is a positive constant and n > 1 • 

a) Find the oscillation period of the particle in the potential for general n. [It is 
enough to set up the integral which will give this answer.] 

b) If the energy of the particle is increased by 2, by what factor does the period 
change? [Flint: Get the integral in dimensionless form.] 


13. A one-dimensional potential is given by 


U (x) = 


1 , 2 

— kx , x > 0 
2 

1 2 

-kx , x < 0 . 

2 
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Is a point particle placed at the origin stable or unstable? Explain carefully. 

14. A one dimensional potential is given by, 

1 , 1 

U (x) = - -kx 2 + -dx 4 . 

2 4 

a) Find the points of stable and unstable equilibrium. Tell me which are stable, 
which are unstable. 

b) Find the turning points for a particle of energy E>0. 

15. A point particle is restricted to motion along the x-axis. The potential energy of 
the particle is given by 

U ( x ) = U 0 ( 1-cos ( — ) ) • 

x n 
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This looks like: 



a) Find the points of stable, unstable equilibrium. 

b) If the particle is started at x=0, what minimum velocity is necessary to reach 
the region x > Jtx 0 ? 

c) Find the frequency of small oscillations about the point x=0. 


16. The force as a function of position on a particle of mass 1 kg is given by 


F 


-2, x > 1 
-2x,-1 < x < 1 
2, x < -1 


a) Find the potential, U (x). 

b) Find the period of motion the particle in this potential if the particle starts at 
Xq = 2 with zero initial velocity. 

c) Find the period of small oscillations of this particle. 


17. A particle of mass m follows the trajectory y = Kx 2 (K is a constant) in the xy 
plane under the influence of gravity as shown. The amplitude of oscillation of 
the particle is small. Find the period of oscillation of the particle. 
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18. A weak unstretched watch spring is attached to a pendulum of length E with 
a mass m at the end. Gravity is present and the equilibrium position is at 
4>=0. The potential energy of the mass is given by U((j)) = mgC( l-cos({)) 

+ -k d) 2 ’ where k s is the (rotational) spring constant. Assume k s « mgE. 

a) Show that the squared angular frequency of small oscillations about the 
equilibrium angle, (|)=0, is 





position 


b) The spring is twisted approximately 2jt radians and let go. Show that the new 
stable equlibrium angle, 4>o, is given by 


- 2it (1 - . 

mql 

c) Show that the frequency of small oscillations about the new equilibrium 
position is approximately 

2 g k 2 jt 2 k 

rn = _ + __ _ l 

new i mf 2 gm 2 f 

19. The (squared) characteristic frequencies of a system, oo 2 , are determined by the 
system of equations, 


A t . a. = ma) 2 a,, , 

(understood sum on j) where the a ( are complex quantities but the A^j are real and 
symmetric in E and j One may prove that the square of the characteristic frequencies, 
0 ) 2 , are real quantities by completing the following steps: 
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1) Multiply by Sl f on both sides of the above and sum on £. 

2) Take the complex conjugate of the result of (1). 

3) Now complete the argument that (w 2 ) = 0 ) 2 . 

20. A particle is subject to the two dimensional potential, 
U(x,y) = jk(x 2 + y 2 ) + Hyx 2 , 
where k and H are constants (k>0) . 


a) Find the equilibrium positions of the system (there is more than one). 

b) Examine the equilibrium positions in (a) for stability by calculating the 
characteristic angular frequencies of the system. Which positions are stable? 


21. A particle is subject to the three dimensional potential, 


U( x,y, z) = -k (x2 + + *») + Hxy + Gyz , 


where k, H and G are constants (k>0). Find the characteristic angular frequencies of 
the system. Show that k > n/g 2 + H 2 is necessary for stability. 
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3 LINEAR OSCILLATIONS 


GENERAL RESTORING FORCES IN ONE AND TWO DIMENSIONS 

We already looked at a spring example in the last Chapter. Let us imagine setting up 
such a system in a gravitational field, as below. (Assume the motion is approximately one¬ 
dimensional.) 


+x 


f / / 


Equilibrium 

position -► 

(origin) 


/ / / / / 



-mg 


Displace the mass from equilibrium or give it an initial velocity. Newton’s law is 


ma = F tot , 


or 


mx = F sp (x) - mg. 


(3.1) 


in this case. Do a Taylor series expansion on F sp (x) : 
F sp (x) = F sp (0) + x 


+ X 

(dF \ 

sp 

x 2 

+ — 

KM 


l J 0 

2 

l a x z j 


(3.2) 


The above coordinate system has been defined so that the total force vanishes at the 
equilibrium position, 


F sp (0) - mg = 0. 


Thus 


mx 


/ dF 

\ x 2 

( a 2 F 'i 

sp 

i + — 

sp 

\ ax 

Jo 2 

l a x z ) 


(3.3) 


(3.4) 
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Defining (“spring constant” again) 


k - - 


dx J' 


(3.5) 


and neglecting higher-order terms in x in Eq.(3.4), we have the approximate, small oscillation, 
equation of motion (assuming k > 0) 


x + — x = 0 , 
m 


(3.6) 


which we have seen and solved previously. The most general solution involves two undetermined 


constants and can be assumed in the form ( co 0 = — ) 

V m 

x(t) = A cosu) 0 t + Bsino) 0 t, 


(3.7) 


or 


x(t) = Csin(oo 0 t - 6) 


(3.8) 


It’s easy to find the relationship between the coefficients in (3.7) and (3.8). We have 
Csin(o) 0 t - 6) = C cosh sin (o) 0 t) - C sin 6 cos ( oo 0 t) , 


fA = -C sinh 
IB = C cosh 


2 2 2 

A + B = C . 


(3.9) 


(Note that we may always choose C>0 . ) Of course, D cos (o) 0 t - h' ) is also a possible 
general solution. 


As usual, the first integral of the force equation with respect to x gives work done and yields 
the energy conservation equation: 


J dt x [x + WqX = oj , 

=> — x 2 + — (On x 2 = const. 
2 2 

Multiply by m and identify the terms: 


1 -2 1,2 
— mx + — kx = const. , 
2 2 


or T + U = E . 


(3.10) 


(3.11) 
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Usual connections (x is the period, v 0 is frequency): 


oo 0 T = 2 k 



v 


o 



1 

X 


Let’s consider forces in two dimensions. 


(3.12) 



F 

y 


y 


equilibrium 
/ position 


x 


Let’s assume 


f x ( x ' y) = F x(x) 

F y (x, y) = F y (y) ,1 


(3.13) 


so the motions in the x and y directions are independent. Again, for small amplitudes 
around equilibrium 


F x (x) - FJO) + x 

"IT 

F v (y) z F v (o) + y 


l dx ) ' 

l ay) l ^°' 


Define the independent spring constants 


k = - 


dK 

dx 


V- — _ 

ic=0 ' y ~ 


d A. 

, 8y 


ly = 0 f 


and assume they are positive. Then we have the equations of motion 


x + 0 )^ x = 0 , 
y + 0 ) y y = 0 . 


(3.14) 
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where 0 ) 2 




The general solutions are 


x(t) = A cos (w x t - a),] 
y(t) = B cos (a) y t - |3) .J 


(3.15) 


Special case: 0 ) x = a) y = ( 0 , A, B>0. Then we get the following behaviors as a function of 
6 = a - (3, the phase difference. 




8=0 


6 


jt 

2 


(a) 


(b) 



(c) 


(d) 


Notice the (b) and (d) trajectories are the same except the sense of rotation is opposite. In 
the equal frequency case, one may show that 

2 2 A 2 A 2 2 

x + y — - 2xy — cos 5 = A sin 5 . (3.16) 

B B 

Now by introducing new rotated coordinates x' and Y' r 

x = x 1 cos y - y 1 sin y, 
y = x 1 sin y + y 1 cos y. 
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and substituting in (3.16), one may now show with the choice 
2AB cos 6 

tan 2y = —--— , 

A 2 - B 2 

that (3.16) can now be written as 


(3.18) 


Cx' 2 + Dy' 2 = 1 , 

where 

B 2 cos 2 y - A 2 sin 2 y 
A 2 B 2 sin 2 6 cos 2y ' 

A 2 cos 2 y - B 2 sin 2 y 
A 2 B 2 sin 2 5 cos 2y 

(3.19) is the equation of an ellipse, if sin 6 ^ 0. In the case 5 



Let’s say now that 0 ) y = 0 ) x + Aw. Then 


(3.19) 


(3.20) 

(3.21) 

= ± — , we get for example 
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X(t) = A COS (0) x t - a) , 

y(t) = B cos (oo x t + (Aoot - (3) ) . 


( 3 . 22 ) 


Effectively, 6(t) = a - (3 + Acot , a time-dependent phase factor now. If Aco « 0) x 
one can think of an evolution in time from diagrams (a) to (d) above as 6 changes from 
0 to 2jt. This evolution occurs more quickly for larger Aco . 

Of course, many more types of motions exist when Wx^Wy. Consider two cases: 

1. —— = rational number. We are guaranteed in this case to generate a closed 

« y 

trajectory. Take the case, 

co^ _ 199 _ 199 

to y ” 200 ^ T x 200 ' 

which means it will take 200 repetitions of the T y period to equal 199 exact repetitions of T x , 
after which the motion repeats. These are called “Lissajous curves.” Some of these are 
quite beautiful. 

2. —— = irrational number. In this case, the motion will never repeat and every 
Wy 

point in the range —|A| <x<|a|, — |b| <y < | B | will be reached by the 
motion in principle. 

Damped oscillations 

Next, imagine making the one-dimensional situation more realistic by adding a damping 
force. 
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There are now two horizontal forces acting on the mass, the spring, F sp , and the frictional 
damping force, F d . We notice that the damping force must reverse sign when the velocity 
reverses. We will set 


F tot K *) = F sp ( X ) + F d ( X ) f 


(3.23) 


F sp = -kx , k > 0, 
F d = -bx , b > 0. 


(3.24) 


Notice that a force proportional to x' would not reverse sign as we desire. Our equation 
of motion is now of the form 


mx = -kx - bx , 


x + 2|3x + oo 2 x = 0 , 


(3.25) 



(3.26) 


Because of the damping of the motion, we no longer expect the energy (mass kinetic plus 
spring potential) to be conserved. However, we can still calculate, as usual, the work done 
on the system in the usual two ways (W 12 denotes the work done in changing the system 
from state 1 to state 2): 



W,. 


1 /'2 
-m(x 2 



(3.27) 


or 


W,. 


2 T 

= -J (kx + bx)dx =-k (x 2 2 - x 2 ) - bj dt x 2 . 


ti 


(3.28) 


If we use t=0 as the starting time of the motion, (3.27) must equal (3.28), and we find 

t 2 tq 

E (t 2 ) + bj dt x 2 = E (tj + bj dt x 2 , (3.29) 

0 0 

where E (t) = — mx 2 + — kx 2 gives the instantaneous system energy. (This procedure is 

equivalent to integrating the equation of motion, (3.25), as we did before to get (3.10).) 
Eq.(3.29) informs us that 
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E d (t ') = b f dt x 2 , 


(3.30) 


is the energy lost during the motion because of frictional damping, and the total energy, 
including heat, is conserved. Thus, 


E(t) + E d (t) = const., 

and the instantaneous change in the system energy (kinetic plus potential) is 

-b*U 

dt dt 

The right hand side is negative forb>0, as it should be. 

How does one solve the differential equation (3.25)? We will substitute 


rt 

x = e , 


x = re rt , x = r 2 e rt , 


(3.31) 


(3.32) 


giving the “characteristic equation”, 


r 2 + 2|3r + oo^ = 0. 


The above is a quadratic equation, the roots of which are given by 
r = -(3 ± Jfi 2 - (a 2 0 . 


(3.33) 


(3.34) 


The factor in the square root can be positive, zero, or negative, leading to the following 
three cases. 


(1) w 0 > (3 (“underdamped”) The roots of (3.34) are complex. Define 
>i = V w o - P 2 > 0 f 


0 ), 


r = -|3 ± io) 1 . 


(3.35) 


(3.36) 


The two solutions are: 


_-flt in>it _-fit _ -io)i t 

e e , e s 


or 
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e pt cc^ccqt) , e pt sin^t) , 

When x ( 0 ) > 0 and x ( 0 ) = 0 , The motion looks like: 



The system has no precise frequency since the motion is not strictly periodic. However, there 
is a type of period involved since aspects of the motion reoccur, although the amplitude 
decreases continuously. In particular the points where x = 0 seem to reoccur periodically. 
Where does this happen? Let’s build in some B.C.’s to examine the motion in detail. Use: 

x(0) = x 0 , x(0) = v 0 . 


Assume 

x(t) = e _pt [A 1 e i “ lt + A 2 e- iMlt ]. 
x (0) = x 0 => x 0 = A, + A 2 , 

x (0) = v 0 => v 0 = e- pt {-(I [a^ + A.e^] 

+io) 1 [A 1 e iMt - A 2 e“ imt ] } | t = 0 

v 0 = -P[a x + A 2 ] + icoJX - A 2 } 


Solving (3.37) and (3.38) for A1 and A2 gives 
v 0 + |3x 0 + ib) 1 x 0 


Ai = 


A 2 = 


2ica 1 

-v 0 - ftx 0 + ico 1 x 0 
2ia), 
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Introducing e 1 11 = cos oc^t ± i sin oOjt, the solution with the above B.C.’s reads 


(v 0 + |3x 0 ) . , 

x(t) = e 1 [--- sin(o) 1 t) + x 0 cos(o)jt) ] . 


O), 


(3.41) 


Now investigate when x = 0 : 


_ -pt"P (v 0 + Px 0 ) 


0 = e“ p [- 


ca, 


sin oc^t - |3x 0 cos oc^t 


+(v 0 + |3x 0 ) cos cc^t - XgiDj sin oc^t] , 


pv 0 + |3 2 x 0 + x 0 w 2 „ _^ 

-sin (Dj^t = v Q cos • 

(D, 


(3.42) 
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When v 0 = 0 , then (0 1 t — nrt, n—0,1,2,..., just like the undamped ( P — 0 ) case. 
The picture is: 



Motion is “quasi-periodic.” When v 0 * 0 the motion is still quasi-periodic (Ajnqt ) = Jt) , 
but we now get oqt = mt +const., n = 0,1,2 , . . . 

(2) (3 = o) 0 (“critically damped”) Now one has the repeated root, 

r - (3.43) 


and so the solutions are 


e 


-pt 


r 


te 


-Pt 


Solutions are non-oscillatory, although x (t ) is permitted one change of sign under certain 
circumstances. It’s fairly easy to see what the condition is. Assume x 0 > 0 (|>, t already are 
also). The only way to make x(t) ^ 0 is if we choose v 0 negative in the general solution, 

x (t) = [ x o + ( v 0 + pXoh] e_pt • (3- 44 ) 

If we choose v 0 ^ -|3x 0 , we see from the above that x(t) is always positive or zero. 
However, for v 0 < -px 0 we must have x(t) < 0 for some values of t since 

lim x(t) -* (v 0 + px 0 )te“ pt . (3.45) 

t—*oo 

Thus, by continuity the minimum magnitude of downward v 0 which results in a change 
in sign of x (t) is | v 0 | min = P x o r toward the origin. 
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(3) o) 0 < p. (“overdamped”)- Define 
0 ) 2 = -^p 2 - oo 2 > 0, 

=> r = -p ± oj 2 , 

so there are two real roots. Solutions are: 


(3.46) 


Also non-oscillatory; U) 2 is not the frequency of any motion (as 0 )! was in case (1)). x(t) 
is again permitted one change of sign under certain conditions. I will refrain from writing 
down the general solution in terms of x 0 and v 0 because this will be a problem. The 
general form is 

x(t) = e“ pt [A ie “ 2t + A 2 e~“ 2t ] . (3.47) 

Circuit/Oscillator Analogy 

We will now point out a useful analogy between mechanical oscillator and electrical circuit 
combinations. This will be deduced from the differential equations that describe these 
different situations. 


The analog of 



(no damping) 
mx + kx = 0 


is the LC circuit (define positive direction of current to be counterclockwise) 
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C 


KirchofFs law is: 

V L + V c = 0 , 

1 r 

LI + - I dt = 0 
r J 


L q + — q = 0. 
C 


Another example: 


(I = q) 



Compare with: 



V R : RI . 


(3.48) 
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Lq + Rq + — q = 0. 
C 


(3.49) 


Here are the analogous quantities: 


Mechanical 

Electrical 

force ,F 

voltage,V 

x (mass or spring 
displacements) 

q (inductance or 
capacitive charge) 

X = V 

<3 = 1 

m 

L 


1 

k 

C 

b 

R 


By “displacements” or “charges” above I mean values measured from equilibrium. The above 
Table now gives: 


( 0 , 




1 

Vlc 

R 

<-> — . 

2L 


r 


Let’s find the circuit 


analog of some mechanical situations. Consider: 


(3.50) 



fXj : extension of spring 1 , F = -k 1 x 1 . 
]x 2 : extension of spring 2 , F = -k 2 x 2 . 

The same force, F, is exerted on both springs, which implies 

kiXi = k 2 x 2 . 
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Since 


X = X! + X 2 , 

we have 



k,k, 

mx + -—— x = 0 . 

kj + k 2 
= k s 


(3.52) 


(3.53) 


"k s " is the equivalent spring constant of k 2 and k 2 , attached in series. Using the above 
analogy, this equation becomes 


Lq + 



0 , 


(3.54) 
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which represents a circuit with capacitors in parallel: 

L 




C P = 


C x + C 2 


What about springs connected parallel to one another? 



Here we have 


Ft = -kiXi , F 2 = -k 2 x 2 , x x = x 2 = x. 


F = F i + F 2 = ~( k i + k 2 ) x * 


mx 


x + + k 2 ) x = 0. 


= k r 


(3.56) 


Again, substituting from the mechanical/electrical analogy, the above is equivalent to 


Lq + — + — q = 0 , 

yCi C 2 J 


(3.57) 


which represents a circuit with capacitors in series: 
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L 



V„ 



f 





+ 


^2 


Let’s see if we can find the electrical analog of one more case, this time with two masses: 



equilibrium equilibrium 
position for position for 
mass 1 mass 2 


I will derive the equations of motion by deriving the potential energy of the system and 
using Newton’s second law. Since there are two springs in the system, the potential energy 
has the form 

u = J k i x i 2 + \ k 2 x 2 2 • (3-59) 

However, notice that while x 1 and x 2 are the spring displacements, the corresponding mass 
displacements from equilibrium are x x and x : + x 2 . Let us define 

x, = x, , ] 

1 1 • (3.60) 

x 2 = x 1 + x 2 ,J 

leading to 

U = ^ k ix' + ^ k 2 (x 2 - x,) 2 . (3.61) 

The reason for the redefinition is that Newton’s laws refer to particle, not spring, displacements. 
Therefore, we have 
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m i X l 


m 2 x 2 


dU 

dx 1 

dU 

dx 2 


-k 1 x 1 + k 2 (x 2 - x x ) , 
-k 2 (x 2 - Xl ) . 


(3.62) 


These are our coupled differential equations. We will not attempt to solve them but will 
simply use our mechanical/electrical analogy to turn Eqs.(3.62) into 


Mi + Ti - (q 2 - qj 


L q + (q _ q ) = 0 . 

C, 


= 0 , 


(3.63) 


Using Kirchoffs law, this is easily seen to be equivalent to two circuits across a capacitor: 


Li C x 



If we rewrite Eqs.(3.62) in terms of the spring diplacement variables instead of mass 
diplacement variables, we would have found the same circuit but with individual currents 
running across the capacitors rather than the inductors. 

Driven Harmonic Oscillations 

I now wish to discuss driven oscillations. In the electrical case, we will consider a driven 
LRC circuit, 
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C 



V(t) (generator) 


=> Lq + Rq + — q = V(t) . 

C 

We already know the mechanical analog is a driven, damped oscillator: 


mx + bx + kx = F(t) , 


(3.64) 


or 


x + 2|3x + 0 )qX = A(t) , 


(3.65) 


where A( t) 
total force. 


F(t) 


m 


. Notice that "F(t)" is the external force on the particle, not the 


We will talk about two methods to solve Eq.(3.65): 

A. Fourier Series 

B. Green functions 

Method "A" will work for periodic functions, F(t), while method "B" works in the 
periodic or non-periodic case (but is more difficult to carry out). 

Let’s begin by talking about a special case for F(t), which, however, will lead us smoothly 

_ F o 

into a discussion of Fourier series solutions. The equation we will consider is (Ao— where 

F(t)=F 0 cos(cot) ) 

x + 2px + o)qX = A 0 cos (cot) . (3.66) 

Assume a solution of the form, 

x p (t) = D(co 0f a), (3) cos (cat - 8(a) 0 , oo, (3) ), (3.67) 


Download free eBooks at bookboon.com 


107 












MODERN INTRODUCTORY MECHANICS 


LINEAR OSCILLATIONS 


and substitute into Eq.(3.66). After expanding the sine and cosine in the quantities cat 
and 5, and equating the coefficients of sin cat and cos cat, we have the equations, 


cosine: 


D(cao - ca 2 ) cos 6 + 2Dpca sin 6 = A 0 , 


(3.68) 


sine: 


D(oao - ca 2 ) sin 6 - 2D|3ca cos 6=0. 


(3.69) 


From (3.69), 

tan 6 = 

from which we find 

sin 6 = 


2|3ca 


2 2 
(D 0 - CD 


2(3ca 


^/(2(3ca) 2 + (ca 2 - ca 2 ) 2 


cos 6 = 


2 2 
O) 0 - CO 


2pca) 2 + (ca 2 - ca 2 ) 2 


(3.70) 


(3.71) 


(3.72) 


where I have chosen the signs in (3.71) and (3.72) such that 6 >0. From (3.68), (3.71) 
and (3.72), 
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D(CO) = 




(a 2 ) 2 + (2|3to) 2 

Therefore the particular solution can be written very simply as 
x p (t) = D(to) cos (art - 6) , 


6 = tan 


2(3 a) 


2 2 
CD 0 - (0 


(3.73) 


(3.74) 

(3.75) 


Our choice 5 >0 means that we do not use the principle branch of the inverse tangent, 
but instead use: 



Notice that the particular solution, (3.74) is completely periodic; it represents the steady- 
state solution after all initial transients have damped out. We can now add to (3.74) the 
complementary solution to meet a given set of initial conditions. (The complementary 
solution corresponds to one of the three types of solutions for the damped oscillator studied 
previously.) 



values of 5 are in I. 
(We required 5>0 

and we have 5 —* 0 as 

to 0.) 


5(a)) increases as 
o) increases 


Download free eBooks at bookboon.com 


109 











MODERN INTRODUCTORY MECHANICS 


LINEAR OSCILLATIONS 


When w 0 < to => sin 6 > 0 but cos 6 < 0, 



values of 6 are in II. 
since 6(to) must be a 
continuous function of to. 


6(oo) continues to increase 
as go increases 


In any case we see that 

0 < 6 < JT . 


Comparing the cosines of the driving and resulting motion: 

driving: cos(tot), 

motion: cos (wt - 6). 

We see that 6 is the phase difference between the applied force and the steady-state motion. 

A quick and easy way of solving harmonic electrical/mechanical problems is to use complex 
numbers. Let’s redo the solution in this faster manner. Consider, 


mx + bx + kx = F(t) . 

Let’s assume (real parts understood) 

F(t) = F 0 e irat ,1 
x(t) = x 0 e iMt J 

where F 0 is real but x 0 is complex. Substituting these above gives 
[-mto 2 + ibw + k]x 0 = F 0 , 


(3.76) 


(3.77) 
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or ( 0) 


2 

0 



b 

-) 

m 


x o = 


-Lb 


(w 2 - oo 2 ) + (2|3(o) 2 


5 = tan 


J 2(3oa 


2 2 
K (X) Q - 0 ) 


The full solution is thus (again, take the real part) 


x(t) = D(co) , 


the same as (3.74) above. 


(3.78) 

(3.79) 


(3.80) 


An important aspect of this solution is contained in the amplitude, D ( 0 )). We will see a 
graph of it shortly. This graph has a maximum where 


dD(a)) 

da) 


= 0. 


(3.81) 


This is called “amplitude resonance”. It happens when (algebra is a problem) 


(Or — -Ji o^ 2(3 . 


(3.82) 


Thus, unlike the undamped oscillator in Chapter 2 which was seen to have an infinite 
amplitude (as t —» o°)at(0 = (O 0 , the damped oscillator displays a more physical behavior. 


What is the value of D((o) at 0 ) = ( 0 R ? 

A 0 

(o 2 - a) 2 ) 2 + 4(0 2 |3 2 

_ A 0 

^4(3 2 (a) 2 - (3 2 ) P <<a) o 2(3(o 0 ' 



(3.83) 


Thus, for fixed driving force, F(t), when |3«(O 0 , the amplitude at resonance is inversely 
proportional to the damping constant, (3. 
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Now let’s look at D ( 0 )) at the shifted frequencies, co = 0 ) R ± (3 • To the same approximation 

2 

(leaving out terms |3“ or higher), 

® R - <*)„, (3.84) 


- D K ± P) = -r-- 

((">,' - K ± P) ) + 

Ap _ Ap _ 

- V2 (2(3o) 0 ) ' 

Comparing (3.85) and (3.83) shows that 

D (to R + (3) - -^j= D(to R ) , (3« O) 0 . 

The graph we obtain for |3« u) 0 is then: 


(3.85) 


(3.86) 



Thus when |3« o) 0 , the width of the amplitude response graph is given approximately by 
2(3. This is called the “full width” of the resonance response curve. 


An important figure of merit for such oscillating systems is the so-called Q (quality) factor, 
defined as 


Qa 


on 


(Aco) 


Full 


(3.87) 
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where (Aoo) Full is the full width. When Q » 1 the resonance curve displays a very high, 
narrow peak. In a problem we will see that Q A is proportional to the average total energy 
divided by the energy lost during a cycle of driven oscillation, when oo = 0 ) R and Qa >:> 1 • 

An additional type of resonance is called “velocity resonance”. The velocity response in the 
steady state is 

x p (t) = -toD(co) sin (cot - 6), ( 3 . 88 ) 

from which we define 


D v (to) = coD(co) , 


and 


dD v (co) 


dco 


co = co v 


0 . 


(3.89) 


(3.90) 
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One can show that 


W v W 0 


(3.91) 


so that velocity resonance occurs at the undamped angular frequency. In addition, the equality 

oa„ 


Q„ - 


(Acd) f 


2(3 


(3.92) 


CD r 


is exact in this situation. (One only has Q A - —- , Q A >> 1 , for amplitude resonance.) 

2(3 

When discussing driven motion, it is often convenient to introduce time-averaged quantities. 
We define 


c ••• > = - J dt (•••) 


2 JI 


oa 

2tc 




(3.93) 


Therefore, 


f il UJ 

sin 2 (cot - 6)) = — / dt sin 2 (cot - 6) , 
2 JT 


2jt:-5 


-I —u 

-f dx si: 

2jt _ 6 


.2 1 [X 1 . „ 

in x = —-— sin 2x 


2it \ 2 4 


2jt-S 


-6 


l 

2jt 


2n - 6 


- i sin (4jt - 26) - f — - i sin (-26) 


— . (independent of 5) 


Similarly for cos 2 (oat - 5) . Therefore, we have for the kinetic and potential energies, 


(t) = ^ moa 2 D 2 (oa) , 

(u) = ^ moa 2 D 2 (oa) . 


(3.94) 
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Notice from the above that (t) » (u) for w » (D 0 , so that most of the energy is kinetic. 
Anyone who has pumped their feet too quickly on a swing will understand this statement! 
Also notice that since (t) °c D 2 (at) , the kinetic energy resonates at 0) = (D v (velocity 
resonance) while potential energy resonates at o> = cd r (amplitude resonance). 

Fourier Series Methods 

Our solution, Eqs.(3.74) and (3.75), to the differential equation, (3.66), now affords a 
general method of solving (3.65) when the driving force, F ( t ), is periodic. Our differential 
equation can be written in operator form as 

ff x(t) = A(t) , (3.95) 


where 


d d , 

ff = —- + 2(5 — + on . 

dt 2 dt 0 


A linear differential operator is such that 
ff(x 1 + x 2 ) = ffx 1 + ffx 2 . 


(3.96) 


(3.97) 


Eq.(3.97) holds for any derivative operator of any order, which can be proven by induction. 
Thus, (3.96) is a linear operator. This would not be the case, for example, for a damping 
term proportional to x 2 , since 

2 


( d ( x l + x 2 )l 

I 2 ( dx A 

i 2 1 ^ dX2 "i 

dt ] 

i * ra 

1 + 


The linearity of (3.96) now implies that the steady state solution to ((() is an arbitrary phase) 


F(t) = ^ a n cos (a) n t - (|)) , 


(3.98) 


is just 


X r M = t: 2 

m 


ol 


= tan 


n V (®0 - w n) + 4 co^(3 2 
- i ( 2(D n (5 


cos (co n t - <t» - 5 n ), 


(3.99) 

(3.100) 


The reason this is important is because of the Fourier series theorem. Any periodic function 
F(t) = F(t + T) (subject to some nonrestrictive continuity conditions) can be represented 
by the infinite series, 
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F(t) = — + ^ (a n cos (neat) + b n sin (neat)), 


( 3 . 101 ) 


where 


2 ? , , . ( 2 iint ' 

a n = — | dt F(t ) cos - 

x _t \ x 

2 


2 ? 


b n = - | dt 1 F(t ' ) sin 


( 2jmt 


X _x 
2 


\ T 


( 3 . 102 ) 


One can actually integrate over any time interval of length X in doing the integrals (3.102). 
The particular solution to (3.65) is then 


x p(t) = 


where 


1 °° 

+ — ^ D (not) (a n cos (ncot - 6 (nco)) + b n sin (neat - 5 (noa))) ,(3.103) 


2ro 0 m m n=1 


D (noa) = 


(w 2 0 - n 2 oa 2 ) + 4n 2 oa 2 (3 2 


6(noa) = tan 1 


2noa(3 


2 2 2 * 
(o 0 - n to 


(3.104) 

(3.105) 


Let us do a Fourier series decomposition as an example. 


Example 1 

F(t) 


F 0 , 0 < t < | 

—F 0 , - ^ < t < 0 . 



F(t) 

F 0 

T 

2 


X 

2 

- 




0 

-F 

r 0 
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Solution: 

Since this is an odd function of t, we have 

a n — Of 


for all n. For b n we have 


T 


2F n , . / 2nnt 

dt sin - 


\ x 


2F 0 } ^ , . / 2jtnt 

+ —5- I dt sin - 

To \ T 


b n = 


2F„ 


jtn 


[■ 0 , n even 


(l - cos (Jtn)) = <{ 

l 


4F, 


0 ,n odd . 


jtn 


The explicit solution to the damped oscillator equation with this force is 


x p(t) = 


4 Fa 


D(nu)) 


Jtm „ 1 n 

n = 1 


sin(noot - 6(nco)), 


where 6(nco) is given by (3.105) as usual. 


Green Function Methods 

As mentioned above, there is another, more general method of solving the driven mechanical 
circuit equation (3.65). We will start the explanation of the method by asking what happens 
when we “kick” an underdamped oscillator. Qualitatively, we get: 



T 

We have “kicked” it at a time beginning at t --. The above is in response to an 

idealized force, 
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F T (t) 

F„,_ 


2 t 


(3.106) 

V 

Let’s idealize this situation further to a very short, sharp kick so that the resulting motion 
is independent of the details of F T (t). Define 

f = J dt F(t), (3.107) 

which we will call the impulse. It is momentum transfer. Clearly, from our form (3.106) 
we have 


This is given by 


F T (t) = 


T 

2 0 


T 

t <- 

2 

T T 

- — < t < — 

2 2 

T 

t > — . 

2 


I 


F 0 T. 


(3.108) 


What happens to the oscillator during this time, T? Under the conditions 
1 2jt 

T « — and T « — , the oscillator is approximately “free”. So, we simply have 

(3 O) 0 


mx - F 0 , 

T 

x - — | F 0 dt 

m _t 
2 


F T 

D- 1 - 

m 


=> 


_ 1 F 0 2 

x(T)-- T 

2 m 


I 

r 

m 


1 I 

-T . 

2 m 


(3.109) 


Therefore, in the limit T —» 0 + ( "0+" means to approach zero from positive numbers), 
I = const. , our impulse function serves only to give the oscillator an initial velocity, 
given by (3.109). This limit produces an increasingly narrow force profile: 
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Let’s define the “kicking function”, 


3(t) = lim F T (t). 

T—»0 + 

I=const. 

We have ( a < 0 < b) 

T 

b b 2 

f ^5(t)dt = f lim F 0 dt = lim J F 0 dt= I, 

** ^ T-»0 + t—>0+ t 

I=const. I = const. - — 

2 


(3.110) 


(3.111) 
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where I have assumed the limit and the integral can be interchanged. Also notice for an 
arbitrary but continuous function A (t), 


f 3(t) A(t)dt = lim F 0 f A(t) dt = I A(0) , 


( 3 . 112 ) 


T—»0 + *p 

I = const. - — 


- A(0)T 

where I have made the same assumption. In addition, 


S(t) = 3(-t) , 


(3.113) 


by construction. 


Now let us say we kick the damped oscillator at time t = t ' instead of at t = 0 . Here 
is the original situation: 


S(t" ) 


0 (origin) t' 1 

Here is the same situation described with a shifted time axis: 


S(t-t') , 


0 (new origin) 


Therefore, t - t ' is the appropriate shifted argument of the kicking function, Sf. 
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Let us solve 


x + 2|3x + 0 )qX = — S(t - t ' ) , 


m 


(3.114) 


subject to the initial conditions, 


b. 


x(t') = 0, 
lim x(t) = 


m 


(" lim " means t approaches t' from the positive side.) The second condition is just 

t —* t ' + 

supplying the initial velocity condition, Eq.(3.109). Whenever t ^ t ' the right hand side 
of (3.114) is zero and so is easy to solve. (Let us assume (3 < oo 0 ) : 


x(t) = 


A ' e pt sin (ccqt + a ') , t > t 
0 , t < t ' 


(3.115) 


Notice the two undetermined constants, A ' and a ' . For convenience, let’s take 
a ' = —CDjt '+ a , 

A' = Ae pt , 

since t' is a constant as far as our differential equation is concerned. Using these redefined 
parameters in (3.115) becomes 

\ Ae“ p<t “ t 1 sin (oadt - t’ ) + a) , t > t’ 
x(t) = \ (3.116) 

[ 0 , t < t'. 

Now apply initial condition "a" above when t=t' : 

=> A sin (a) = 0 , 

which means we can choose a = 0 . The remaining parameter, A, is determined by condition 

"b" . To see how this condition comes about, let’s integrate both sides of the differential 

T 

equation, (3.114), over a short time interval, 2e (e > ^ ) , centered about t=t' : 

1 


lim f 

F. —* 0 4- J 


x + 2px + o) 0 x = — is(t - t ' ) dt 
m 

lim (x (t '+ e) - x(t e) ) = — . 

e-o+ v-v--' m 


v o 


o 


Using the explicit t > t ' solution then shows 


lim x(t '+ e) = AO), => A = 

E^0 + 


men. 


Download free eBooks at bookboon.com 


121 








MODERN INTRODUCTORY MECHANICS 


LINEAR OSCILLATIONS 


Everything is now determined. Using a new notation for our solution, x(t) —■> x(t, t 1 ), 
we have 


x(t,f ) 


I -P(t-f) . r /. 

- e sinkodt 

mo) 1 


0 


t')] , t > t' 


r 


t < t' . 


(3.117) 


Tltis solution looks like we thought it would (assuming 


I > 0 so that v 0 > 0 ) : 



Notice that x(t,t') is continuous across t = t' (condition "b") but that it’s first 
derivative is not. The velocity profile is 



In solving (3.114) we have done much more than solve a rather isolated, special example. 
Just like the solution we formed for the harmonic force, we can now build on this solution 
to solve a wide variety of other problems. The comparison with the Fourier series solution is 
quite enlightening, and I present analogous quantities in the two situations in the following 
Table: 
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Fourier Series 


Green Function 


Initial condition 
of oscillator 


periodic motion 
(x(t + T) = x(t)) 


quiescent 

(x(t,t' ) = 0, t< t' ) 


"primitive" force 
"primitive" 
solution 
Nature of F(t) 

decomposition 

F(t) representation 

General solution 
for x(t) 


F 0 cos (tot - 0) 

D((0) cos (COt - <j) - 8) x (G t 1 ) | t’ fixed 


discrete 

^ a n cos(ntOt - <|>) 
n 


continuous 


dt’ 


F(t') 


(t-t 1 ) 


^a n D(nto)cos(ntot —(|)-5(nco)) ? 

n 


In either case we are decomposing the given F (t) into either a discrete series of harmonic 
functions (Fourier) or a continuous density of “kicks” (Green function) and using linearity 
and the known “primitive” solutions to get the response, x(t). I have deliberately left a 
blank in the table where the Green function general solution is supposed to be. We can guess 
the form of the solution based upon the analogy. The table informs us that the following 
quantities or operations are analogous: 

00 

2 ^ 1 ' 

n —oo 


a - 


D (not) cos (neat - 4> - 5 (noa)j <=> x(t, t ' ) . 


Based upon this, we will hypothesize the form of the general solution using x (t, t' ) : 

, (F(t ' ) 


x(t) = J dt 


x(t, t ' ) . 


(3.118) 


F( t' ) 

That is, the full solution is just given by integrating the force per impulse, -1- L over 

I 

the solutions, x (t, t' ) For convenience, let us introduce (finally, the “Green function”!) 

x(t, t' ) 


G(t,t' ) = 


(3.119) 


Since G(t,t' ) — 0 for t < t' . we can write the general solution as 
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X(t) = J dt' G( t, t' ) F(t' ), 


(3.120) 


Let us verify the correctness of (3.120) by substituting it back into (3.65). First, to make 
stonger contact with standard terminology, let us introduce 


6(t - t ' ) = — £s(t - t ' ) , 

I 

called the Dirac delta function. We now have solved 


( 3 . 121 ) 


dt 


+ 2p — + (Oq I G(t, t ' ) = — 6 (t - t ') , 

dt / m 


where the delta function has the property (A ( t ) is some arbitrary smooth function and 

a < 0 < b) 


f 5(t) A(t) dt = 


A(0) , a < 0 < b 
0 , otherwise . 


To verify we have a solution, putting x(t) = J*dt' G(t, t' ) F(t' ) into (3.65) gives 
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, 1 

X + 2|3x + O) 0 X = — F(t) , 
m 

- F(t), 
m 


We can now automate the solution to (3.65): 

(1). Solve the linear differential equation, 


f 


dt' 


d d , 

— + 2P — + to 2 0 


dt 


dt 


G(t, t' )F(t' ) = 


= —8(t- t') 
m 


- F(t) = - F(t). 
m m 


ft G(t, t ' ) = — 6(t - t ' ), 
m 

for some second-order linear differential operator, 

a. G(t,t' ) = 0, t <t' , 

b. lim G(t, t ' ) = — . 

t^f+ m 



( 3 . 122 ) 


. . . , subject to the conditions, 


(2) , Do the integral 

t 

x(t) = Jdt' G(t, f ) F(t' ). (3.123) 

—oo 

(3) . One can add the complementary solution to reach some desired initial conditions. 

Although our “kicked” solution, x (t, t ' ) represents an initially quiescent oscillator, this is not 
necessarily true for x ( t ) in (3.123) because we can let F(t' ) be defined all the way to t 1 —> -oo. 
If we do this with a periodic F(t ' ), the resulting x(t) from (3.123) will correspond to 
the particular solution. In addition, in step (3) one can add to (3.123) a complementary 
solution which will then match any desired initial conditions, quiescent or not. Thus, the 
Green function method works in essentially all cases; the price to be paid for this generality 
is a sometimes difficult integral to perform. 

Let us do some examples to illustrate this method. 

Example 2 

A driven undamped oscillator’s equation of motion is given by ( 0 ) * O) 0 ) 
mx + kx = F 0 sin oat . 
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Find x (t) for an initially quiescent oscillator at t = 0 : 
x(0) = 0 , x(0) = 0 . 

Solution: 

This is actually Example 1 from Ch.2. The appropriate Green funtion is just the p 
limit of the G(t,t') in (3.117). We will use 


G(t,t 1 ) = 


—sin [oo 0 (t - t')],t > t' 
0 , t < t' , 


moo r 


Since we are not interested in the solution before t = 0 , we may take the force profile to be: 



Thus 


x(t) = J dt 1 G(t, t 1 ) F(t' ) 


moo 


i u 

— J* dt ' sin u) 0 (t - t ' ) sin oat ' . 


0 0 


We can write this as 


x(t) = —— <Jsin (OotJ dt' cos oa 0 t' sin oat' 


moa r 


- cos oa n tJ dt' sin oa n t' sin oat'! 


Some useful integrals are, a * b: 
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J 

f 

f 


sin ax sin bx dx 


sin ax cos bx dx 


cos ax cos bx dx 


sin (a - b)x sin (a + b)x 
2(a - b) 2(a + b) 

cos (a + b)x cos (a - b)x 
2(a + b) 2(a - b) 

sin(a + b)x sin(a - b)x 
2(a + b) 2(a - b) 


a = b case: 


/ sin 2 ax dx = — (x - sin 2ax), 


f 3i 


in ax cos ax 


dx = sin 2 ax. 


J cos 2 ax dx = ^ (x + sin 2ax) . 


Using these above and doing the algebra gives 


x(t) = S- 


m to? - to 2 


to 

sin tot - — sin co 0 t 
C0 n 


This is the same as we had previously. We are “cracking a walnut with a sledgehammer.” 


Example 3 

Redo the harmonically driven, damped oscillator using the Green function technique. That 
is, solve the differential equation 

mx + bx + kx = F 0 cos (tot) , 


where the right hand side defines the force for all t, -oo < t < °° . 

Solution: 

Since the force has been defined for t —^ -oo , we should recover the particular solution, 
Eqs.(3.74) and (3.75) above. We have 

t 

x(t) = J dt ' G(t, t ' ) F 0 cos (tot ' ) 

—00 

t 

= f dt' e“ l>(t_t 1 sinfoodt - t' )] cos tot'. 

moo, J 

i —00 
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This integral looks nasty. However, notice 

cos cot' sinfw^t - t' )] = |sin[(ot' + ot^t - t' )] 

- sin[art' - to^t - t' )]j , 

so that 

x(t) = —^— | J dt' e“ |!<t “ t ’ sin [(a) - wjt' +0)^] 


- jf dt' e P(t 11 1 sin [(a> + cojt o^t] 


One can now do the two integrals above with the use of 

Q aX 

f e ax sin xdx = —-(a sin x - cos x). 

J a 2 + 1 

After changing variables, the result is 


x(t) = 


2m(0, 


(a) - 0)J 


P 

CO — CO, 


(0-0), 


sin art - cos art 


+ 1 


P 


(a) + atj / j 3 \ 2 \a> + a) 


sin art - cos art 


+ 1 


a) + a>. 


x(t) = A- 
2m 


4co(3 


(a) 2 - O) 2 ) + 4|3 


sin a)t 


22 

0) 


2 |a> 2 - a) 2 ) 


(a) 2 - a> 2 j +4(3 


cos art 


22 

0) 


Download free eBooks at bookboon.com 


128 



















MODERN INTRODUCTORY MECHANICS 


LINEAR OSCILLATIONS 


This can be written 


17 


\ 



2oop 

- oo 2 ) + 4(3 2 oo 2 
= sin 6 


sin art 


x(t) = 


- <■>;)’ 


+ 4P 2 o) 2 


cos (cat - 6) . 


This is precisely our previous solution to the problem; the phase angle 5 having the same 
meaning as before. 
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CHAPTER 3 PROBLEMS 

1. Substituting Eqs.(3.17) into (3.16), show that Eq.(3.19) (an ellipse) results, 
where C and D are given by Eqs.(3.20) and (3.21), and tan 2y is given 
by Eq.(3.18). [Hints: After substituting (3.17) into (3.16), require that the 
coefficient of the x'y ' cross term vanish. After this, use (3.18) to eliminate the 
cos 6 term and do the algebra leading to Eqs.(3.20) and (3.21).] 

2. Show that the coefficients C and D in Eqs.(3.20) and (3.21) are both positive. 
Do this as follows: 


a) Show that requiring C>0 is the same as 


B + A > 


- B 


cos 2y 

b) Similarly, show that D > 0 implies 


B + A > 


B - A' 


cos 2y 

c) Show that (a) and (b) together imply the true statement, 
cos 2 6 < 1, 


when the explicit form for cos 2y is used. [The other three cases, C>0 D<0, C<0 D>0, 
and C<0 D<0 may be eliminated one by one.] 


3. For the out-of-phase two dimensional oscillator, study the special case where 

A=B and show that in this case we can choose 

1 - cos 6 1 + cos 6 

C — —^, D — 


A 2 sin 2 5 

71 


A 2 sin 2 6 ' 


and that y = ± —. [Note: It is not correct to simply put A=B in Eqs.(3.20) and (3.21)] 

4. (a) The general solution of the damped oscillator, 
x + 2|3x + o) 0 2 x = 0, 

when (3 > O) 0 (overdamped case) is given by 
x(t) = e - ! 31 A 1 e“ 2t + A 2 e _ 2 

where W 2 = ^/p 2 - 0) 2 and Ai f 2 are constants. For initial conditions 


x ( 0 ) = x 0 , x(0) = v 0 . 
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find Ai and A 2 . 

b) Find conditions on Xo , Vo which permit a change in sign of x(t). 

5. In the analogy developed between mechanical oscillators and electrical circuits we 
used KirchofFs circuit law, which says the sum of voltages around a closed loop 
is zero. Which of Newton’s laws of motion is analogous to this circuit law? Can 
you formulate circuit statements analagous to the other two of Newton’s laws? 

6. A mass is attached to a wall with the three springs shown. The system is without 
friction. 



Using the circuit/mechanical analogy, find the equivalent circuit to this mechanical 
oscillator. 


7. Consider the system with two masses: 



Xi f 2 locate the positions of masses mi f 2 relative to mi’s equilibrium position, Xo. 
Assume there is a frictional force proportional to velocity between mi and m 2 , but that 
there is no frictional force between mi and the ground and that the oscillations are small. 

a) Write down Newton’s equations of motion for this system. 

b) Find the electrical circuit analogous to the mechanical situation in (a). 
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8. Consider the system of two masses shown. One mass rides on top of the other 
and the two springs have different spring constants, and k 2 - There is a 
friction term linear in velocity between mi and m2 but not between rrq and the 
ground; the oscillations are small. 


/ 

/ 

/ 

/ 


k l 

m 2 

k 2 

— mrv|~~ 




/ 

/ 

/ 

/ 


a) Write out Newton’s equations of motion for this system. 

b) Find and draw the equivalent electrical circuit to this mechanical problem. 

9. Consider the system of two masses and three springs shown. 


/ 

/ 

/ 

/ 


~nmr l 1 


"lim 


'Tmr 


/ 

/ 

/ 

/ 


a) Write out Newton’s equations of motion for this system. 

b) Find and draw the equivalent electrical circuit to this mechanical problem. 

10. Consider the driven, damped harmonic oscillator: 
x + 2px + o) 0 2 x = A 0 cos art. 

a) Show that the instantaneous energy of the oscillator 
(E(t) = — mx 2 + —mo) 0 2 x 2 ) satisfies 

dE(t) , 2 . 

—-— = -2fimx + mAo x cos oat. 
dt 

[Hint: Take the derivative of E (t ) and use the above differential equation.] 

b) Give an interpretation of the two terms on the right. 
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11. a) Show that “amplitude resonance” occurs at angular driving frequency 

(0 R = - 2|3 2 . 

b) Show that velocity resonance takes place at the system’s natural angular 
frequency, 

Q) v = Wo. 

c) Find the angular frequency of accleration resonance. That is, find the driving 
frequency, w A , such that the acceleration amplitude of the oscillator is a 
maximum. 

12. Show that the velocity quality factor, defined by 

CDv 

Qv " 2P 

(w = co 0 from problem 11 (b)) can also be written exactly as 
WO 

Qv " Aw ' 

where Aw = full width of the velocity resonance curve at —j= X maximum velocity. 

V2 

13. Can a driven, overdamped oscillator undergo amplitude resonance? What about 
velocity resonance? Explain why or why not. 

14. Show that the Q A factor for light damping for the driven oscillator may be 
written as 

Q A - ftN e , 

where N e is the number of oscillations of the free oscillator which occur in the time it 
takes for the amplitude of the undamped oscillator to decrease to ^ of it’s initial value. 

15. (a) Show that 

< E > w + w 0 
AE 8 jt|3w ' 

where <E> is the time averaged energy of a driven oscillator and AE is the energy lost 
during a cycle of driven oscillation. 

b) Using (a), prove that for light damping the amplitude quality factor, Q A , is thus 
given by 
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Qa - 


2n 


< E > 
AE 


f 


16. Referring back to Prob.10 above, explicitly compute the time average of both 
terms in the energy change of the damped oscillator over 1 cycle, < dE ^ >. 
Does the answer you get surprise you? 


17. Obtain the Fourier series representing the periodic function, 


F(t) 


It looks like: 


0, - A 

,2jt 

cos (— t) , 
A 


< t < 0, 

0 < t < A. 
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"A" is an arbitrary constant. 

18.A damped, driven harmonic oscillator with natural frequency Wo and damping 
constant p is subject to the force (for all times times, t), 

F(t) = a : cos (cat) + a 2 cos(2oat). 


where a x and a 2 are constants. Find the average kinetic, <T>> and potential, <U>, 
energies of this system. [Some of the integrals given in Example 2 above may be useful.] 


19. A underdamped harmonic oscillator (damping constant b and natural angular 
frequency Wq) is subject to a force as a function of time given by 


F(t) = 


F 0 cos cat, t < 0 
0, t > 0 


This looks like: 



Find the motion of the oscillator, x ( t ), for t>0 . [Flint: There is a simple way of doing 
this problem considering that you should know x ( t ), for t<0 . ] 
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20. An initially quiescent underdamped harmonic oscillator, with damping constant 
p amd natural frequency Wo, is subject to an infinite, periodic square wave of 
amplitude A and period T, starting at time t=0. It’s differential equation is: 

d 2 d 7 1 

(dt2 + 2 Pdt + w o ) x (t) = m F (t)- 



Find the response of the system, x (t), for all times t. [Hint: Use Example 1 above, 
making sure the initial conditions x(0) =0, x(0) = 0 are met at t=0+.] 


21. a) Solve (Green’s function for a free particle) 


subject to: 1. G (t, t' ) = 0, t < t' , 


dG(t,t' ) 
dt 


1 

m 


for t 


t' + . 


b) Using 

x(t) 


f dt ' G(t, t ' )F(t ' ) , 


find x (t) for 


F (t' ) = 


F 0 , f > 0 
0, t' < 0 


where 


Fq is a constant. The result should be very familiar! 
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22. Solve the harmonic oscillator equation, 


mx + kx 


0, t < 0 

F 0 sin(o)t) , t > 0 


for an initially quiescent oscillator using the Greens function technique for the case 
2 k , 

00 = 0)0 ( w 0 = — ) . Use the known Green’s function for this problem from the text. 

m 

[Hint: The answer for x ( t ) is given somewhere in the text.] 


23.An initially quiescent, undamped oscillator with mass m and natural angular 
frequency oo 0 is subjected to a force, 


F(t) 


0, t < -T, 

—F 0 , - T < t < 0, 

F 0 , 0 < t < T, 

0, t > T. 


This looks like: 


F(t) 


-T 


0 


T 


t 


-F n 


a) Using the Green’s function method, find x(t) for all times, t. [Hint: 

Consider the regions -T<t<0, 0<t<T and t>T separately.] 

2 jcn 

b) Show that x (t) = 0 results for t > T when oo 0 = , n = 

12 3 

-1 L f j m • • • 


24. Consider the following force profile on an undamped (damping constant 
P = 0 ) oscillator of natural angular frequeny, co 0 . The driving angular 
frequency is ca ^ (jOq • 
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71 


The force is zero for t > and for <- _ —. Find the response of the system 


0) 


0) 


71 


71 


1 o 2o) 


x (t) for t > . [Ans: For t > , I get x (t) = - 

' ' a) a) ' 5 ' ' moo n (D^-w 


2 -m2 


Jt TC 

sin( ) cos (oo n t).] Extra Credit: Also find x(t) for - < t < 

' 0) v o ’ 0) 


0) 


25. The Green function for a damped oscillator satisfies 

+ 2 PcE + “o )G(t,t 1 ) = 

where 6(t-t' ) is the Dirac delta-function. Solve this equation and construct the 
appropriate Green function, G ( t , t ' ) , for a critically damped oscillator, O) 0 = [> 
(The solution in the notes and in the book is for the underdamped case, o) 0 > (5.) 

26. One can use the Green function to also build in arbitrary initial conditions for 
an oscillator. Newton’s equation for an undamped oscillator is, 

d 2 2 1 

(dt 2 + “0 )x(t) = “ f (t) , 

m 

where f (t) is the external force. It’s Green function is 


G(t, t ' ) = 


— sin [oo 0 (t - t ')] , t > t ' 
0 , t < t ' . 


mo) 


For an oscillator in an arbitrary initial state at t=0, show that x(t) is given by 

(Xq and Vo are the initial position and velocity, respectively, at t=0), 

" / dG(t,t' ) 

x(t) = f dt'G(t,t')f(t') + m Ig( t, 0 ) v 0 - x 0 -- 

o 

[Fiint: There are special techniques to show these sorts of results. For our purposes, 
it’s enough to show that x (t) is a solution of the above differential equation with 
the correct initial position and velocity, x 0 and v 0 . One must aso show that 
G (t, t' ) = G(-t' , -t) .] 

Download free eBooks at bookboon.com 



138 












MODERN INTRODUCTORY MECHANICS 


NONLINEAR OSCILLATIONS 


4 NONLINEAR OSCILLATIONS 


THE ANHARMONIC OSCILLATOR 

Let us go back to how we started Chapter 3. 


+x 


Equilibrium 0 

position _►. 

(origin) 


/////// 


FI 


F.„(x) 


-mg 
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We have 


mx = F tot = F sp (x) - mg. 


F sp (x) = F (0) + x 


dF 


sp 


x 2 ( d 2 F 


sp ' 

mg 


ax ) 0 2 \ d *‘} 0 


= -k = 2 e 

So, a more accurate description of the motion for a general spring is 


mx = -kx + ex 2 ; “anharmonic oscillator” 


(4.1) 


(4.2) 


this is a non-linear differential equation. It is still second order (requiring two initial 
conditions) and still conservative. 


Even though (4.2) is non-linear, we can calculate t(x) for portions of the motion not separated 
by turning points as before. The potential energy function can be taken as 


U(x) 



and looks like (assuming £ > 0) 


(4.3) 



Thus, the potential is no longer symmetric about the equilibrium position. In fact, the magnitude 
of the force is slightly greater when it compressed (x < 0) than when it is stretched. 


The formal solution to t(x) is given by Eq.(2.42): 

dx ' 


x(t) 


(t -1 0 ). ±f 


- ( E - u < x ' >) 
m 


(4.4) 
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with U(x' ) given by Eq.(4.3). Although we have a formal solution, turning t(x) into 
x(t), which much more useful, is quite difficult, except in special cases. Even finding 
closed form expressions for the turning points, x 3 , x 2 , is awkward since one is examining 
the cubic equation, 


U(x) - E = 0, 


(4.5) 


=> 


3 


X 




= 0 , 


(4.6) 


for the real roots x 3 , x 2 and x 3 


We can relate the anharmonic oscillator equation of motion to a special function called 
an elliptic integral. If we agree to choose t 0 = 0 when the spring is released at rest after 
compression, we will have. 



where is the negative root, which is a turning point of (4.6), assuming £ > 0. The 
three roots of (4.6) can be pictured as the positions where the E=const. line cuts the 
potential energy: 


U( x) 



The integral in (4.7) can be related to an elliptic integral of the first kind: 

c dx 

■ -f V(1 - x’Xi - kV) • 

Using the integral, 


(4.8) 
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“_dx^_ 

{ V( a - x'Xb - X')(x'-c) 

where a > b ^ u > c and 



In our case, we have 

a = x 3 , b = x 2 , c = x lf u = x(t). 


Va - 


f(y, q) 


(4.9) 


(4.10) 


(4.11) 


which gives 

t = 2 


/ 3m' 

r_ 

F (Yr q) 

Us, 


xl 

CO 

1 

X 

ted in the period ol 

r two: 

/3m\ 

1/2 i 

F (!' 9) 

UJ 

v 

X 

co 

1 

X 


(4.12) 


(4.13) 


Sometimes the expression F ( 1, q ) is called a “complete elliptic integral of the first kind.” 


Exact expressions for the roots of (4.6) are available; however, to first nontrivial order in £ 
, one can show that 


x i ^ 


x, - 


2E 
k 

2E 
k 

3 /k 

2 L 


1 

( 2E 

f 2 


1 - — 
3 

(i? 

E 

i 

\ 




1 


y% \ 


1 + 3 1 

V 


/ 

t 

(. 8 

E \ 

2 ^| 


( si 


'I' 



(4.14) 
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Using (4.13) and (4.14) and an approximate form for F(l,q) (see Eq.(4.27) below and 
the HW problem), we can find an approximate expression for x for small amplitude motion, 
which will be a power series in £. I find that 4 


X 






= x„ 


(4.15) 


so that the period is slightly in excess of x 0 , for fixed E and K. 

The Plane Pendulum 

Let us examine another intrinsically nonlinear system, the plane pendulum. 



Assume planar motion of the mass m and take the origin of coordinates (cylindrical) to be 
at 0, the point of attachment, (z-axis is out of the page.) Note the attachment is such that 
the pendulum can swing a full 360° and go “over the top”. 

For angular variables, the appropriate concept is torque: 

L = N , 

N = r x F = mg^ sin (fi (-e z ) , 

L = r x p = mr x v . 
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In cylindrical coordinates, 

remember 


v = pi p 

+ p4>e^ + 

ze , 

V z 

0 


0 

=> L = (m pe p ) 

X 

-©-• 
CD > 

-e- 

= m p 2 ^S z , 

=> L = m i 2 § 

e z (p = 

£ here). 

The equation of motion is 

thus 



ntf 2 4> = -mg^ sin 4> . (4.16) 

For small oscillations sin cj) ~ (() and we have harmonic motion 

$ + - <l> - 0 , (4.17) 

t 

2 g 

with angular frequency to 0 = — . 
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Going beyond this approximation, we get from (4.16) the energy conservation equation 
by integration 

=> — ntf 2 (j) 2 - mg£ cos d) = const. (4 18) 

2 

We can always add or subtract a constant in defining the total energy, E. We will choose 
const. -» E - mg^ , 


which corresponds to choosing the zero of potential energy at the point (j) = 0 . We have 

(4.19) 


E = — m/ 2 (j) 2 + mg^ (l - cos (j>) . 


= T 


= U(<J)) 


The periodic potential looks like: 



The diagram tells us that for energy values such that E < 2mg/ ( or (j) 0 < K) the 
pendulum motion is bounded in 4>. this would correspond to an energy E 2 in the above 
plot. We can get the relation between t and <j) by again relating the motion integral to an 
elliptic integral. From (4.18) we have 

* - * (^) ( E - U «M)‘ /2 ' (4.20) 


Choosing t 0 = 0 at (|) = 0 now gives 


t = 


\ V2 V d(J) ' 


J 


\l/2 


, 2 ; 0 (e - U ((j) ')) 

where we can integrate only up to the first turning point, 4> = 
(() = §o > the energy, E , is given by 

E = mg^ (l - cos 4> 0 ) = 2mg^ sin 2 . 


(4.21) 


At maximum extension, 


(4.22) 
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Also, we can write 


U(4>) = 2mqt sin 2 — , 
2 


which now gives 


(4.23) 


, 1/2 




d(t> 


J, 4 

0 I • 2 TO • 2 Y 

1 sin - sin — 


1/2 


(4.24) 


Defining 


x s 


• <t> 

sin — 
2 


sin 


4> 0 


4*0 

sin — r 
2 


k = sin — 


we 


have 


dx = 


cos — d(|> 

2 sin — 


V1 - k 2 x 2 
2k 


dcfi , 


, !/2 k sin - 


-ii 


J 


2kdx 


(1 - kVf (*■ - kV)‘' ! 


(e) 

1/2 

| „ / 


1 F ( 


• <t> 

in — 
2 


(4.25) 


where we again see the elliptic integral of the first kind. The period, T, is now given by 
setting 4> = <j) 0 in (4.25) and multiplying by four, 

/ \ 1/2 

x = 4 I c ~\ F(l, k) , (4.26) 

We may get an approximate expression for t for small oscillations, k « 1 , by expanding 
F(l, k): 


F(l, k) - f . , = = , 

d Vl - x i/l - k x 

= f dX [l + - k 2 x 2 + - (k 2 x 2 ) 2 

j Vi - X 2 V 2 8 ' J 


(4.27) 


Download free eBooks at bookboon.com 


146 















MODERN INTRODUCTORY MECHANICS 


NONLINEAR OSCILLATIONS 


Doing the integrals in (4.27), we find 


T = 4 


= 2jt 


Also 


(e) 

1/2 

( 71 


9jt 4 \ 



- 

+ — k + 

-k + . . . 

UJ 

l 

\2 

8 

128 ) 


(L) 

1/2 

1 ( 

u 



1 + — k 2 + k 4 + . . 
4 64 


= 


k = sin — - — - — + 


sin 

2 2 48 

k 2 1 + 

4 48 


k *. « + .. 

16 


Collecting terms, and to order (|) 0 , we thus find 

4> 0 2 


T = T„ 1 + —— + 


11 ad ^ 

(Po "I" • • • | 

16 3072 j 


Again, the period is slightly increased. 


(4.28) 


One can understand why the plane pendulum’s period is nearly independent of amplitude 
from the smallness of the (|) 0 2 , 4>o coefficients. For example, choosing (j) o = — we get (from 

a numerical evaluation of F (1, —= ] ) 

l V2j 

x ~ x 0 (1.1803) , 
or about an 18% change. 


Phase Diagrams and Nonlinear Oscillations 

An important qualitative technique for the understanding of both linear and nonlinear 
systems is the construction of phase diagrams. Let us return to the harmonic oscillator. 
It’s energy is 


E 



(4.29) 


Plotting x(t) and x(t) simultaneously gives a series of ellipses for various energies: 
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This is called “phase space.” The arrows show the evolution of the system in time. 

Of course, the trajectories will never intersect for E ' * E This is true for any linear, second 
order, system. This can be seen by the following argument. Assume such a system had a 
phase diagram where the trajectories crossed, as in 


x(t) 



Imagine starting a particle at the above intersection point. If the above could happen, 
this would mean the particle at the intersection has two possible solutions. However, the 
solutions of linear, second order, differential equations are unique, given initial conditions 
on x and x. We will see shortly that a special type of intersection is actually allowed for 
nonlinear systems. 

Another type of behavior not allowed for linear, second order systems in phase space is a 
bifurcation, as in the following. 
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Actually, it is clear that no deterministic system can display such behavior. Although 
bifurcations in phase space do not occur, bifurcations in parameter space descriptions can 
and do occur, as we will also see shortly 

Since the underdamped oscillator looses energy as time progresses, it’s phase diagram will 
look like: 
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It spirals in toward the point of stable equilibrium taking an infinite amount of time to 
come to rest. 


Now turn around the sign of k. Take 

1 2 

U(x) --kx , k > 0 . 

2 


(4.30) 


The potential is a downward opening parabola: 



The lines represent trajectories of different energies. The energy equation is given by (4.29) 
with the sign of k reversed: 


E 



(4.31) 


This is the equation of a hyperbola. The phase diagram is 
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X(t) 



Let’s try to follow the trajectory of the particle shown with x 0 > 0 , x(0)<0 and E = 0. 

We have 

(4.32) 



which is the equation of a straight line. By integrating, we have 



(4.33) 


This means it takes an infinite amount of time to reach the origin. Notice, however, that 
x (t) = 0 is still a (trivial) solution the equation of motion, 


k 

x - — x = 0 . 
m 


(4.34) 


This leads to the picture for E = 0 trajectories: 
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X(t) 



I am suggesting that there is an isolated E = 0 solution at the origin, which can be reached 
only after an infinite amount of time in any continuous manner. I will symbolize this by 
drawing a circle around the limiting point of the motion. Thus, the equation of motion, 
(4.34), which is a linear second order equation, avoids the possibility of an intersection 
by making the time necessary to reach the intersection point infinite! Such points always 
correspond to unstable equilibrium positions. 


Nonlinear second order systems can have intersections in x , x phase space. Consider 
the potential (k > 0) 


U(x) = - — k | x | 3/2 . 
2 


(4.35) 


It corresponds to the (highly nonlinear) second order differential equation, 


3 x 
mx - — k —— = 0 • 


x 


|l/2 


(4.36) 


For E = 0 trajectories, 


Ik | | 3/4 

X = ± - |x| 


(4.37) 


The situation looks like: 
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The time necessary to reach the origin on path (1) or (2) is now finite. Integrating (4.37) 
in the two cases gives 


path (1): x 


x, 


path (2): x = - 


/ 


1/4 


_ 11 / 4 


X n 


if 

_ Ie t 

\m 4 


Setting x = 0 gives the time necessary to reach the origin, 


path (1): 


T n = 4 x 1 / 4 I 7) 


m 


1/2 


1— ll/4 1 

( m\ 

| x o| 1 

s 


1/2 


path (2): X 0 = 

Notice that at the origin on either path, we have 


(4.38) 

(4.39) 


(4.40) 

(4.41) 


path (1): 

x Lt 0 = 0 , 

(4.42) 

path (2): 

* u = 0 • 

(4.43) 


Thus the two paths do touch at the origin, 0, which is now reached in a finite amount of 
time. Note that this does not say such equations are not deterministic. It simply means, 
mathematically speaking, that the solutions, given initial x, x data are not unique. In fact, 
one can show that path (3) is the continuation of path (1) and path (4) is the continuation 
of path (2). 


Let’s now examine the phase diagram for the plane pendulum. Looking back at the potential, 
u(4>) , we find the qualitative picture: 
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Notice the periodicity of the phase diagram. Notice also that when the full range of is allowed, this 
phase diagram has both the bound motion seen in the harmonic oscillator phase diagram ( E = E 2 ) 
as well as the unbound motion of the potential with the downward opening parabola ( 
E = Ex)- The trajectory which separates these two types of motion has E = E 0 and is 
called the separatrix. 


Let’s again examine the trajectories for this special energy to see whether there is an intersection 
of trajectories at (f> = ±mt, n = 0, 1, 2, 3 ... Setting E 0 = 2mg in (4.19), we have 


i W 2 § 2 + 2mg^ sin 2 


2mg^ , 


(4.44) 


=> (j) = ± 2 




(4.45) 


Setting ([> = n + 6 (6 « 1) , we have 


cos 


(n + 6 ^ _ 5 

\ 2 ) 2 


r 


and the equation of motion in the close vicinity to <|> = % is given by 


5 



(4.46) 


This is the same type of equation of motion for the unstable system with U(x) --kx 2 , 

Eq.(4.32). Thus, although the plane pendulum is nonlinear, in the vicinity of the point 
4> = rut, the system behaves increasingly linearly. Thus, the picture near 4> = JC is essentially 
the same as before, and the system with energy E 0 takes an infinite amount of time to 
reach this point. 
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The Logistic Difference Equation 

The purpose of this Chapter is to introduce you to aspects of nonlinear systems. This is a 
vast subject and so we will only get a flavor of some of the new phenomenon and concepts 
which emerge. The best way to introduce this material is through a study of some properties 
of simple one-dimensional difference equations. These mappings are analogous to dynamical 
systems in ways we will discuss. 

A one-dimensional difference equation mapping is given by 

x n + ! = f ( x n ) (4.47) 

Let us consider the example (the “logistic map”) 

x n + ! = 1 - pxj; (4.48) 

for m in the open interval (0,2) and xn in the closed interval [-1,1]. A fixed point, x , of 
such a mapping is a value of x n such that x n + 1 = x n : 

X = 1 - [1 X , 
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(4.49) 


However, a very interesting thing happens when p is equal to or slightly larger than .75. We 
have a bifurcation in the fixed point diagram. The map after the first bifurcation looks like: 



We encounter other bifurcations for larger values of p as indicated: 


Pi - .75 

first 

p 2 - 1.2 5 

2 nd 

P3 = 1.3680989... 

3 rd 

p 4 = 1.3940461... 

f+ 


This process of bifurcation continues until at p 


bifurcation 

JJ 

ii 

u 

- 1.401155... an infinite number of 


bifurcations occurs. This result is very surprising because we are seeing an infinite amount 
of structure from a simple quadratic mapping. 
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* 2 

The bifurcation seen at p =.75 (where x* = —) is easily understood from the function 

3 

we are evaluating, a graph of which is below. 



- 1.000 x 1.000 


Notice that there is a perpendicular crossing of the function and the 4 5° line for this value of P. It 
is in this situation that a bifurcation occurs. Some iterative lines, which use x n as an input and show 
x n + 1 as an output, are shown in the above figure. Notice how they tend to spiral inward toward 
x , the fixed point. We will see in the context of a problem below that such a crossing 
becomes unstable if the magnitude of the slope of the line becomes greater than 1 at the 
fixed point. This is in fact what occurs for p slightly in excess of .75. Then the crossing 
point between the function and the 4 5° line becomes unstable and a stable square of iterative 
lines develops around this unstable point; see the below which is for p =.80, for which 
x* = 0.345 and 0.905. 



- 1.000 x 1 .000 
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There are certain universal features of quadratic mappings, one of which is the fact that 
(discovered by M. Feigenbaum) 

A 

Fn - , (4.50) 


for n»l where 6 - 4.669201609... (called the Feigenbaum constant). One can show 
that this implies that 


lim 

n -> oo 


On Fn-1 

O' n+ 1 - Fn 


6 . 


(4.51) 


which says that the ratio of bifurcation intervals tends to a fixed limit, although the ratios 
themselves are shrinking to zero. The reason these things were not discovered previously is 
because computers were not available to “explore” this relation, and because nobody thought 
such interesting structure would emerge from such a simple mapping. 


I do not mean to give you the impression that all is chaos for p > p<x> Within the range 
[ Pa, ,2] are also regions of order. For example, for m slightly in excess of 1.75 is a 
region of period 3. The structure of the fixed point diagram is both infinitely complex and 
astoundingly simple. The fixed point map is illustrated below. 


■ .000 


■I .000 


G.LOQ 
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Fractals 

An elementary example of a fractal is given by the “Cantor set of the middle third”, in 
which the middle third of a series of line segments is removed, as follows: 


step segments size 
OIL 

1 2 L/3 

2 4 L/9 

3 8 L/27 


This is an example of a fractal set of zero length. 

As another example of a fractal, consider the complex generalization of the quadratic mapping 
considered above, the so-called “Mandelbrot set”. It is generated by simply iterating 

z n + i = + c , |c| < 2, (4.52) 

in the complex plane. For each value c, one chooses z 0 = 0 and iterates the result. We 
then implement the algorithm: 

iteration test 

Nk + c | > 2 

k 

|zj( + c| < 2 

If, after N iterations we still have |z£ + c| < 
other color denoting how many iterations have been done. The surface of this infinitely 
structured object is the Mandelbrot set; it is a connected set. It is an example of a fractal 
which has infinite length! In the below picture the real axis is horizontal and the imaginary 
axis is vertical. 


result 

stop 

continue 

2 , one colors the point black; if not, some 
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Since the Mandelbrot set is generated from a difference equation which is just a generalization 
to complex numbers of the logistic equation, we might expect that they would be closely 
related. They are related in the following way. As one cuts through the Mandelbrot diagram 
along the real axis (going from right to left in the above diagram), one is looking at stable 
numbers on the real axis. The big bulb above corresponds to the stable single period line of 
fixed points in the bifurcation diagram. The first bifurcation corresponds to the appearance 
of the smaller blub to the left, the next bifurcation the the next smaller bulb and so forth. 
The point of infinite number of bifurcations corresponds to the point in the above axis 
where an infinite number of bulbs are stacked. Thus going along the real axis from right to 
left in the Mandelbrot set there is a one to one correspondence with the bifurcation cascade 
seen in the logistic equation fixed point map going from left to right. The Mandelbrot set 
has been called the most beautiful and complicated object in mathematics — one gets a 
glimpse of the infinite on a sheet of paper! 

Chaos in Physical Systems 

Why talk about fractals and maps in a course on particle dynamics? Because the phase 
space of driven, dissipative systems can be thought of as an iterative mapping, like our 
difference equation, and the trajectories are often fractal objects. Let us consider a system 
of first order equations, 

dx-, _ dx 2 _ dx, 

— = v i( x )'“7- = v 2 ( x )'— = v 3 (x), 

dt dt dt 

or 

dx _ 

— = v(x). (4.53) 
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Because these equations are the equations for particle motion, usually the x 3 , x 2 , etc. 
turn out to be generalized velocities or coordinates. The first such system which displayed 
chaotic motion was studied by Edward Lorenz, a meterologist, with an early computer 
called a “Royal McBee.” His system of equations were 

^ = 6(y - x) , 
dt 
dy 

— = rx - y - xz , 
dt 

dz K 
— = xy - bz. 

dt 

Here 6, r and b are constants and obviously Xi=x x 2 =y and x 3 =z • The trajectories 
in the three dimensional phase space of this system result in chaotic motion and possess a 
fractal structure. 
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The damped driven oscillator may be characterized as in (4.53). One may write 


dx 


dt 1 
dy 
dt ' 
dz 
dt ' 


= -cy - sin x + F cos z, (c > 0) 
= 0 ). (a constant) 


(4.54) 


The physical interpretation of each term is 

z: driving phase 

x: angular position of pendulum 

y: dimensionless angular velocity 


and the above represents: 



The above first order equations are in convenient, dimensionless form, with c, F, t ' 
and co related to m, i , b (damping constant), N d (external torque constant) and 0) d 
(driving frequency) by 


C = 


mgf 


C0O = 


l 


F = 


N, 


mqf 


, t' = 


0 ), 


co 0 t, 0) = — 
0) n 


The three equations, (4.54), are equivalent to the single second order differential equation, 
mf 2 x = -bx - mgf sin x + N d cos (o) d t + cp) , (4.55) 

which is just Newton’s law for the damped, driven pendulum. 
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One necessary condition for chaotic behavior is that the equations in the above first order 
form contain a nonlinear term coupling one variable to another. This is seen in the -sin 
x term in the second equation. Surprisingly, a system of three such dynamical equations 
is already sufficient to display chaotic behavior. 

The true phase space of the damped, driven pendulum is three-dimensional because of the 
three first order equations above. Physically, this is because it is no longer enough to give 
x ( t ) and x(t) to specify the initial state of the system. It is necessary also to specify 
the phase of the driving term, which is given by (co d t+ (|)) in (4.55). However, all the 
trajectories in this three-dimensional space may be projected back on the old x(t), x(t) 
plane. This projected version of phase space can display many types of motion for a given 
set of parameters C, F and 0 ) and a given set of initial conditions. Analytic methods are 
of limited usefulness here and computer simulations must be used in order to gain insight 
into the types of motions that occur. 

There are motions where the pendulum moves in a completely periodic manner, although 
the period of motion may not be the same as the period of the driving force. These are 
called closed trajectories, an example of which is the harmonic oscillator phase space above. 
There are also motions which display no periodicity at all, although there may be submotions 
which continue to reoccur. This is chaotic motion. After initial transients die out the damped 
pendulum evolves to a steady state trajectory called an attractor. The simplest example of 
an attractor is the point of stable equilibrium in the two-dimensional phase space of the 
damped, free pendulum. An attactor can also be a closed trajectory. 

For chaotic motion in the driven, damped pendulum, the projected x(t), x(t) phase 
space diagram is usually shown as a filled-in portion of phase space. Not all velocities, x(t) 

, are available to such a system because the system is dissipative (has a positive damping 
constant.) Of course, it can never reach every point of projected phase space in a finite 
amount of time. Actually, what one sees after a simulation is done is a tight, packed set 
of intersecting and non-intersecting trajectories. (These lines do not intersect in the full 
three-dimensional phase space.) The infinitely structured geometrical object that the system 
tends toward in the steady state when chaos is involved is called a strange attractor, and 
is another example of a fractal. 

Because of the intricacy of the projected phase space plots and because the full three- 

dimensional phase space is hard to show, one needs tools to help uncover structures present 

in the phase space. One such tool is called a Poincaire section. Such a plot shows an overlap 

of a large number of “snapshots” of the x(t), x ( t ) state of the system. If we agree to take 

2 jc 

these snapshots at a period T d = — , where w is the driving angular frequency, we can 

0) 

see at once whether the system is periodic or chaotic. 
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The simplest thing one can observe is a finite number of points. This immediately implies 
that the motion is a closed, periodic trajectory. Of course, the location of the points in the 
Poincaire section is a function of the phase, <t>, of the projection, so this does not have any 
particular meaning. If X points appear, the resulting motion could have been completed 

in a period, X, given by nx d , — X d , — x d , • • • , where — , — ,. . . itself is not an 

2 3 2 3 

integer. If in fact the closed trajectory has X = 2 x d , this is called “period 2 motion”. The 


transition from period one to period two motion as the parameters of the pendulum are 
changed is called period doubling and is the analog of the bifurcations seen in the fixed-point 
difference equation mapping. A Poincaire section can also be regarded as a mapping, using 
particle dynamics instead of a quadratic formula, from a state n to n + 1 The trajectories 
that the system settles down to in the steady state for the pendulum are the analog of the 
fixed points, x , in the difference equation. Such mappings, we have learned, can involve 
chaotic results and fractals in either context. 


Since chaotic dynamical systems are still deterministic, true bifurcations in phase space do 
not exist. However, this bit of wisdom is essentially useless for true chaotic motion since 
neighboring trajectories, no matter how close to begin with, eventually diverge exponentially. 
(More on how to characterize this divergence in a moment.) However as I said above, 
bifurcations do exist in their parameter space descriptions. For example, if one measures 
x (angular velocity) as a function of the driving torque, F, on a Poincaire map, one sees 
bifucations appear as F is increased. These bifucations will be present either because a change 
in the initial conditions produces a new steady state or because a true period doubling 
trajectory has evolved. Eventually the period doublings merge into chaos, just as in the 
one-dimensional mappings. This period doubling route to chaos is again characterized by 
6, the Feigenbaum number and has been seen in physical systems. 

Dissipative Phase Space 

How can chaos in dynamical systems be identified? Let us first try to understand the 
evolution of dissipative systems in phase space. Consider a simple case of two-dimensional 
particle trajectories in phase space as shown below: 



v(x lf x 2 ) 


volume," V(0) 
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dx 

Each particle has been assigned a “velocity”, v(x lf x 2 ) = -, which is a function of 

dt 

position. The general formula for the instantaneous flux out of the (static) volume V ( 0 ) is 
f = $ ds p(x lf x 2 ) v • n , (4.56) 

S 

where the closed “surface”, S, is simply a rectangular perimeter. We have introduced in 
(4.56) a position dependent density of points, p(x lf x 2 ). Setting the density of points to 
one per unit volume, p(x x , x 2 ) = l and using Gauss’ law in two dimensions we obtain, 


f = J* dV V • v . 


(4.57) 


The result (4.57) applies to an arbitrary number of dimensions if we generalize the meanings 
of the surface " S" and the volume “V" i n (4.56) and (4.57). We will allow a dynamic volume, 
V ( t ), to expand or contact to contain the same number of particles. The instantaneous 
change in this volume is just given by 


8V(t) = 8t f(t) , 


dV 

=> - = fdV V • v. 

dt v 

For example, for the system (4.54), and using dimensionless time, t' 


(4.58) 

(4.59) 


v 


= (x 2 , -cx 2 - sin x : + F cos x 3 , to) , 


V • v = -c , 


dV 


dt 


; = -cV(t ), 


=> V(t' ) = V(0)e“ ct . (4.60) 

From this we learn that the phase space is exponentially contracting ( c > 0 ) This is a 
general characteristic of dissipative systems. V • v = 0 characterizes a nondissipative 
system whose phase space volume can stretch and change, but whose volume is fixed. One 
can show that systems which possess a Hamiltonian (Ch.7) behave in this manner; this is 
called Liouville’s theorem. 
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Lyapunov Exponents 

Even though the phase space volume is contracting for these systems, chaos emerges as a 
stretching of the phase space volume so that nearby phase space points eventually exponentially 
diverge. We can characterize the divergence of trajectories in a phase space plots by the 
Lyapunov exponents. They are related to the previous phase space considerations via the 
expression (£i is initial separation in the i— direction between x.^ and x.^ at t=0) 


lim |x ± (t) - x ± (t)| = s i e >_it (4.61) 

Since the overall volume is given by (imagine a rectangle) 
lim V(t) - Yl | x 'i M " x i Ml 

i 

<n is the product symbol over i-1,2,3) we find 
i 

2 m 

lim V(t) - Yi e i e > ' it = ([\ e ± ) e 1 , 

i i 

^(oT 

and the total volume grows with an exponent which is a sum of the individual Lyapunov 
exponents times the time. Thus, for the damped oscillator, 

3 

2 X i = “ c ' (4.64) 

i = 1 

and the sum of the exponents is just the negative of the damping constant. Although the 
sum of exponents is assured negative, only one exponent has to be positive to indicate an 
exponential stretching of the phase space volume, V (t). Thus, all of the exponents must be 
examined before a possible chaotic motion can be identified. (The calculation of the individual 
Lyapunov exponents for this system is much more difficult than getting the sum as above.) 

The Lyapunov exponent for one-dimensional mappings can be numerically calculated. Lor 
small enough £ one assumes 

f(x 0 + e) - f(x 0 ) - £ ^ | Xo , (4.65) 

where x n+1 = f (x n ) Next, one has 

f (2> (x 0 + e) - f <2> (x 0 ) - £ | • (4.66) 

dx 

where f* 2 \x 0 ) is called the second iterate of f (x 0 ), 


(4.62) 


(4.63) 
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f<2> ( x o) s f ( f ( x o)) * 


(4.67) 


Let us call x-l = f(x 0 ), x 2 = f(x x ) = f(f(x 0 ) , etc. Taking the derivative of f* 2 \x 0 ) 
now gives by the chain rule, 

df (2) , df (2) 


dx 

After n iterations, 


_ df | _ df | df . 

dx lf(Xo) dx Xo_ dx Xl dx x ° 


(4.68) 


f (n> (x + e) - f <n) (x) = E 


df , df 
dx Xn ' 1 dx 


df 

dx 


(4.69) 


The definition of the Lyapunov exponent is essentially the same as before, with the iterate 
number, n, playing the role of time (x 0 = x 0 + e), 


lim 

n-*oo 


X n “ X n 

n n 


A, = lim — In 
n 


£(n) 

( x o + £ ) - f<n) 

( X o) 


E 


' df 

1 df 1 


dx 

lxn - 1 dx 1x0 

i 


= e 


nA. 


(4.70) 
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= lim 

n-»oo 


1 

n 


n-1 


5 > 



(4.71) 


This can be numerically calculated as a function of the parameters of the mapping and 
chaotic regions can be identified. 

If one of the Lyapunov exponents is positive, this is a necessary but not sufficient condition 
for dynamical chaos to emerge. For example, if we (somewhat unrealistically) changed the 
damping constant to be negative for the harmonic oscillator, we would have a system which 
is exponentially sensitive to initial conditions, but which is linear and so is not chaotic. On 
the other hand, even though the sum of the Lyapunov exponents for the damped pendulum 
is negative, only one has to be positive for this nonlinear system to exhibit chaotic motion. 
It is a combination of the nonlinear behavior and the exponential growth in at least one 
direction which are characteristic of dynamical chaos. 

The Intermittent Transition to Chaos 

We have learned about one path or way of transitioning from simple, predictable periodic 
motion to chaotic motion. It was called period doubling and is illustrated in the logistic 
fixed point map above. However, there are other ways systems can transition from simple 
periodic motions to chaotic ones. A system may appear to be periodic or well behaved for 
a long time and then diplay “bursts” of nonperiodic behavior. As one changes the system 
parameters closer to the chaotic regime, these bursts take place more frequently until fully 
chaotic behavior results. 

This intermittent transition to chaos can also be understood from our one dimensional 
logistic map. In the above fixed point map near the region of the large period 3 window 
for m slightly above 1.75, the system is in a three period, that is, after a few iterations it 
settles down is a stable system with three fixed points. 



too 


1.682 P 1.897 
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However, notice the wall of chaos to the left of this periodic region. When we consider 
values just slightly inside this chaotic region, we almost have a three period. That is, the 
system behaves for a long time as if it is in a three period, but with occasional bursts of 
chaos. The further we imbed ourselves in this chaotic region, the less we see of the three 
period and the more chaos we see. How can this sort of behavior be understood? 

Amazingly, this behavior may be understood very simply. The function f(x n ) = 1 - px 3 
with p= 1.7 5 looks like the following. 



-1.000 X 1.000 


However, the third iterate of this function, f* 3 \x n ) = f(f(f(x n ) )), looks like the following: 



- 1.000 x 1.000 
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Notice the four values in x where f* 3 \x n ) touches the 45° line. One of these is unstable 
(which one?); the other three are fixed points of the mapping, x , and correspond to the 
three lines of bifurcation in the nonchaotic window. Let us zoom in on the middle these 
fixed point regions near x=0 when p has a value of 1.7495, slightly into the the chaotic 
region. Notice that a “tangency channel” has developed where the function is close to, but 
not actually touching the 45° line. 



-0.014 X 0.058 


Also shown in this figure are some iterates of f* 3 \x n ) going into the tangency channel. 
Notice the the movement of the iterates through the channel becomes quite slow in the 
bottleneck part of the channel. This corresponds to the situation of having a nearly periodic 
three cycle system. However, eventually the iterates tunnel through the channel, reaching 
the chaotic region outside, corresponding to the “bursts” of chaotic motion. Put another 
way, the closer we are to true tangency (at p= 1.7 5), the narrower the channel and the 
more time the system displays regular, predictable motion. The bursts occur more often the 
further we go into in the chaotic region because the tangency channel is wider. 

The number of iterations that the system stays in the channel, displaying quasi-periodic 
motion, can be understood from the above difference equation. In the vicinity of the 
tangency, the third iterate, , may be expanded in a Taylor series to read 

f (3) (x n ) = x n + a c (x n - x*) 2 + b c (p c - p) , (4.72) 

where p c = 1.75 and x* is one of the three fixed point values. After the transformation 
y n = (x n - x )/b c , the recusion relation x n + 3 = f^ 3 \x n ) in the vicinity of the 
channel may be written 
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y n+ 3 - y n = + e, 


(4.73) 


where a = a c b c and E=p c -p>0. When the number of steps in the channel is 
large, (4.73) may be approximated by a differential equation, 

— = ay 2 + e . (4.74) 

dn 


This differential equation may be solved by separation of variables for n: 


n (y 0 ut. y in ) 



(4.75) 


where y in and y out are the entrance and exit values of y in the channel. Notice that for 

y Q ut' yin fixed ’ n(yout'y±n) g° es to infinit y like -p as e -> o. 

v£ 

Tliis describes the transition to chaos that one is seeing in the logistic map near p=l. 7 5 
and many other regions. The transition to chaos in a period doubling manner and in an 
intermittent manner are very closely related. We have seen that period doubling occurs when 
there is a perpendicular crossing of the 4 5° line. Intermittent chaos emerges when instead 
the crossing is a tangent to the line. In another sense they are also related. Notice in the 
above graph of the 3 period we have been discussing, that the other side of the periodic 
window displays a period doubling cascade to chaos. Such windows always displays these 
two, complimentary, transition routes. 


CHAPTER 4 PROBLEMS 

1. Verify that the change of variables 



converts Eq.(4.8) into (4.9). 

2. By making the substitution x = sin0 in (4.27), show that the elliptic integral 
F ( 1, k) may be approximated by, 

Jt / 1 o 9 4 \ 

F(1,k) = — 1 + - k 2 + — k 4 +...), 

2 \ 4 64 I 

for k«l. 

3. A free (not driven) underdamped oscillator is released or set in motion under the 
following initial conditions: 

a) x ( 0 ) = x o , x(0) =0 ( x 0 > 0 ) 

b) x(0) = 0 , x(0) = x 0 ( x 0 >0) . 
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Sketch phase diagrams of the resulting motion in the two cases. Make sure your phase 
trajectories have arrows. 

4. a) The potential function (A,X>0) 

U(x) = A(x 4 -X 2 x 2 ) , 

looks like 



Sketch a phase space diagram showing qualitative trajectories for energies Ei f 2 , 3 , 4 * 

(E 3 = 0 . ) (These trajectories must have a direction.) Identify points of stable and unstable 
equilbrium in your sketch. 

5. A one-dimensional potential has the form, U(x) = Uoe -ax2 



0 x 0 x 


a) Sketch a qualitative phase space diagram corresponding to system total 
energies Ei, E 2 and E 3 (E 2 = Uo» ) 

b) Starting from x=Xo with total energy E = Uo, find the time necessary to 
reach the origin, x=0. Is it finite or infinite? 
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6. A one-dimensional potential energy, U(x) is given by U(x) = -K|x| n , where 
| x | is the absolute value of x, K is a positive constant and n > 1 For what 

values of "n" does a particle, starting at Xq^ 0 take a finite time to reach the 

1 

top of the potential hill when the total energy, E — 2 mx 2 is zero? 

7. The steady-state response of a damped oscillator to a periodic force, 

mx + bx + kx = F 0 cos (art). 


is 


x(t) = D(a))cos(a)t - 6), 
h „ k 

where ( 2 (3 = 

D(co) 


, O) 0 , Aq 

m 

A„ 


0 r “U m 

m m 111 


yf( 0 )q - a) 2 ) 2 + (2a)(3) 2 


, tan 5 = 


2|3a) 


0)o - 0) 2 


a) Show that the trajectories in (projected) x , x phase space are ellipses. 

b) Show that the driven oscillation total energy is constant in time if 0 )=(Oq. 
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c) We know that the (unprojected) phase space of the same oscillator when it 

is free (i.e., undamped, nondriven) is also an ellipse. Show that these two 

2 2 

ellip ses coincide only if w = co 0 and their total energies are equal. 

8. Given the one-dimensional mapping (all x n > 0), 

x n+l — "\/ x n i 

a) Find all of the fixed points (there are two) and show that one is stable, the 
other is not. 

b) Calculate the unique Lyapunov exponent for this mapping. 

9. We investigated the mapping (p E ( 0,2 ) , x E [-1,1]) 

x n+ i = 1 - uxg 

for fixed points. We found period doubling at p=.75,1.25,.... Show that this 
equation is equivalent to the mapping 

y n+ i = ay n ( 1 - y n ) • 

for a E ( 0,4 ) and y E [ 0,1 ]. At what value of a does the first period doubling 
occur? 

10. When n becomes large, the period doubling values, fin/ for the mapping in 
prob.9 obey 

A_ 

fin ’** fipQ - gn r 

where U ffl , 6 and A are constants. From this show that 
lim ~ M ' n ~ 1 = 6. 

fin + l - fin 

This says that the ratio of intervals tends to a fixed limit (although the intervals themselves 
are shrinking to zero). 

11. a) A real linear mapping is given by 

x n +i = Px n + c. 
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where | (3 | < 1 and c is a real constant. Show that the single fixed point of the mapping 

Q 

is given by x = --— • What happens for | (3 | > 1 ? 

b) Explain the relevance of the result in (a) for a linear mapping to the idea of 
bifurcations seen in nonlinear mappings. 


12. As shown in problem 11, the real linear finite difference equation, 


x n +i = |3x n + c. 


0 

with | (3 | < 1 and c a constant has a single (stable) fixed point at x* = -— . Let’s 

say we decide to find x by brute force iteration. Picking an initial value Xq , we iterate 
the difference equation until we obtain | x N -x* | < £, where e is a small positive number. 
Show that the number of iterations needed to achive this accuracy is given by 


N = 


In 


x o ~ x 


In I b 


showing convergence is worst for | (3 | — 1. 


13. Consider the non-linear difference equation given by 


x n+ i = F(x n ) = a x n ( l-x n ) 2 . 


for 0<x n <l and ct >1. It looks like: 



o.ooo x 1.000 

x 
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a) Find the location of the fixed points of the mapping for 0 < x n < 1 • 

b) At what value of a does the first “period doubling” bifurcation occur? 

14. Consider the “tent map” given by (0<cc<l) 

x n+ i = 2ax n , 0<x n <l/2, 
x n+ i = 2a(l-x n ), l/2<x n <l. 

a) Find the values of the fixed points of this map. (There are two.) 

b) Examine each point in (a) for stabilty. For what values of a are the fixed 
points stable? [Hint: Think about slope.] 

15. a) Show that the transformation, 

2 i— 

Yn = ~ sin -1 (y/x n ) , 

in x n+ i = a x n ( l-x n ) for a-4 converts it into a tent map, 

Yn+i = 2y n , 0<y n <l/2, 

Yn+i = 2(l-y n ) , 0<y n <l/2 . 

b) Therefore calculate the Lyupanov exponent of the above transformation for 
a=4. (Ans: Cn (2) .) 

16. Consider the three coupled first-order differential equations that Edward Lorenz 
considered (same as in this text), 

dXi 

dt" = 8(X2- Xl) ' 
dx 2 

dF = rxi - x 2 - xix 3 , 
dx 3 

dF = x i x 2 - bx 3 , 

where 6, r and b are constants. The rate at which a volume of the configuration 
space changes is given by 

V(t) = V(0)e Xt . 
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Find A,. Under what conditions is this system dissipative (does the configuration 
space shrink in time)? Given the values 6=10, r=-2, b=-10, is the system 
dissipative? Can you tell if the system is chaotic? Answer the same two questions for 

6=10, r=-2, b=-12. 

17.Show that when the three equations in (4.54) are combined, they give the single 
equation, (4.55). 

18. Evaluate the Lyapunov exponent of the linear mapping, 
x n +i = (3x n + c. 


Show that it is negative only for | P | < 1 • 

19. By numerical calculation, find the values of the three stable fixed points of the 
mapping f (3) (x n ) = f(f(f(x n ))), f(x n ) = 1 - \ix 2 n , for E=1.76. 

20.Show that when Eq.(4.72) and the transformation y n = (x n - x )/ b c are 
used in the recusion relation x n + 3 = f* 3 *(x n ), the recursion relation in the 
vicinity of a tangency channel may be written 

Yn + 3 - Yn = a Y 2 n + E- 


21. (a) Use the relation Eq.(4.71) and the result (4.73) to show that the Lyapunov 
exponent for a single pass through the tangency channel is approximately 
given by 


y out 


^(you^ yin) s 


n (y out ' y in ) Y: 


f d y 


ln(l + 2ay) 
ay 2 + e 


[Flint: Make the same sort of approximation in going from a sum to an integral that 
we did in the text in solving (4.74).] 


b) Make the approximation that 2ay is small inside the integral above and do 
the integral. From this predict the behavior of A,(y out , y in ) on E as £ —* 0 
( y±n ' yout ' held constant). 
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5 GRAVITATION 


NEWTON'S LAW OF GRAVITATION 

Newton’s law of gravitation for point masses is (force on m due to M): 
mM ~ 

F = -G — e r . (5.1) 

r 

Picture: (M usually more massive than m) 




> Apply now 


REDEFINE YOUR FUTURE 

AXA GLOBAL GRADUATE 
PROGRAM 2015 


redefining /standards » 
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Obviously an idealization. (Masses are not points. Also, forces are not transmitted 
instantaneously .) More general form for a mass distribution: 



where 


r" 


= r - 


r' 


(can not be taken outside the integral in general) 


Even this is an idealization for, say, two planets: 



dM -»• pdr'ldv' , dm p 2 (r)dv , 

F - -G/Jdv'dv p . (Oft (*) ^ . 


(5.4) 
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Gravitational Potential 


Things are beginning to get complicated. We need some simplifications. First simplification: 
the idea of a gravitational potential. Go back to point masses case. Notice 

v r G —) = _GMm v t)' 


t 


\ 


3x, 


V V? + X 2 + X 3 / 


2x, 


2 (x; + X* + x;) 3 ' 2 


Likewise for x 2 , x 3 • 


- / GMm\ GMm „ 

V-1 = —r- e = -F . 




Mm 


Therefore, with U(r) = -G - (makes U=0 at r=oo), we have 

r 

F = -VU 


(5.5) 


(5.6) 


u 


Define = — and call it the “gravitational potential”. It represents the work/mass to 


m 


move a test mass from one point to another. In the case of M being an extended mass, it 
is easy to see that 


r P ( r ') 

U(r) = -Gm I — L dv ' , 

r p (r ') 

=> <F(r) = -G f ^ dv ' . 
Verification: 

-VU = Gm f p (r ') V — dv ' . 


(5.7) 

(5.8) 


We have 


dx, 


(x; - x,) 2 + (x 2 - x 2 ) 2 + ( X ; - x 3 ) 2 


( x i - x i) 


f I" '_ F ^ 

-VU = Gm f dv ' p (r ') a- —L 


r* "2 
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The point here is that is much simpler to calculate than F. One usually defines 

g = — => g = -V<f> . ( 5 . 9 ) 

m 

Near the Earth’s surface: 

|g| - 980 cm/sec 2 direction: “downward” 

Modifications for Extended Objects 

We can recover the usual form of gravitational potential energy near the Earth’s surface, with 
the help of a result we will prove momentarily The result we need is that for a spherical 
mass distribution 


, „ r P ( r ') , , GM 

O = -G EA— L dv =-, 

J r " R 


U = - 


GMm\ 

R ) 


(5.10) 


where R is the distance to the center of the sphere. Pictorially, 



M 


<fi(R) 


R 


Then we have (R — a + h , h « a) 
1 1 


1_ 

R 


a + h 


a(l + * 




GMm ( h\ GMm 

U --1-=-+ mgh 

a \ a a 


GM / , 2 \ 

g = —- (= 980 cm / sec ) 


5.11) 


f (assuming a perfectly spherical Earth) 

Let’s now prove the result that the potential outside a spherically symmetric distribution 
of matter is independent of the size of the distribution and acts as if all the mass were 
concentrated at the sphere’s center. (This is the problem which spurred Newton to the 
creation of calculus. Note that r-> R here.) 
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density: p( 



0 : origin of coordinates 


(r) = -Gj dv ' P ^ 1 , (dv ' = r' 2 sin 0 d0 dcp dr' ) 
=> O(R) = -2nG f dr ' p (r ') r' 2 f d0 . 


( 5 . 12 ) 

(5.13) 


Now 


R = r + r , 

=> r " 2 = r' 2 + R 2 - 2Rr' cos 0, 

/ \ r p (r ') r 12 ~ 

=> 0>(r) = -2kGJ dr ' 1-2—-J 


2r ' R sin 0d0 


2r ' R *£ ,2 + R 2 - 2Rr ' cos 0 

In the later form, the 0 integral is a perfect differential. We have 


<D(R) = -2jtG J- 2^r' 2 + R 2 - 2Rr' cos 0 

o 2R 

Notice that the above now involves absolute values 
2ttC a 

$(r) -- f dr' r' p(r')(|r' +R| - |r - r'|). 

R o 

If the point P is outside the mass distribution, we have |R - r'| = R 
above becomes 

<b(R) = - — fdr'r' 2 p(r'). 

R J o 

However, the total mass of the distribution is 

a 

M = 4jtJ dr' r' 2 p(r') , 


(5.14) 

(5.15) 

(5.16) 


(5.17) 


(5.18) 
r 1 , and the 

(5.19) 


( 5 . 20 ) 
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so therefore 


<£(r) = - 


GM 

R 


and thus 


GMm - 



(5.21) 


( 5 . 22 ) 


which says that the distributed mass just behaves as if it were concentrated at a point at the 
sphere’s center. Interestingly, the same essential form holds even if the point P is located 
within the mass distribution, r'< R. Then, one may show (HW problem) that the same 
form holds, 

GM(R)m - 

F --TT" e R , (5.23) 


where 

R 

m(r) ■ 4tt J dr' r' 2 p(r') , 

0 


(5.24) 


is the total mass within a sphere of radius R. See the following picture of the situation. 



Eotvos Experiment on Composition Dependence of Gravitational Forces 

I would now like to pick up the thread of the argument concerning the supposed equality 
between gravitational and inertial mass (Chapter 2). This equality can be tested. Let me 
describe the Edtvds experiment. (Follows J. Weber, “General Relativity and Gravity Waves”). 
This experiment was done in 1905 and established the equivalence of these two quantities 
for specific materials, to about one part in 100 million. 

Here is the basic idea of a torsion balance: 
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(Also fictional 
"forces ", ) 


(Gravitational mass ) 1|2 = M 1;2 
(Inertial mass ) lj2 = m lj2 


We will define the ratio of gravitational to inertial mass to be 


a 


1,2 


M 

m 


1,2 

1,2 


Origin of “fictional” forces is the spinning Earth: 
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|force| 


Apply to the Earth: 


mv 2 

- = moo a 

a 



Net torque on torsion balance: 

T - | x (G. + i.) - | x (G 2 + I 2 ), 


F = G 1 + G 2 + + I 2 , (along wire) 

F • T 


T = 


, (crucial thin+g to calculate) 


T = 


( G 1 + G 2 + I i + I 2 ) * ~ x ( G i + I i) 


bJ 

( e . + 

1 X 1 

( e , + *01 

Ii 

p| 



(5.25) 

(5.26) 


(5.27) 


The second term above just doubles the result from the first. We have the reduction: 

( G i + ij x ( g 2 + i 2 ) 

= (-g M x e r + m 1 ca 2 R cos 4>e m ) x ( _ 9 M 2 e r + ™ 2 a) 2 R cos <t>e m ) > 
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= -g oo 2 R cos 4> (MjIHj - x e m ) . (5.28) 

Putting the pieces together yields the component of torque causing a rotation, 

-g o) 2 R cos 4) m 1 m 2 (a 1 - a 2 )b • (e r x e m ) 

J-n — i-1 • 

F 

Since we are working only to first order in - a 2 ) ,, we may use 

|f| ~ g(M, + M 2 ), (5.30) 

- ID (5.31) 


m 1 m 2 


+ M 2 


yielding finally 


Tn = -pw 2 R cos 4» (ai-a 2 ) b • ( e r x e m ) , 


(5.32) 


where b * ( e r x e m ) is a max in E-W direction. 


Description of the actual measurement: 


The torsion balance is brought to equilibrium by turning it until it points in an E - W 
direction. One can not tell yet if there is a nonzero component of torque, T 11; since it is 
in equilibrium. Now turn the apparatus 180° reversing the sense of b. If ct 1 ^ a 2 there 
will be a torque which will rotate the rod relative to the frame which supports the balance 
by an angle twice the (unobserved) initial deflection. See the diagram below. 


position 


looking down 
on apparatus: 



Best to do the expt. at ^ ~ 4 5° since (with b oriented E - W) 
cos {f)|b • (e r x e m )| = b sin 4 > cos 4 > . 
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Later experiments: Dicke, Roll, Krotkov established equivalence between gravitational and 
inertial masses, again for specific materials, to about one part in 100,000 million (1 in 
10—11!!! See the October ’99 Physics Today article from Clifford Will for more updated 
information on gravitational composition dependence experiments. In the literature, the 
postulate that there is no composition dependence to gravitational forces is called the “weak 
equivalence principle.” 

CHAPTER 5 PROBLEMS 

1. Consider two point masses which are stationary in some inertial reference frame. 


I* 

-•- 

Mi 


R 


♦ 


M 


2 


a) Evaluate and sketch as a function of x (measured along the distance between 
the planets) the gravitational potential, <f> (x) , a test mass would experience. 

b) Find the value of x (other than x =±oo) for which the test mass feels zero 
gravitational force. Is this point stable? 


2 . Show, for a spherical mass distribution, that the gravitational force on a small 
test particle of mass m located inside a planet of radius " a" and a distance R 
from it’s center (a>R) that 


F 


GM(R)m 


R z 


r 


where 

R 

m(r) ■ 4jt J dr' r' 2 p(r') , 

0 

is the total mass within the sphere of radius R. See the following picture of the situation. 
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mass within sphere 
radius R = M(R) 


of 


3. Find the force of gravity on a small spacecraft of mass m, distance x" from the 
middle of a long, straight space station (Babylon 5) of uniform density per unit 
length, p, and length L. Take your origin of coordinates, 0 , in the middle of 
the space station, as shown. 



P m 


_ L 0 

2 

L x" 

2 


4. A small rock of mass m is orbiting the Babylon 5 space station. It orbits in a 
plane perpendicular to the station, in a circle of radius R away from the space 
station’s midpoint. The station has a length L and a constant mass/length, p. 


R 


L 

T 


o 


p 


L 

T 


Find the period of the orbit, T, in terms of R, L and p. [Hint: Get the gravitational 

2 

acceleration and set it equal to centripetal acceleration, — . I have a table of integrals.] 

R 

5. Scientific Superman digs a hole completely through planet Krypton (which has 
a spherically symmetric mass distribution) and measures the force of gravity as a 
function of the distance from Krypton’s center, R. He finds 

F= - H(R)e R , 
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where e R is the usual unit vector pointing away from the planet’s center. Show that this 
means Krypton’s mass density, p ( R ), is given by 


p (R) = p 55 (R2 H(R)) • 


[Hint: Remember that 

GM(R)m - 

F = - r 2 e R , 

where M ( R ) is the mass contained within a distance R from the planet’s center.] 

6. Scientific Superman finds that Krypton’s mass density, p, is a constant in R. He 
now drops an evil villain down the hole he dug through Krypton. Show that the 
motion of the evil villain through Krypton is simple harmonic. Find the villain’s 
period of motion in terms of the assumed constant p. 

7. Two tiny planets of equal mass, M, attract one another and move along the line 
connecting them (no rotational motion). They are initially stationary and located 
symmetrical distances Xq and -Xq from the origin, 0, as shown. 



Show that the time necessary for the planets to crash into one another is given by 
2 


T 2 = 


71 


2G M 


(x 0 f • 


[Hint: I had to look up an integral. This result is actually a special case of what is known 
as Kepler’s third law of motion which says that the square of the period of motion is 
proportional to the cube of an orbit parameter.] 
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6 CALCULUS OF VARIATIONS 


EULER-LAGRANGE EQUATION 

In ordinary calculus, we are used to finding stationary values (max’s, min’s or saddle points) 
of functions f ( x ). We now consider the problem of finding a function to make an integral 
stationary 

Simplest form: 

x 2 

J = f f [y(x) , y ' (X) ; x] dx . 

X 1 

dep. variable ft t indep. variable 

dy 

- an indep. dep. variable 

dx 


Some examples of f ( y, y ' ; x ) are: 


b) 


•/i 

+ 

' 2 

y 

M 

+ 

' 2 

y 


c) a (x) y + b (x) y 2 


Just an in differential calculus, the vanishing of a first derivative is a necessary, but not 
sufficient condition for a max or min, so in the calculus of variations one speaks of 1 st or 
2 nd variations of J. Here we will only work with the first variation and rely on geometrical 
or physical reasoning to decide whether we have found a max, min or inflection point. The 
point of learning this is the next chapter, where it will help us formulate a new, more general, 
view of classical dynamics. This chapter, however, is largely mathematics and geometry. 

We will often need the following theorem in justifying statements in this Chapter: 

If F(x) is a continuous real function for Xi ^ x £ x 2 and if 


x 2 

J dx r|(x) F(x) 

XI 


0 


( 6 . 2 ) 


for all real, continuous, once-differentiable functions T| (x) such that 
11 (xi) = ti(x 2 ) = 0 , 


(6.3) 
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then F(x) vanishes identically. 


I will not prove it, but here is an “heuristic” 5 argument. 


F(x) 


F(x) 


x 



In each region (say a ^ x ^ b above) we have 
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b 

f dx v|(x) F(x) > 0 . (6.4) 

a 

But we have that 

x 2 a b 

J* dx ri(x) F(x) = J dx r|(x) F(x) + J dx r|(x) F(x) + ... . (6.5) 

Xj_ x x a 

The only way to reconcile (6.2), (6.4) and (6.5) is if F (x) = 0 everywhere in the interval 
x x < x < x 2 . 

We can now go ahead with the rest of the development of this chapter. 

We will need to develop some formalism first. Consider the two paths shown, one of which 
is assumed to be the extremum path. 



y(a,x) is a neighboring or companion function. 

y(a,x) = y(x) + ap (x) ; y'(a,x) = y'(x) + ap'(x). (6.6) 


a is a continuous real parameter and p (x) is an arbitrary, once differentiable function for 
which p(x 2 ) = r|(x 2 ) = 0. That is, we are keeping the end points fixed. We can then 
form the comparison integral (assuming y (x) is twice differentiable) 

x 2 

J(a) - f f [y(a, x),y' (a, x) ; x] dx , (6.7) 

x i 

where J ( 0 ) is the extremum sought. The condition is then just 


dJ(a) | 

da l“ =0 


( 6 . 8 ) 


Using the chain rule (I am using simplified notation), 
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df 

da 

dJ | * 2 \ df dy df dy'l 

da ^ a = 0 J [dy to dy' 3a J ' 

where, using the previous definitions 


(6.9) 


dy 

— = iT(x) , 

da 

dy' dri 

— = — 1 = ri'(x). 

da dx 


r 


( 6 . 10 ) 


df df 

We will call objects like ~— r ~—: “functional derivatives” since we are taking derivatives 

dy dy 

with respect to (unknown) functions, not variables. 


We have now 

dJ | 


* 2 [df df 1 

f I — TKx) + — 11 (X) | dx . 

XI L ' 


da '“= u J [dy " ' dy* 1 ’ J 


( 6 . 11 ) 


Let us concentrate on the second term on the right hand side. By integration by parts, 


f V ( x ) dx = fi(x) | 2 - f j- f~~] r l(x)dx . 
^ dy dy *i { dx (dy / 


( 6 . 12 ) 


df 


However, since r| ( x 1 2 ) = 0 and assuming - is nonsingular at x 1 2 we have that the 

dy' 

boundary term vanishes. Therefore, the stationary condition is 


2 

J 


df 

d 

Sf 1] 

dy 

dx 

{dy')\ 


r\ (x) dx = 0 . 


(6.13) 


From our initial theorem, we now conclude that 

df d { df\ 

-— — = 0. 

dy dx \ dy ') 


(6.14) 


This is the Euler-Lagrange equation. It represents only a necessary condition for an extremum. 
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Example 1 

Consider the problem of finding the shortest distance between two points in a plane. 


y 



(* 2 ,y 2 ) 


ds 


*Jdx 2 + dy 2 


x 


Solution 

The distance along an infinitesimal distance ds along a presumed path can be written 
ds = ^/dx 2 + dy 2 = ±Jl + y' 2 dx . 


We will choose the + sign to go along with the above picture where x 2 > x 2 . The total 
length of the curve is 

x 2 x 2 _ 

L = ds = f ^/l + y'^ dx . 

XI XI = f 

d £ 

Notice — = 0 ,, so the E-L eq- is simply 
dy 

dx \dy y 


or 


df 

— = C, , 
dy' 


where C7 is a constant. Therefore 


C 2 = y 


1 + y 


y = 


± c. 


l / 1 - c l 2 


= c. 
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This is just the equation of a straight line of slope C 2 
y = C 2 x + C 3 , 

where C 3 is another allowed constant. In order to find this specific line, set 
y 1 = C 2 X! + C 3 , 


y 2 = c 2 x 2 + c 3 , 


x + za - y*. 

\X 2 - xj x 2 - x x 


Of course, the extremum corresponds to a minimum since there is no maximum distance. 
The choice of y as the dependent variable and x as the independent one above was 
completely arbitrary. 


In the above example, you may have noticed a slight indiscretion. I overlooked the possibility 
of certain paths in formulating the integral. Thus, consider the problem of finding the 
shortest path for the situation: 


y 



x 


Download free eBooks at bookboon.com 


195 






MODERN INTRODUCTORY MECHANICS 


CALCULUS OF VARIATIONS 


Our formalism assumes that the stationary path is single-valued and can be written as 
y = y (x). We can not characterize any path in this case in such a manner. One’s first 
impulse is to simply reformulate the problem so that x = x(y). However, a more general 
solution to the problem is to consider parametric representations of the paths, x (t) , y (t) , 
where t is a parameter. 5 We will consider such path representations later, which we will see 
involve characterizing x and y as dependent variables and t as the independent variable. 


Integrated Form of Euler’s Equation 

3 f 

It often happens that — = 0 . In this case there is an integrated form of the E-L equation 

dx 

that is easier to use. In applications to mechanics, this integrated form will yield the energy 
equation rather than Newton’s equations. Consider 


£ Y 


df 

3y\ 


= y 


df 

dy 7 


, d 
+ y — 


df 


dx \dy 


(6.15) 


However, by the definition of the total derivative 


df 

dx 


df 


df df 
— + — y + — 

dx dy dy 


7 y 


(6.16) 


Thus 


A y.il = A + y . _A + A Ml 

dx \ dy y dx l dy dx l^dy 


(6.17) 


However, if the E-L eq-s hold we have 


df 


dx \ dy 


df 

dx 


df 

f - y - = const. 

3y' 


(6.18) 


This yields a first order differential equation whereas the E-L equations give a second order 
form. The form of this equation when there is more than one dependent variable will be 
investigated in a HW problem. 


Brachistochrone Problem 

As an application of this, we consider a problem which originally gave rise to the calculus 
of variations, the brachistochrone. (Johann Bernoulli posed it in 1696; it was solved by 
Newton within 12 hours of his receiving it!) 
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y 



We wish to find the curve that minimizes the time necessary for a bead attached to a wire 
to go from ( 0 , y i) to (x 2 ,0 ) under the influence of a constant gravitational field. First, 
we have to construct the integral we wish to vary. 


dt 


ds 

V «- mg2 of velocity 



Brain power 


% 


By 2020, wind could provide one-tenth of our planet's 
electricity needs. Already today, SKF's innovative know¬ 
how is crucial to running a large proportion of the 
world's wind turbines. 

Up to 25 t of the generating costs relate to mainte¬ 
nance. These can be reduced dramatically thanks to our 
^sterns for on-line condition monitoring and automatic 
lunhtation. We help make it more economical to create 
clearfflfcdneaper energy out of thin air. 

By snaring our experience, expertise, and creativity, 
industries can boost performance beyond expectations. 

. Therefore we need the best employees who can 
meet this challenge! 

The Power of Knowledge Engineering 


Plug into The Power of Knowledge Engineering. 
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T s r ^ = r Vdx 2 + dy 2 = J dx VT+ y 

J -rj J TT J 


V 


V 


V 


Conserve energy: i w’ + mgy - mgy, , 


v - V 2 9 (yi - y)' 


1 7 11 + y 

=* t -— S' 


' 2 


v 1/2 


V 2 g o Ui - y, 


dx . 


=f[y,y 1 ;x] 


df 


Notice — = 0 .. Use 2— form of the E-L equations: 
dx 

df 

f - y — = C, , 
dy' 1 


1 + y 2 

ji - y. 


v 1/2 


((y. - y)(i + y 2 )) 


- r 

1/2 '"'I 


b ~ yi + y 

' 2 t-j 

y = - 

yi - y 


y = ± 


2a - yi + y 


r i 

f l) 

I • 1 

r ■ V 


v yi - y 

f +: positive slope, negative slope 


Separate variables: 


,1/2 


dx = ± 


yi - y 


,2a - Yi + y) 

{yi - y) d y 


dy , 


/2a (y x - y) - (y x - y) 

Define z = Yi ~ Yr 


v zdz 

x = + J 


0 V 


2az - z 


(6.19) 


( 6 . 20 ) 

( 6 . 21 ) 


( 6 . 22 ) 


(6.23) 


(6.24) 
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The top sign still implies a positive slope, the bottom means a negative slope on the trajectory, 
Notice the denominator is only real if 


0 ^ z < 2a , 

Yi - 2a < y < y, . 


(6.25) 


Also notice that y can turn negative (go below the x-axis) if 2 a > y 1 . Let us introduce 
z = a(l - cos 0') , ( 6 .2 6 ) 


which respects the above inequality in z. We then have 

±1 

x = +a f d0' S ~‘~ n ^ (l - cos 0'). 

J sin 0' 


(6.27) 


Notice the absolute value above and the fact that the upper limit satisfies yi - y = a(l - cos0) 

If the upper limit, 9, is such that 0 < TC, then sin0'> 0 and we choose the root for 

the negative slope, l^Z < (and positive x): 
dx u 

x = a(0 - sin 0). 0 > n 

Note that if 0 > TC then sin0'> 0 for jt < 0'< 0. However, it is at 0'= n that an 
“undershoot” in the trajectory occurs (see later discussion) and the path must have positive 
slope, ^Z > q which means we have two changes in sign. Thus, requiring x > 0 for all 
0 valuefrcsults in the representation: 

z = y 1 - y = a(l - cos 0) ,] 

. . . [ (6.28) 

x = a(0 - sin 0). | 

This is a parametric representation of the path. 0 is just a convenient parameter in this 
problem; however, we will get a physical interpretation of it in a moment. 


To build in the endpoints, we need to solve 

x 2 = a(0 - sin 0) ,] 
y 1 = a(l - cos 0 ) .j 

to find 0 (the final value of 0) and "a". One way to do this is to take the ratio, 
0 - sin 0 


(6.29) 


x. 


Yi 


1 - cos 0 


(6.30) 
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solving for 0 by numerical means and then substitute in either of Eqs.(6.29) to find "a" 
. Eqs.(6.28) are the equations of a cycloid, which is generated by rolling a hoop of radius 
a without slipping along the underside of the line y = Yi- 



I mentioned above the possibility of an undershoot along the path, where y < 0 . This is rather 
nonintuitive behavior, given that we are looking at minimum time paths to travel from one 
given position to another. In a sense, we have traded a longer path length for a greater average 
velocity along the path caused by gravity, the extra velocity “kick” occurring on the initial part 
of the motion. The same sort of phenomenon happens to high thrust rockets launched from 
Earth into low circular orbits. Then optimal paths often produce an overshoot in the trajectory 
and actual path rises above the final desired orbital radius, the extra “kick” happening on the 
final part of the motion. In this case also we are trading a longer path length for a gravity 
“boost” (free from atmospheric drag) given by the downward inserion into orbit. 



You will examine in a HW problem the precise condition under which an undershoot occurs 
in the brachistocrone problem. It may be helpful at this point to examine the undershoot 
cycloid path in more detail. 
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Example 

Compare and contrast the time necessary to reach a point x 2 » Y 1 (causing an undershoot) 
two ways: a) a straight line path and b) a cycloid path. 

Solution 

The two paths look like the following. 


y 



x 


On the straight line path, we have 


1 . 2 

x = — a v t , 
2 x 


a 


x 


g sin a cos a , 


so the total time is 



2 


r 


sin a 



ITT r cos a 




The time necessary for the cycloid path is given by 



However 
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From (6.28) 


dx = a(l - cos 0)d0 , 

(yi - y) = a(l - cos 0) , 



From this we conclude that increments in the parameter 0 and time, t, are related by 


d0 = /— dt. Since x 2 » yi, we notice that the cycloid is almost complete, so 0 — 2 jx 

V a 

and x 2 - 2jta , 



Thus, we have 


T c x 2 »Yi 



» 1 . 


and the cycloid path is clearly superior. Thus, although the undershoot path goes “out of 
it’s way” in dipping vertically, it acquires a larger average x-component of velocity by the 
extra fall. 


The Case of More than One Dependent Variable 


Often problems have more than one dependent variable. We can characterize such integral 
problems in the following form. 

r i 


j - ft 

xi 


Yi( x ), y 2 (x),... , Yl (x), y 2 (x),... ; x 


dx . 


n variables 

One constructs a set of comparison functions ( i = 1, . . • , n) 


(6.31) 


Yi(a, x) - y ± (x) + ar] i (x) ,] 

0i( x i) = 0i( x 2 ) = 0- J 


(6.32) 
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We proceed as before using the chain rule: 


dJ | 

da l“=° 


A2 11 

- D* 2 


df dy ± df dy ± 


dy ± da dy L da 
T ii( x ) ’1i( x ) 


(6.33) 


Do the integration by parts in the second set of terms: 


dJ 


— = fdx y 

da |a = 0 J A 


df 


d df 


~ — T-rhi(x). 


\dy± dx dy i 


(6.34) 


Given that the r|i(x) represent n independent arbitrary functions, one can use reasoning 
similar to the above to conclude, 


— - — -0 

dy L dx \ dy L 


i = 1,... , n. 


(6.35) 


Example 

Find the shortest path between two points on a plane again using both x and y as 
dependent variables. 
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Solution: 

We reformulate the problem as 
ds = ^/x 2 + y 2 dt , 


where x = - and "t" is a parametric variable like 0 in the brachistochrone problem. 

dt 

It is not related, to time. (T here are no dynamics here.) 

L = J yj x 2 + y 2 dt . 

fcl =f[x,y;t] 


The E-L eq-s are 


df d_(df 

dx dt \ rJx 


dt 


x 


'J 


x 2 + y 2 


= 0 , 


= o => A 

dy dt \dyy dt 


V 


x 2 + y 2 


= 0 


From the above 


x 




x 2 + y 2 


Vk 


x 2 + y 2 


= c lf 
= c 2 . 


2 2 

where Ci i2 are constants. From combining them we find that Ci +C 2 = 1 and therefore 

that 


x = (cos 0)t + C 3 , 
y = (sin 0)t + C 4 , 


is the general parametric solution. This way of doing the problem is as easy or easier than 

JT 

our first solution. By taking 0 = — we see that the solution mentioned earlier, x = x o r 
is indeed an extremum path. 

The Case of More than One Independent Variable 

One can encounter a problem also with several independent degrees of freedom: 
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f 


y ( x a x 2 ) 


dy 

dx 1 

t 



x 





t 


partial derivatives 


(6.36) 


Here, I am showing two independent degrees of freedom, but the argument will be general. 
Then 


J 


ff dx, dx 2 f[...]. 


(6.37) 


Now introduce 

y(a, x lf x 2 ) - y(x lf x 2 ) + ar|(x lf x 2 ) , 


(6.38) 


where r| ( x i f x 2 ) — 0 everywhere on the boundary C : 



As usual 


dJ | 

da l“=° 


ff dx, dx 2 

A 


df dy y 
dy da — 


df 




da 
_ dn 
5x t 


(6.39) 


This can also be written more compactly as 
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— I = ff dXj dx 2 

da la=0 


Of - - 

- T| + V ■ V T| 

ay 


(6.40) 


0 f 

where v, = — - -, v = —- and V i s the gradient operator from Ch.l. The 


dy 

Ox, 


dy 

dx. 


*■1 / \ “"2 ) 

two-dimensional form of the Gauss law is 


f dx : dx 2 V • v = f dd v ■ n , 


(6.41) 


where n is the outward unit normal and d i is an infinitesimal element of the perimeter 
(see the above figure; this is actually equivalent to Stokes’ theorem in two dimensions). This 
gives the two-dimensional form of integration by parts if we let v —> r| v , 


f dXj dx 2 v ■ V r| = -f dXj dx 2 r| V • v + (fj df r| v ■ n . 


(6.41 A) 


Apply (6.41 A) to the second term in (6.40) to obtain 

'df 

dx 1 dx 2 

A 


dJ 

da 


n = ff dX l dX 2 

!a=0 


dy 


T| - T| V ' V 


(6.40A) 


since the function r|(x 1 ,x 2 )=0 everywhere on C. More explicitly, 


dJ | 

da l“ =0 


ff dX l dX 2 

A 



/ \ 


Of Y d 

Of 


dy " dx, 




i V 3x JJ 



(6.42) 


Since ri (x x , x 2 ) is arbitrary in A, by setting 

doc 

of our initial theorem we may conclude 


0 and with a suitable generalization 


/ \ 


Of 

dy 


2 



df 



0. 


(6.43) 


For n independent variables, this argument easily extends to show the upper limit in the 
second term is n. 
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Example 

Extremize the integral, 


J = ff chq dx 2 


dy 

dx, 


dy 

dx. 


=f 


dy_ dy ■ 

Y 'dx 1 'dx 2 ' x ^' x 2 


Solution: 

The E-L equation for one independent and two dependent degrees of freedom reads 


df 

dy 


d 

dx 


df 


df 


• if 

\dx 1 


dx 2 1 
2 d 


= 0 


dy 

\dx 2 


We have 


df 

— = 

dy 


df 


dy 
dx , 


= 2 


dy 
dx. 


df 


dy 

,3x. 


= 2 


3y 

dx. 


We obtain 


dy dy 

Hr + -^r = °* 

dXj dx 2 

Notice these terms form the two dimensional Laplacian, so 


V y = 0, 


is the E-L equation. 

Obviously, one can have cases where a combination of more than one independent as well 
as dependent degrees of freedom occur in a problem. Then 


f = f 


Yif y 


2' 


n dep. variables 


dy x ' dy± 

dx, dx. 


dy 2 dy 2 


dx, dx. 


p of them 


x lf x 


2 < • • • 


p of them 


p indep. variables 


n of these groupings 
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In the case of n dependent and p independent degrees of freedom, the E-L equations are 
simply 


df 

dy ± 



j-i 


df 


dx , 


9 (^, 

dx 3 


)= 0, i=l,2 


• r 


n. 


(6.44) 


In the above, “ i" is a free index and “ J" is a dummy (summed) index. In applications 
to mechanics, there is usually only a single independent variable, the time. 


Constraints 

It often happens that we wish to minimize the value of some integral in a situation where 
there are more than one dependent degrees of freedom, which, however can not vary 
independently. In other words, there is some constraint connecting variations of two (or 
more) dependent variables that appear in f. There are three ways of dealing with this situation: 

1. Algebraically eliminate one variable using the constraint condition. 

2. Reformulate the problem so that the variations in the dependent variables 
become truly independent. 

3. Use Lagrange multipliers. 

Situation 1 occurs when one can actually solve the equation of constraint (when it’s linear, 
for example) and eliminate one or more variables. Thus, the constraint is rather trivially 
removed. 

Situation 2 usually means to change to a new coordinate systems. An example would be 
the motion of a bead on a hoop. 
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Regarding x and y as dependent variables, it would be foolish to use a rectangular coordinate 
system since 


2 . 2 

x + y 


R 




and variations in x and y are not independent. It would be far superior to change to 
cylindrical coordinates. Then the only dependent variable is the angular coordinate (|), for 
which there is no constraint. (It is understood that the independent variable in this situation 
is t, the time.) 


As an example of situation 2, in the case of geometry rather than dynamics, let us find the 
shortest distance between two points on a sphere. 



X 


In general it takes 3 coordinates to locate a point in 3 dimensions, but since our point is 
given to be on the surface of a sphere of radius R, we will need only 2 coordinates, which we 
will take to be the spherical coordinate angles 0 and 4>. The constraint has been eliminated. 

The length element in these coordinates is 

ds = R (d0 2 + sin 2 0d(j) 2 ) . (6.45) 

Without loss of generality we can write 

ds = R [0 ,2 + sin 2 0 ] | dcp | , (6.46) 
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d0 

where 0' = — • (|) is now the independent variable, 0 is the dependent. Identify 
d(|) 


:[0, 0' ; (j)] = (0 ,2 + sin 2 0) 1/2 . 


(6.47) 


df 


Since — = 0 , we may use the second (integrated) form of the E-L equations: 


(0 ,2 + sin 2 0) 1/2 - 0' 


( 


0 ' 


\ 


(0 ,2 + sin 2 0) 


= a , 


sin 


0 = a (0' 2 + sin 2 0) ' , 


(^_ = _ 

\dd) sin 2 0(sin 2 0 - a 2 ) ' 

d0 

d(j) = 


sin Qy-j sin 2 0-1 


(6.48) 
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We will carry the ± sign with us to see the effect. We can now write 


d(j) = 


sin 2 9 


d0 


esc 0d0 


sin 0 



t 

csc 2 0 = 1 + cot 2 0 


d(j) = 

t 


P" 1 — cot 0 
d0 

yf: 1 - p“ 2 COt 2 0 


d . _i 

^cot 0^ 

d0 

1 P JJ 


d0 . 


(6.49) 


d cot 0 
d0 


= esc 2 0 , p = - 1 • 


In this form we see that the right hand side is a perfect differential, and thus 

_J cot 0^ 


4» = ± sin" 


P 


+ b , 


(6.50) 


where "b" is another constant. Notice that there are now two constants, b and P, 
available to fit two points on the sphere. To make this result look more familiar, take the 
sine of both sides and use 


sin (4> - b) - sin <j> cos b - sin b cos 4> , 

cot 0 


sin 4> cos b - sin b cos (j) = ± 


P 


=> R sin 0(sin 4> cos b - sin b cos 4>) = 

Identify ^ cos 0 = z , 

R sin 0 cos (J) = x , 

R sin 0 sin cj> = y , 


R cos 0 

P 


y cos b - x sin b 



(6.51) 


which is the equation of a plane passing through the origin (x, y, z = 0 ) . The sign on 
the right of (6.51) may be absorbed into a redefinition of P and adds no extra generality in 
defining the plane. The intersection of the plane and sphere forms a “great circle” connecting 
the two points on the sphere. There is still a choice of path (long way or short way) which 
connects the two arbitrary points along this great circle. This is a case where both extreme 
values are uncovered by the variational principle. 
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Lagrange Multipliers 

Let me begin to explain what is meant by the use of Lagrange multipliers, as in situation 3 
above. I’ll start in a context of functions rather than integrals. Let us say we wish to extremize 
a function of more than one variable, such as f (A,B) ; but subject to some constraint 
connecting the variables A and B so that their variations are not independent. We assume 
the constraint can be written as (called “holonomic”; we’ll see more of this next chapter) 

g (A, B) = 0. (6.52) 

Examples: 

1. A 2 + B 2 = R 2 => g(A, B) = A 2 + B 2 - R 2 

2. A = CB + D => (A, B) = A - CB - D 

The conventional approach to this problem is to try to eliminate one variable in terms of 
the other, as in 

g(A, B) = 0 => f(A, B(A) ) , 


and then apply the zero slope necessary condition 


df | 
dA l A ° 


0 . 


The Lagrange method avoids the explicit elimination of variables by using the constraint 
equation as a subsidiary condition. In the following, it is understood that all variations are 
being evaluated at the point or points of extremum. 


The total variation of f is given by the chain rule, 


df 



dA + 



'df 

dA 


df dB 
dB dA 


dA . 


The extremum condition is 


(6.53) 


df = 0, 

df df dB 

dA dB dA ’ 
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Now from g(A,B) =0 we have 


0 = dg = 



dA + 



r 


dB = _/ag\/3gV 
dA UaJUb/ 


df _ _ df (_ ( dg\ (dg)- 1 ' 
dA dB [ UaJUbJ )' 

df /dgT 1 I _ 3f (dgY 1 i 
dA\dA,J I a o' b o dB^i l A o' B o 


(6.54) 


Again, the variations are being carried out at the extremum, and to emphasize this I am 
indicating this point explicitly as A 0 , B 0 in (6.54). I have also made the conventional 
definition of the Lagrange multiplier, A, in a variable-symmetric way in the above. It is 
just a constant. Writing out the two equations contained in (6.54) separately and explicitly 
stating the constraint, viewed as a subsidiary condition, we have 


df . dg n I 

— + A — = 0 

dA dA 

df . dg 

dB dB ' 

g(A,B) = 0 , 


(6.55) 


which is a complete characterization of the problem. Notice we can write these equations 
a bit more elegantly if we define 

f * - f + Ag . ( 6 . 56 ) 

Then the conditions 


df* 


dA 

df* 


dB 



= o, \ 


g(A, B) = 0 . 


(6.57) 


look like a normal, unconstrained, extremization problem, except for the last explicit constraint. 
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Example 

Find the rectangle of greatest area which can be inscribed in a circle of radius “a.” 


y 


x 



Solution: 

The area is 4xy and the constraint is x 2 + y 2 = a 2 • Therefore 
f(x,y) = 4xy, 


g(x,y) = x 2 + y 2 - a 2 = 0 


We have 


df dg 

— + X — = 4y + 2X x 

dx dx 

df dg 

— + X — = 4x + 2X y 

dy dy 


2.2 2 

x + y = a . 


= 0 , 
= 0 . 


(E.l) 

(E.2) 

(E.3) 


Combining (E.l) and (E.2), we find x 2 = y 2 . Then using (E.3), we obtain 
x = y = . This just gives a square. 

•J2 

This method can be generalized to m variables, x i5 and n constraints as 
df*(x lf ... , x m ) 


dx, 


= 0 , i = 1,... , m 


gj(x lf x 2 ,... ,x m ) = 0 , j = 1, ... ,n 

n 

where f* = f(x lf x 2f ... , x m ) + J ••• * X J • 


(6.58) 

(6.59) 
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How does this general method apply to the situation with integrals and functionals rather 
than functions and variables? Let us consider the simplest case of 2 dependent and one 
independent variable (like the x,y position of a particle at time t). We have 

f = f[yir y 2 r y'h y 2 ; x ]- 

As usual, 


Yi (a, x) = y, (x) + arp (x) , 

y 2 K x ) = y 2 ( x ) + «h 2 ( x ) f 

hi ( x i, 2 ) = h 2 ( x i, 2 ) = o. 

Then, skipping some hopefully familiar steps, we have 


(6.60) 


dJ 

da la= 


A 2 

= f dx 
0 j 


{ Of 


df 


\\ 


6y x dx \dy i;/ 
But now assume that the variations 


hi ( x ) 


df 


df 


dy 2 dx 1 dy 


h 2 ( x ) 


2 // 


d Yi 

da 


hi 


dy 2 

da 


h 


2 r 


• (6.61) 


are no longer independent because of a constraint of the form 


g(yi,y2,*x) = o . 


(6.62) 


However, we can write 

-I -?<* 

da la=0 z 


df 


df 


dy x dx \5y, 


( 


Of 


df 


dy 2 dx l dy 


h 2 


2 // 


hi 


hi( x ) 


Requiring this to vanish and since r\ x (x) is arbitrary, we have 


df 


df 


df 


dy 1 dx \3yj7 \dy 2 d.x\dy 2// 
Since g(y 1 (a,x),y 2 (a,x);x) = 0 , we have 


df\\ T| 2 (x) 


hi(x) 


(6.63) 


dg 

da 


dq dy ± 
dy 1 jhx^ 

=hT 

h 2 ( x ) 

hi( x ) 


dq dy 2 
dy 2 ^da^ 


3g V dq 


= 0, 


dyj\dy. 


(6.64) 
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Used above to eliminate this ratio, we find 


r dt 

d ( df \1/ dg ) 

r df 

d ( df \1 

( 

|_3Yi 

>1 

CO 

-1 

<T2 

1 M 

1 T3 

1 

[dy 2 

dx Uy 2 J_ 

UyJ 


(6.65) 


Set both sides of (6.65) equal to -X(x) , an unknown function of x. Our three equations, 
representing a complete characterization of the extremum problem, are now 


df 

dYi 


df 

dy 2 


_d_ _df_ 
dx {dy[ 

d ( df ' 
dx 1 dy 2 


V ' dy, 

+ x(*)*L 
v ' 3y 2 


g (yif y 2 ; x) = o. 


o, 

o, - 


( 6 . 66 ) 
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This represents 3 equations in 3 unknowns: YirY 2 r and X. Again, we may simplify these to 


«1.A «1, .o. 

dyx dx \dy 1 

df* _d_ ( df* _ 
dy 2 dx \dy ' ,f 


(6.67) 


g(yif y 2 ; x ) = o , 


where f* = f + kg. The general case of m dependent variables, yi(x) , and n 
constraints can be written as 


df* _ d (df* 
dYi dx \ dy ± 


i = 1,, m 


( 6 . 68 ) 


gjfyir ••• f y n ; x ) = 0 »■ j = 1 > ^ , 


(6.69) 


where 

n 

f * s f + 2 M x )g k • (6.70) 

k = 1 

The above represent m + n equations in m + n unknowns. 


We will see in the next Chapter that the generalized E-L equations, with an appropriate 
choice of f, are equivalent to Newton’s force equations on a particle. Therefore, it is not 
surprising to learn that the Lagrange multipliers, L k (x), often have the physical interpretation 
as the forces of constraint. We will leave examples of the use of these equations to the next 
chapter. 


Isoperimetric Problems 

I can not leave the subject of the calculus of variations without mention of the subject of 
so-called isoperimetric problems, even though we will not be using this in the Chapters 
which follow. Isoperimetric problems are also constraint problems, except the constraint is 
in the form of another, fixed, integral rather than a relation between dependent variables. 
In it’s simplest form, the problem is to find an extremum of 

x 2 

J = Jdx f[y,y' ;x], (6 . 71) 

Xl 
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subject to the integral 

x 2 

L = f dx g [y, y ' ; x] , (6.72) 

X 1 

having a fixed value. Typically, L represents a perimeter on a given curve. Here, the use of 
Lagrange multipliers is also useful. The formalism is quite similar to the use of Lagrange 
multipliers in the function/variable case studied earlier. 

Let us introduce a path varied about the extremum with independent arbitrary functions 
Y|(x) and k( x ) using a lf and a 2 as parameters. Then define 

y(a lf a 2 , x) - y(x) + a 1 rtfx) + a 2 £j(x), (6.73) 

where at the endpoints 

x 1 , 2 ) = ?( x i, 2 ) = 0 r ( 6 -74) 

as usual. Then J and L become functions of cti and a 2 , 

J -» j(a lf a 2 ) , 

L -»• L(a lf a 2 ) , 


and the sought extremum is at = a 2 = 0. At this point, the formalism becomes 
identical to that studied before, with a function of two variables connected by a holonomic 
constraint between the variables. We can follow the same reasoning as before to arrive at 

dJ ( dL \ _1 I _ dJ_ l dL \ 

da 1 \3aJ l a i,2=o da 2 \3a 2 J 

where I have defined the Lagrange multiplier, A,, which is again a (unknown) constant. 
Writing part of (6.75) out in full tells us (the other part gives the same information) 


= -A 


(6.75) 


dJ 

3a 1 


+ A 


dL 

da 1 


= 0 , 


(6.76) 


=> 



A(2L]) + 

dx ^<9y ' )J \^6y 


d 

dx 



T|(X) = 0 . 


Since (x) is arbitrary, we conclude that 


(6.77) 


T_jL7L).o, 

dy dx \dy ') 

f* s f + Ag . 
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X 2 

After y (x) is found, the condition L = J* dx g then fixes one of the parameters in the 

Xl 

solution. There is an interesting duality in such problems in that if we extremize L and fix 
the value of J, reversing the roles of f and g we would get the exact same equations as 
above and thus the same solution. 

Thus the constraint in isoperimetric problems is built in simply by adding to f the constraint 
condition multiplied by X, and varying f * in the E-L equations rather than f . The general 
rule we are seeing is: to build in a constraint, add the constraint condition, multiplied by 
an unknown multiplier, and vary the result as if there were no constraint! 

Example 

Find the largest area enclosed between the x-axis and a string of fixed length attached at 
it’s ends to points x x and x 2 on the x-axis. (“Dido’s problem”) 



Solution: 

We have for the area and the constraint, respectively, 


X 2 

J = J dx y , 

Xl 

X 2 

L = J dx -J 1 + y 2 . 

Xl 


We identify 

f * = y + X-Jl + y 2 , 


and we write 
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— (y + X-/T 

dy ' 


+ y’ 2 ) - 4: i(y + W 1 + y 2 ) = o, 


dx dy' 


i/l 


— w- + y / 

df* 

(One can also use the integrated form of the E-L equations here since - 

work necessary either way is almost the same.) Integrating (E. 4 ) , we find 


1 

X 


(E.4) 
0 , but the 


x 


V 1 + y 2 *■ 


= T + C, 


where C is a constant of integration. Now, solving for y 


one finds 


x 

X 


+ c 


11 ' 1 f + c 


where a ± sign from taking the square root is being kept. Separating variables and introducing 
x 

z = — + C , 

X 


we 


find 


. X r Zdz 

y - ±A J 


i 

which integrates immediately to 


x 


Y = ± Kl 1 ~ 1^ + c l + K , 
where K is another constant of integration. This can be rearranged, after squaring, to read 
(y - k) 2 + (x + CX) 2 = X 2 . 


Thus, the solution is always a section of a circle. In the case above, it might look like 
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The three constants in the solution, K, C and X, are fixed by x lf x 2 and the length 
of the string, L. Obviously, in this case our extremum corresponds to a maximum. Also, 
obviously, there is no solution to the problem when L < |x 2 - xj ! 

Variation of the End Points of Integration 

The last thing I want to discuss in this Chapter is the subject of varying endpoints in our 
integrals. Suppose the limits of the integral 

x 2 

J = f dx f[y, y' ; x], (6.79) 

Xl 

were not fixed, but were allowed to vary as well as the “interior” points, as shown. 



Obviously, T| (x) no longer vanishes at the endpoints. Although this seems like an unlikely 
situation, we will see in the next Chapter that endpoint variations yield conservation laws. 
Note that the infinitesimal endpoint variations occur on given curves at the two boundaries 
of the integral. 
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Let us evaluate the variation in J. We have 6 


dJ | 

ckT 1“-° 


d 

da 


x 2 (a) 

f dx f[y(a, x), y 1 (a, x); x] , 

xi(a) 


(6.80) 


where we are making the endpoints change with ct, as in the above picture, and where it 
is understood the right hand side terms are evaluated at a = 0. We find (making use of 
the Leibnitz rule since there is a dependence in x lf2 ) 


dJ | dx X2 X2 r df 

- _ n = f - + I dX - , 

da ' a da J da 

X 1 XI 


(6.81) 


The second term on the right in (6.81) is the usual term we get when the endpoints are 
fixed. We have 


df _ df dy df dy' 

da dy da dy ' 

V "rf" 


(6.82) 
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and doing an integration by parts on the term involving r|' yields 


dJ 

da 


- f ^ + „^ | 

a=o ( da dy ') Xi 


,/dx 


df 

d 

1 df ^1 

dy 

dx 

-1 

>1 

<“C> 


il(x) . 


(6.83) 


We now require 


dJ 

da 


a = 0 


= 0 , 


(6.84) 


and assume that the E-L equations are satisfied in the interior so that the integral in (6.83) 
vanishes. This results in the endpoint condition: 


' dx df 

f- + in- 

v da dy ' 


= 0 . 


(6.85) 


The endpoint requirement looks arbitrary because h appears there. However, the value of 
T at the endpoints is actually determined through the given endpoint curves, 


so that 


dy (a, x) . / dy(x) dx dii dx 

---- = — + r\ + a —- — 

da la=0 \ dx da dx da 


dy dx 

q --y — 

da da 


Multiplying by da, the endpoint requirement above assumes the form 

= 0 . 


( 6 . 86 ) 


(6.87) 


/ , df \ 

df , 1 

f - y — 7 

dx + -dy 

dy ') 

dy' j x 


( 6 . 88 ) 


This is called the “transversality condition.” The slope at the endpoints determines the ratio 
of the variations of dy to dx; see the figure below. 
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In the case of more than one dependent variable, the second and third terms in (6.88) 
become summed over all the dependent variables, y ± . For example, when there are two 
dependent variables, y and y 2 , we obtain 


<*-2 



)dx + 



0 


(6.88a) 


There are many special cases and re-writings of (6.88) that can be considered. In cases where 
the first point is fixed and only x 2 is varied, we have 


, df \, df , 

f - Y —- dx + —- dy 


= 0 . 


dy 'j dy 

If in addition one has dx=0 at x 2) the condition becomes simply 


(6.89) 


df 

dy' 


= 0 , 


(6.90) 


x 2 


and if dy=0 at x 2 , 


f - y' 


df \ 
dy 'J 


= 0 


(6.91) 


Example 

Find the path of least time for a sliding bead from the point ( 0 ,yi) to any point along 
the line x—x 2 . 
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Solution: 

First of all, we know that the fastest path from the given point to any given point on x=x 2 
is a cycloid. The question is, which of these cycloids yields the smallest time to arrive at 
x=x 2 ? Since a single endpoint, x 2 , is being varied for which dx=0 , the endpoint condition 
is just (6.90) above. We have from our earlier example that 


f = 


1 + y 2 

vYi - y, 


, 1/2 


and we obtain 


df y ' 

dy ' ,J{yi - y) (i + y 2 ) X2 
=> y'(x 2 ) = o. 


Therefore, the path which has zero slope at x=x 2 i s the winner. Paths above this one have 
a shorter distance to go but acquire a smaller a smaller average velocity. Paths below this 
have greater magnitude of velocity by dipping lower in y, but have longer path lengths. 
Here, an “undershoot” solution never occurs, or rather, is always on the verge of occurring. 


CHAPTER 6 PROBLEMS 

1. Let’s say we had an integral of the form 

x 2 

J = = J dx f[y", y', y; x] . 

x i 

Find the equation that f must satisfy for J to be an extremum. Only consider those 
paths y (a , x ) = y ( x ) +a r|(x ) such that 

il(xi) = ri(x 2 ) = 0, 


and 


H ’ (xi) = r| ' (x 2 ) = 0. 
(assume r|(x) is twice differentiable.) 
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df 


2. For = 0 , show that the integrated form of the Euler-Lagrange equation 

becomes (C is a constant) 


df 


f - 2 yi dy[ = C ' 


i=l 


when f is now a function of n dependent variables, 


f = f lYi'Yi' • ■ ■ 'Yn'Y'i ’Y: 


■y; ; x ] 


[Hint: Proceed as in the text, but use the chain rule.] 


3. Given 


-*•2 

J = f dx f[y, y';x], 


show that the Euler-Lagrange equation may also be written as 

(df) 


d 2 f 


+ Y 


d 2 f 


K dy U J “ \dydy’ 

4. a) Extremize the integral 


+ 


d 2 f 

dxdy’ 


,dy t 


= 0. 


J = J dx [ay' 2 - by 2 ] , 

0 

where a and b are positive constants, to find y(x) such that y(0)=0, y( 1) =1. 
b) Find the value of the integral J using the stationary path y(x) found in a). 


5. Find the y(x) that gives a minimum value of 


I = 


f dx(y ' ) 2 

_0 _ 

jdxy 2 


with y ( o )-y ( 1 )-0 . [From Ch.2 problems, “Variational Mechanics for Engineers” 

by M. W. Wilcox (my Father). Hint: What would you do if you took the derivative of 
g (x) 

f ( x > = h(x) ? ] 
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6. a) 


What is the initial slope, 


, of the brachistochrone solution? 
dx 



X2 

b) Under what conditions on the ratio — will the cycloid dip below 

(undershoot) the point at (X 2 ,0 ) before reaching the point (X 2 ,0 ) ? [Hint: 
The lowest point on the cycloid occurs when 0 = it and an undershoot will 
occur if the y value is negative there.] 


7. A space station which is rotating at an angular velocity 0) about an axis through 
it’s central hub has a inter-spoke mail delivery system in which packages are sent 
down enclosed frictionless tubes from a central hub to various stations along the 
rim of the station. 
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All packages start from rest at the hub along (|) = 0, r = a , then curve along a path 
to reach their final destination at (|)=(t>i, r=R. The connection between speed along the 
path, v , and r is v=oor .(This comes from energy conservation since kinetic energy is 
\ mv 2 and potential energy is \ moo 2 r 2 - ) 


a) Show that the time for delivery can be written (use cylindrical coordinates; 

r' = — ) 
del) 


T = - f d(t> ^ 

a) „ 


r ,2 + r 2 


b) Minimize the integral in (a) and show the path is a logarithmic spiral. 


8. Extremize the 3 dimensional integral (4> = cj> ( x ) ) , 


J = I I I dxidx 2 dx 3 [ — (V 4>) 2 - p(x) <f> ] , 
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where p ( x ) is a given function and V is the three dimensional gradient operator. Find 
the partial differential equation for 4> that results. 

9. Find the path of shortest distance between two points on a sphere of radius R 
as in the text. 



This time treat 0 as a dependent variable and f|) the independent variable. 


10. a) Re-do the sphere as in the text, using 0 as a dependent variable and 4* the 
independent variable. Ffowever, use the usual (nonintegrated) Euler-Lagrange 
equations and show that one obtains 


cos 0 d / 1 \ 

sin 0 dtj) y sin 2 oJ 

. d0 

where 0' = T 7 
d(j) 

b) Show that 


d 

/cos 0\ 

dtj) 

(sin oj 


1 

= - 0 ' : jt: 
sin 2 0 


and use this to solve the equation in (a). 


11. One more time with the sphere. Using both 0 and as dependent variables, 
formulate this situation as a parametic problem. It is enough to simply write 
down the appropriate extreemum conditions in this case. Combine the two E-L 
equations and show that one obtains the same equation given in prob. 10(a) 
above. 
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12. Use the method of Lagrange multipliers to find the area of the largest rectangle, 
with sides parallel to the major and minor axes, which can lie inscribed in the 
ellipse: 


x 
( a 


) + 


<b 


1 . 


y 



13. Find the cubical of greatest volume which can be inscribed in a sphere of radius 
R given that x=2y where 2x and 2y are lengths of two sides of the cubical. 
Formulate this as a Lagrange multipier problem. 

14. a) Find the equation of the shorest curve between two points on a cylinder of 

radius R using cylindrical coordinates, [ds 2 = R 2 d(j) 2 + dz 2 .] 
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b) Let us say we are given an unstreachable length of string, L. Find the 

equation of the line on the cylinder’s surface which encloses the greatest area 
on the cylinder’s surface between the string and a vertical line on the cylinder. 
See the below: 



r2 r2 

Ans.: (z-K) 2 + (X(|>-C) 2 = , where C, K and X are constants. 

(This is just the equation of a portion of a circle pasted on the surface of the cylinder.) 

15. An economical mathematician wishes to make a lampshade with as small 
an area as possible by rotating a given curve in the y-x plane passing through 
( x ifYi) and (X 2 , Y 2 )around the Y-axis. A wire form is used to construct this 
curve, which can be bent but has a fixed length, L. 



Find the form of the equation, y (x) , which describes this lampshade. [Hint: Set this 
problem up as a constraint problem using a Lagrange multiplier. Make sure you tell me 
how the Lagrange multiplier is determined.] 
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16. Prove the following: 

The shortest perimeter of an arbitrary path enclosing a fixed area between the x-axis and 
a curve y (x) , as shown, is a portion of a circle. 



[Hint: This may not be as difficult as you may think.] 

17. A wire of fixed length L is used to connect the points shown. A bead slides 
without friction along the wire. 



We wish to fashion the wire’s shape to construct a path leading from ( 0 , y 1 ) to (x 2 , 0 ) 

which uses minimum time, but with fixed length. We of course have L > \j y^ + • Set 

this problem up as a constraint problem using a Lagrange multiplier. It is not necessary 
to attempt to solve the resulting differential equation for y (x). 

18. Another mathematician now wishes to construct a lampshade of least area 
for the surface generated by revolving about the x-axis a line, y (x) , which 
connects the point at x=0,y=a to anywhere on the fixed wall at x=b, as 
shown. 
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y 


y=a 



/ 


x=0 


x=b 

X 


x 


Show that the equation is 

x-b 

y(x) = C cosh( — ), 

where C (the value of y (b) is determined by 

b 

a = C cosh( c ). 

[Hint: You must consider endpoint variations at x=b.] 

19. Going back to the space station mail problem, we found that the time for 

dr 

delivery was (using cylindrical coordinates; r' = -) 


d(j) 



a) Minimize the time subject to the condition that the length of the tube, L, is 
fixed. Simplify the implied differential equation as much as possible. (There is 
no need to attempt to solve the resulting differential equation, which is rather 
difficult.) 


b) Given the solution to 7(b), show that the fastest mail path to the outside ring 
is the shortest path by considering endpoint variations of the final location. 
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20.Find the angle, a, of the straight-line path for which a sliding bead will pass 
from the origin to the wall at x = x 2 the fastest. Note we have a constant 
gravitational field downward in the diagram. 



/ 

/ 

/ 

/ 

/ 

/ 

/ 


dy 

Find the equation ofthe brachistichrone which has zero slope at the wall ( 
and show that the bead’s time of flight on this path is less than the result found above. 


x=X2 


= 0 ) 


21. Find the endpoint condition on the cycloid which provides the fastest descent 
from a fixed point ( 0 ,yi) on the y-axis to any point on a line of slope s, as 
shown. 



x 


slope=s 


[This generalizes the last example of the text. Answer: the fastest cycloid always approaches 
perpendicularly to the given line.] 
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7 LAGRANGIAN AND 

HAMILTONIAN MECHANICS 


THE ACTION AND HAMILTON'S PRINCIPLE 

For an n-particle system, Newton’s equations are 
d 2 x a j 

m a 2 “ • (3n eq-s) 

at 

Imagine (no external forces): 


(7.1) 



: force of (3 on a. 


Total force on a: 
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Assume that 

U ap 


(3^a 

II no self-force 
results from a two-body potential 
(no a sum!) 


= -v a u ap . 


- u 


ap 


X a - X p 


(7.2) 


(7.3) 

(7.4) 


Notice that (7.4) implies U a( , = U |la , which gives Newton’s third law from (7.3). Eq.(7.4) 
also means that the force between the particles is along the line connecting them. This is 
not the most general assumption. 


Explicitly, the notation in the following is 


x. 


x„ 


particle no.f fcoordinate 


no. 


m a X ai + 


dx r , 


U 


ap 


^P^a 


= 0 . 


(7.5) 


Another way of writing this equation will now be explained. The total kinetic and potential 
energies are 


T 

U 


\ 2 *ai ' 


a < p 



X a - X p 

p 



f sununing over both a, p such 
that a < p. 


(7.6) 

(7.7) 


The reason for the a < (5 restriction in the sum in (7.7) is to count a given interaction 
between two different particles only once. 


Form L = T - U . Notice 



Beginning to look like the second term in (7.5). 
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Take the special case of 3 bodies (a = 1,2,3 ) to see another way of writing (7.8). We have 

(7.9) 


2 U YP = U 13 + U 23, + U 12 , . 


Y <P 


Y < 3 


y < 2 


Let’s now consider the derivative of this with respect to x 1;L : 

d — d /— — — \ 

2 U YP = — K + U 23 + U 12 ) 


dx 


li Y < P 


dx, 


dx 


d — \ 

(u 13 + u 12 ) 


2 u ip- 


( 7 . 10 ) 


0x u 

Do the same thing for a = 2 (to see a necessary property for the rewriting we will do) 


d y 

p 


d 

dx 2i y 4 


dx 2i 

<3 

(u 23 


u 12 ) 

dx 2i 


d 

(u 23 


_ d 

U 2l) 

dx.. 


2i 



j 

d 

V 

U 2P 


dx 2i 

2 j 

(3*2 



(u 13 + u 23 + u 12 ), 


. U Y p = Up Y necessary 


Generalization for an n-body system: 

d ^ - d 


dx. 


2 u .f - 


■ ai Y < P 

Using this above, we get 


^ X ai P*<* 




dL 


0 ( 

2 

P* a 

3x ai 

^ X ai l 

more, 

d 

(dL) 

d 

' dT ' 

dt 


dt 

UJ 


and 


l 2 P# j 


d 

dt 


t2 

' p,j 

dT 


m p4 


= m a X ai f 




m a X ai 


(7.11) 


( 7 . 12 ) 


(7.13) 


(7.14) 


(7.16) 
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Now, comparing with (7.5), we have 
d ( dL \ dL 


dt \ dx r 


3x . 


= 0. (still 3n equations) 


(7.17) 


Reversing our usual logic, this now implies that the integral 


u 2 

S S / dt L ( X ai' t) F 


(7.18) 


(“the action”) is an extremum, or stationary, when the equations (7.17) hold. "L" is called 
the Lagrangian. Eq.(7.18) represents a vast simplification of the situation. All information 
is contained in a single, scalar, quantity. This gives rise to Hamilton’s principle: 


Of all possible paths along which a 
dynamical system moves, the actual path is 
that which makes the action, 

t 2 

S = J dt (T - U) , 

ti 

stationary. 


Notice we have not proven this principle, we have simply constructed the particular form 
of the action for a special type of system. 

Generalized Coordinates 

Let’s now generalize the systems we are describing by introducing new “generalized” coordinates, 
qi, which uniquely specify the positions of the particles. These coordinates are assumed to 
automatically respect the possible constrained motions of the system. Examples are: 
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planar double pendulum {l ir l 2 fixed) 


z 



m i : 4>i 

m 2 : t )> 2 ' 0 / ^ 


String of length fixed L with masses irq, m 2 at ends slipping through a hole in a plane 
(4h not shown). 

Notice that some generalized coordinates have units of length while others do not. Notice 
also the choice is not unique. In writing the generalized coordinates as q i5 I am dropping 
the double indices on particle number and coordinate axis. In general there are fewer 
qi variables than xaj variables because of the elimination of constraints among the x a j. 
Our new “phase space” for such a system is specified by the generalized positions, qi, and 
the generalized velocities, qi. If there are m constraints, then the number of generalized 
coordinates, s, is 3n - m. (For solids the counting goes like s = 6n - m since 
three coordinates and three angles are needed to specify a body’s position and orientation.) 
If there are no remaining constraints among the remaining generalized variables, they are 
called a proper set of coordinates. 
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Let us assume we know how to write the kinetic and potential energies, which are simply 

scalar quantities under coordinate changes, of the system in terms of the q± and q ± That 
is, we know the transformation equations from the x ai to the qj variables. Symbolically, 
we write 



which implies 

x ai = x ai (q jf q jf t) . 


(7.19a) 


(7.19b) 


Notice that the time can explicitly enter the transformation. We assume that T depends 
only on x ai and x ai and U depends only on x ai , so therefore T(x a i) = T(qi, q i ,t) 
and U (x ai ) =U (qi, t) and we may form 

"^2 

s = f dt L (qif <3i? t) • (7.20) 

tl 

By extremizing this expression, we will find the form of Newton’s equations in the new, 
generalized coordinates. To perform this variation and to introduce a new notation, let us 
again reverse our logic and extremize the action in the form (7.20). This new notation is 
easier to use and more elegant than the notation in the last Chapter. 


Let us imagine varying a particular coordinate q k (t) around the stationary path, as shown 
below. (Note the prime on q k (t) does NOT represent a derivative!) 



q k (t) - q k (t) + 8q k (t) , 

= o. 


t 


tl 


t 2 
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8q k (t) is imagined to be an infinitesimal, arbitrary variation about the extremum solution, q k (t). 
This “frees” us from the clumsy use of a and r| in the previous chapter. The connection 
with the previous notation is 

8( > da r: ( >la- 0 - (7.2i) 


Also 


dq k (t) 

dt 


dq k (t) 

dt 


+ 


dt 


6q k (t) . 


Now 



- ( dq 4 

dq k 

. dq k 

dq k 

V dt J 

V dt J 

dt 

dt 1 

dt 


However t ' = t in the variation, so we have 



( 7 . 22 ) 


(7.23) 


(7.24) 
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Tliis means we can interchange the "5" and the “ — " operations. This makes sense since 
the variation defined in (7.21) is independent of t. 


It’s important to realize that the varied paths in general do not conserve energy or other 
quantities. Think of a particle subject to a potential, U(x). The varied paths change the 
relationship between x and x and, therefore, between kinetic and potential energies. Thus, 
the varied paths generally are not physically possible in classical mechanics. However there is 
a formulation of quantum mechanics which uses a quantity called the path integral. In this 
approach, particles are seen to actually experience these “unphysical” paths, the actual motion 
being a sum of the wave amplitudes associated with each of the paths; the path integral 
represents this sum. One can think of the motion as a tube of possible paths surrounding 
the classical path picked out by the E-L equations. This tube becomes thinner in the limit 
that the motion becomes more classical, i.e. where quantum numbers become large and 
the discrete nature of the quantum world becomes almost continuous. The path integral 
approach to quantum mechanics provides a beautiful physical picture which connects to 
and generalizes the ideas of classical mechanics. 


Under these variations, we have the change in the Lagrangian (sums from 1 to S) 


8L - 2 


(al] 

x ( dL ) 

\ dc 2 k j 

6 q k + L. 

l a< *J 


8S - / dt£ 


Sh 

dk 8q * 

| from (7.24) 


<9L d 


dL x — - * 

dq k <9q k dt 


Do an integration by parts on the second set of terms on the right: 


(7.25) 


(7.26) 




*■ h 




t, k ^k dt 

If the 6q k are independent and arbitrary, when we demand 


-/ dt S A — K- 


(7.27) 


SS = 0. we obtain 


—-— f — ) = 0 

dq k dt (dq k ) 


(7.28) 


One can also show (7.28) follows directly from the original set, (7.5), and the transformation 
equations, (7.19). 
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Examples of the Formalism 

Some examples at this point would be helpful. 

Example 1 

Set up the simple harmonic oscillator in the Lagrangian formalism. 

Solution: 

The simple harmonic oscillator has kinetic and potential energies, given by 

1. 2 l l2 

T = — mx , U = — kx , 

2 2 

1 -2 1, 2 

=s> L = — mx - — kx . 

2 2 

The E-L equation is 

S1 -l(*| - 0, 

dx dt \ dx ) 

=> -kx - — (mx) = 0 , 
dt v ' 

k 

=> x + — x = 0 , 
m 

which is the usual equation of motion. 

Example 2 

Do the undamped plane pendulum in the Lagrangian formalism. As the generalized 
coordinate use an angle. 


/////// 



Solution: 

In rectangular coordinates the kinetic and potential energies are 
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T 


im(k 2 + y 2 ), 
U = mgy . 


The transformation equations are 

x - Esin 0 , 
y = E( 1 - cos 0) . 


Thus 

x 2 + y 2 = f 2 (cos 2 0 + sin 2 0j 0 2 = f 2 0 2 , 

and 

L = — mf 2 0 2 - mg t (l - cos 0) . 

The generalized E-L equation is now 
00 dt { 00 ) 

_n 

=> -mgf sin 0 - — (mf 2 0) = 0 

=> 0 + — sin 0 = 0. 

£ 

Tlris nonlinear equation is also quite familiar from Ch.4. 

Example 3 

Set up the Lagrangian equations of motion for the double pendulum. For simplicity, assume 
m x = m 2 = m as well as = £ 2 = l 
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Solution: 

The kinetic 

T 


and potential energies are 


= i m(x 2 + Yl) + ^ m(x 2 + 


U = mg( yi + y 2 ) • 



The transformation equations are 

Xj = i sin 0 X , Yi = -t cos 0 : , 
x 2 = t (sin 0j + sin 0 2 ) , 
y 2 = -i (cos 0 1 + cos 0 2 ) . 

I’ll leave it to you to show in a HW problem, that the kinetic and potential energies now 
become 

T = m( ;2 0 2 + i mf 2 0 2 + m( 2 Q 1 Q 2 cos (0 : - 0 2 ) , 

U = -mgI? (2 cos 0 : + cos 0 2 ) , 


Needed for the 


dL 

00 1 


E-L equations: 

-m£ 2 QjQ 2 sin (0 X - 0 2 ) - 2mgf sin 0 : 


r 
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dQ- = m^ 2 Q 1 Q 2 sin (0 1 - 0 2 ) - mgf sin 0 2 , 

= 2m f 2 Q 1 + mf 2 0 2 cos (0 : - 0 2 ) , 
dd 1 

= mf 2 0 2 + mf 2 Q 1 cos (0 : - 0 2 ) . 


The equations of motion are now 

0i : 0! + ^ 0 2 cos (0! - Q 2 ) + ~ 0 2 s ^ n (®i ~ 0 2 ) + ~ s: >- n 0i = 0 , 
0 2 : 0 2 + 0! cos (0! - 0 2 ) + 0 2 sin (0 X - 0 2 ) + — sin 0 2 = 0 . 


Do you think you could have derived these equations directly without using a Lagrangian? 


Example 4 

Find the equations of motion for a bead of mass m free to move on a hoop of radius a, 
rotating at constant angular velocity to, as shown below. (Assume the mass can move past 
the point of rotation.) 



Solution: 

We have 

T = |m(x 2 + y 2 ), 


U = mgy . 
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The transformation equations are 

x = a cos cat - a cos (0+ oat) , 
y = a sin cat - a sin (0+ cat) . 

Notice in this case, unlike the other examples, the transformation equations explicitly involve 
the time. Using these, I find 

T = i ma 2 oa 2 + ^ ma 2 (f) + to) - ma 2 oa ^0 + to) cos 0 , 

U = mga sin cat - mga sin(0 + cat) . 


The E-L equation of motion becomes 

0 - ca 2 sin 0 - — cos(0 + cat) = 0 . 
a 

Notice, if we define (J) = 0 + JT and do a rotation in the plane perpendicular to the 
gravitational field, we obtain (g = 0 effectively) 

(fi + ca 2 sin cfi = 0 , 
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which says the bead moves like a pendulum about the far end of the hoop. 

Two Points about Lagrangian Methods 

Let me make several points about Lagrangian methods before we go on. First, the Lagrangian 
of a system is not unique. Consider a new Lagrangian, L' given by 

L ' - L + £ G(q ‘' t] ' (7 ' 2?) 
in one dimension. The action associated with L’ is 


■ - e + J 

tl v------/ 


G q.,t 


(7.30) 


The quantity G(q if t) | is fixed in the variation because (qi)i,2 and t 1(2 are not being 

v 7 i 

changed at the endpoints. Therefore 6S ' = SS and the equations of motion (= Newton’s 
equations = E-L equations) for L' are the same as those with L. 


Another important point is that Hamilton’s principle can be used for systems which have 
potentials very different from Eq. (7.4) above and which are subject to specific types of 
constraints. The crucial assumption in the demonstration that Newton’s equations, (7.1), can 
be written in the E-L form (7.17) is that potentials exist from which the non-constraint 
forces may be derived and which obey Newton’s third law. Such potentials can be time- 
dependent and even velocity-dependent. The form of the relation between such potentials 
and the forces in this case is 


F . = - 


dU 


( d\J ^ 


3x ai 


dt 


dx. 


(7.31) 


Even more generally, there are also E-L- type equations, but not a Lagrangian, for systems 
which do not possess potentials. These are called D’Alembert’s equations. We will pass over 
a description of these types of systems. 


Types of Constraints 

The general form of constraint relations in rectangular coordinates is (j = 1, . . . , m for m 
constraints) 
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2Xi(x,t)dx ai + (fl[(x, t) dt = 0, 


2) W ai( X ' t) X ai + 0)fyx, t) = 0 


Then we may classify constraints as follows: 

df 3 


(1). “Holonomic”: 


0 ) 


: _ 


dx„ 


, to 3 t = 


df ] 

dt 


= 2 


df 3 


df 3 

dt 


df 3 


. Then since 


X ai + - r 

1 r)x ■ (9t 


ai ai 


we may write these simply as 

9 j ( x «i't) = f j (x ai , t) - const. = 0 . 

Also called “integrable” constraints. 


(2). “Non-holonomic”: All other cases. Also called “Non-integrable”. 

As pointed out above, a Lagrangian exists and the E-L equations, Eqs.(7.17), hold for 
systems which have a potentials of the form (7.31). Systems which have holonomic constraints 
and generalized forces derivable from potentials, U(q i ,q i ), have E-L equations have 
the form (7.28). Examples of constraint relations for the three examples which began the 
Generalized Coordinates section are: 


1. Particle on a sphere 

1. g 1 (x lf x 2 , x 3 ) = x 3 + x 2 + x 3 — R 2 = 0 . 
n = 1, m = 1 

s = 3n-m = 2 generalized coordinates (0,4>) 

2. Plane double pendulum 

1- gi = x 13 = 0 

2. g 2 = x 23 = 0 

3. g 3 = x 3 + x 2 - f = 0 

4. g 4 = (x 21 - x n ) 2 + (x 22 - x 12 ) 2 - l 2 = 0 
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n = 2, m = 4 

s = 3n-m = 2 generalized coordinates (0i, 0 2 ) 
3. String of length L through a plane 


1. 

9i 

11 

CO 

II 

O 



2. 

92 

/ 2 2 
— 'V X 21 X 22 

+ X 23 + 

V x n + X 1 2 - L = 0 



n = 2, m 

= 2 



s = 3n - m = 4 generalized coordinates (4>i, §2 r Off) 

All of these are holonomic constraints. A further independent classification is whether the 
constraints involve the time (“rheonomic”) or are time independent (“scleronomic”). It would 
take us too far afield to further discuss these classifications in any further detail. I suggest 
“Principles of Mechanics” by Synge and Griffith for a deeper investigation of these ideas. 


Although I stated the holonomic constraint condition above in terms of rectangular 
coordinates, it is also possible to state such conditions between generalized coordinates 
and to use the constraint formalism of the last chapter. The E-L equations with explicit 
generalized, holonomic constraints (from (6.68)—(6.70)) can be written as 


dU _ _d_ ( 0L^ _ 
dq ± dt (dqj 

9: (9i, t) = 0, 


(i = 1, ...,s) 

> 

(j = 1, ... , m) 


(7.32) 


where 


L = T - U ,U = U+ 2, ^j9j • 

j = i 


(7.33) 


I will call U the “effective potential”. Using the connection between potential and force, 
we see that if one can relate qi to an instantaneous rectangular coordinate Xi by a scale 
factor s such that Xi = qiSjthen the actual force of constraint is 


F = s 3, 


where 


3 


dq ± 


j=i 


(7.34) 
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The constraint formalism is especially useful in revealing the explicit forces of constraint. 

Example 5 

Find the effective force of constraint in the radial direction for the pendulum in Example 
2 above. 


Solution: 

We will write the kinetic and potential energies here as 

T = i mr 2 0 2 , (E.l) 

U = mg (t - r cos 0) , ( E - 2 ) 

where r and 0 are the generalized coordinates and the constraint as (E = constant) 

g = r - i = 0 . 


(Notice how I write the potential in (E.2). If r is changed, this changes only the pendulum’s 
length, not the zero of potential.) The only new equation here is (note g = g(r), so the 
constraint does not affect the 0-equation) 


3L 

dr 


_d_ /Wj 

dt (dr J 


21 


= 0 , 


where I have used L* = T - U* and identified S r from (7.34). Thus (using r = i only 
after the partial derivative is taken) 


3 r = -mg cos 0 - mf0 2 . 


Notice the overall minus sign: the first term counteracts gravity, the second counteracts the 
centrifugal force. This is a general result, applicable even if the pendulum were driven. In 
our case, we know that E = T + U is a constant. Thus 

E = — mf 2 0 2 + mql (l - cos 0) . 

Substituting for 0 2 in S r , we find 

£s r = mg (2 - 3 cos 0) - ^ . 

Now since 
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E = mgf?( l-cos0 o ) , 

where 0 O is the maximum value of 0, we can write 
3 r = mg (2 cos 0 O - 3 cos 0) . 

This is easily seen to be the actual constraining force (due to the pendulum attachment) 
on the mass in the radial direction. Notice S r = -mg for 0 = 0 O = 0, as it should. 

Example 6 

A car accelerates (or decelerates) at an instantaneous rate, x. Find the generalized force on a 
tire at the point of contact with the road if there is no slippage. What does this generalized 
force represent? 



Solution: 

For this problem, all we need know is the kinetic energy of the tire, 

T = — I0 2 + — mx 2 , (E.3) 

2 2 

(I is the moment of inertia of the tire) and the fact that the tire rolls without slipping on 
the pavement: 


x = R0 


(E.4) 


In writing (E.3) I am anticipating the result that a rigid body’s kinetic energy can always 
be written as translational plus rotational energies. The E-L equation in 0 is just 


3L 

d 

( 3L 

30 

dt 

(30 

- 

T - 

lx 

l 

R 


0 


(E.5) 
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where I have used the constraint equation, (E.4), after the partials are taken. Clearly, the 
generalized force, is just the torque exerted on the wheel, the force being given by 
R 'S#, R being the scale factor, S. This simple result is applicable to other situations once 
the instantaneous acceleration x is solved for (often by solving additional E-L equations). 

Endpoint Invariance: Multiparticle Conservation Laws 

I had mentioned in the last Chapter that the endpoint variation considerations there would 
be useful in uncovering conserved quantities. There, we were assuming invariance of the 
action, S, under endpoint variations and finding the consequences. Here we will need to 
directly show this invariance as the first step. This will be done for a model system. 


First, let us examine two types of endpoint variations. Then, from (6.83) and (6.87) we can 
write for the case of s generalized coordinates, qi (using our new "6" notation), 


6S = b L S + 5 e S , 


(7.35) 



(7.36) 


Download free eBooks at bookboon.com 


253 





MODERN INTRODUCTORY MECHANICS 


LAGRANGIAN AND HAMILTONIAN MECHANICS 


5 S 



2 * 



S 


8t + 2 


i = l 



(7.37) 


"SiS" denotes the interior contribution and "S e S" represents the endpoint part. We will 
assume that the E-L equations hold in the interior. Thus, although 5q i (t) ^ 0 one has 
5iS=0 from (7.36). 


Let us look at two specific types of varied paths. 

Case 1: 


i = 1 


r 


(r < s) 



In this case we are looking at the special case of variations, at endpoints as well as interior 
points, 


Sqi = Sq , (i - 1#... f r) 
fit = 0 , 


(7.38) 


for a subset (r ^ s) of the original s coordinates, fiq is just a constant, infinitesimal 
parameter. Then, for this type of variation 


S»’S - Sq 2 


dL 


(7.39) 


The second type of variation is illustrated below. 
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Case 2: 


i = 1 



In this case at the endpoints, 


6q. = 0, 
St * 0 . 


J 


(7.40) 


Therefore, 


S< 2) S = St L - 


l<u 





(7.41) 


Now let now apply these considerations to our model Lagrangian from the beginning of 
this Chapter. It was (a = 1, . . . ,n; i = 1,2,3) 


L(x ai ,x ai ) = T - U, 


(7.42) 


U 

T 


2 u «p(K - x pI) ' 

a < p 


(7.43) 

(7.44) 


We are assuming the E-L equations for this system hold. We notice that L and therefore 
the E-L equations do not change under the translation of coordinates, 


X «i X ai + 8X i ( a11 «)1 

t -*• t, J 


(7.45) 


where the 8x ± are three arbitrary, infinitesimal, time-independent parameters. This is a 
“Case 1” variation with Sx ai = 5x i , 8t = 0 ■ Moreover, the total variation in the action, 

t 2 

6 (1) S = f dt {l (x ai + Sx ± , x ai ) - L (x ai , x ai )} , (7.46) 
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vanishes by the invariance of L. We thus find 


6 (1) S = = 0 , 


(7.47) 


where 


8 “ s - 2 ^ 2 ^ 


<9L 


a = l \^ai / tj 


(7.48) 


Since the are arbitrary, this means 

^ 3L i v» dL i 


a ^ai kl 


-2 


dx . !t 2 

a ai 


(7.49) 


or 


= Pd 


■ti 


(7.50) 


where 


xi 

>*-2 


m a X ai F 


(7.51) 


is the system’s total linear momentum. Linear momentum is conserved as a consequence of 
invariance of L under translations. 


Now let us write this same Lagrangian in cylindrical coordinates (p, 0, z ) The kinetic 
energy is 


T = \ 2 m « (p« + + *») • 


(7.52) 


(7.53) 


This time we notice that L is unchanged under the rotation, 

0 a -** + S 0 f (all a)] 

t -*• t , } 

where 60 is a constant parameter, while keeping p« and z a fixed. This is again a “Case 1” 
variation with 50 a = 50 , 5t = 0 . As before, the total variation vanishes as a consequence 
of the invariance of L, showing that 6^’S = 0 The consequence this time is that 


Y <3L I y 0ij I 
Zj Zj -in It, ' 


<3L 


(7.54) 


a a 


a a 


or 


y m„Pa 0al = y m aPa 


(7.55) 
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expressing conservation of total angular momentum about a particular (arbitrary) axis. 
Angular momentum is conserved as a consequence of invariance of L under rotations. 
Although I used a specific Lagrangian in this argument, any Lagrangian which is invariant 
under spatial translations will conserve total linear momentum, and any Lagrangian invariant 
under rotations will conserve total angular momentum. 


Finally, we notice that the model Lagrangian (7.44) is also unchanged under the substitutions 


x. 


x. 


t + St . 


X ai(t) = x ai(t + , 

(all a, i) 


(7.56) 


where dt is again a constant, infinitesimal parameter, this is a “Case 2” variation with 
Sx t(i = 0. The total variation in S is, 

t 2 + St t 2 

SS = f dt ' L (x ai , x ai , t ') - f dt L (x al , x ai , t) 

t x + St tj 

t2 , 

= f dt { L (x ai (t + St) , k ai (t + St) , t + St) - L (x ai (t), x ai (t) , t)} , 


= st r dt —, 

J at 

ti 


(7.57) 


which vanishes under (7.56) above as long as the time does not appear explicitly in L. 
This means 


5 <2) s = s< 2) s = o. 


and therefore 


( 


\ 


-1 


X ai 


dL 

< 5 X ai/ 


tl 


•-E 

\ i/ a 


we define the “Hamiltonian” as 

dL 


H 




ai 

i,a ai 


dx. 


- L . 


The above implies that H is conserved: 

dH 

- = 0. 

dt 

In our case, not surprisingly 


dL 
dx . 


= m « X ai ' 


dL 

d Ki / 


(7.58) 


(7.59) 


(7.60) 


(7.61) 


(7.62) 
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H = 2T - (T - U) = T + U 


(7.63) 


dH 

(It can also happen that - = 0 but that H * T + U.) Again, although I have used a 

dt 0 L 

model Lagrangian, any L which is invariant under time translations (— = 0) will have 

dt 

a conserved H. Notice that (7.59) says 

dL 


X ai 


dx. 


- L = const. . 


(7.64) 


which is just a statement of the integrated form of the E-L equations, with a change in 
notation from the last Chapter. The condition for this to hold was — = 0, consistent 
with my observations above. 

We may write the endpoint variation for a given Lagrangian, L (q i , q i , t) , as (i = l,...,s) 


6 S = 


-H6t + ^ p L bq L 


(7.65) 


where the general definitions are 

H = " L ' 


Pi s 


dL 

dq ± 


(7.66) 

(7.67) 


(Note: the Pi are generalized “canonical” momenta, which may not correspond to actual, 
physical momenta.) We may repeat the argument in generalized coordinates to show: 


‘Case 1”: variation 

ir dL 

If- = 0 


dq ± 


and that 


‘Case 2”: variation 


p L = const. , 


(7.68) 


If* - 0 
dt 


H = const. 


(7.69) 


These statements may also be shown directly from the E-L equations (representing Newtons 
equations) without invoking the idea of endpoint variations of S. 
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Consequences of Scale Invariance 

Before I go on, I would like to point out that there may be other types of invariances 
in situations that can give us information about the system without actually solving the 
dynamical equations. Let’s take the example 

U(r) = Cr n , (r = |x|) (7.70) 

representing a particle in an external potential. Let’s let (“scaling”) 

r —* ar. (a a real number) (7.71) 

then 

U(cxr) = a n U(r). (7.72) 


Let’s try to make T, the kinetic energy, transform similarly. 

1 — 2 

T = - mr . 

2 


(7.73) 


r -»• ar,t 



2 -2 +n 


a a 


T . 


Therefore L = T - U 


T n T 

L -» a L . 


(7.74) 


Since L is changed only by an overall constant, the equations of motion from 

8 f L dt = 0 (7.75) 

will be unchanged! Let’s look at three cases to see what specifically this implies. (Works for 
any portion of a trajectory) 


(i) Kepler problem, n = -1. 



T °c a 


Kepler's 3— law! 
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(ii) Free fall, n = 1. 



=> T oc Jh . 





(iii) Simple harmonic oscillator, n = 2 . 



T 

=> x independent of amplitude. V 


When Does H=T+U? 

I mentioned above that H is generally not just T + U. Under which conditions does the 
equality hold? We have 

dL _ dT _ dU 

3q± dq ± 5q i (7.76) 
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Now given 


T 


\ 2 ' 


X ai = X a 




we 


have 


x. 


= Y 

dq. 3 dt 


Then the kinetic energy has the general form, 

T = 2 a * b ) + 2 ^ ta't) ^ + c t)' 


j,k 


-t2 


m. 


^ X ai dx ai 


2 H ''Tj dq 

b ) - 2 


m 


: -ik 

^ X ai «5 X ai 


b “ dq, at ' 

1 N' _ ( ^ X cti \ 

■ 9 2 m « v at) ' 


a,i 


Note that ajk — a kj Now we have 


dT 


= 2 a ^fti ^ + % 8 *0 + 


3 <3i j,k 

= 2 2 a ik Tk + b i r 


which means that, 


2 Ir ^ = 2 2 a * + ' 

i j/k j 


and using (7.80), 


2 Jr - 2T - 2 b i % - 2c • 

i i 


Finally, therefore from (7.66), 


H = 2T - (T - U) - 2 I + ^ I ^ - 2c 

tTu q± ’ 


(7.77) 

(7.78) 

(7.79) 

(7.80) 

(7.81) 

(7.82) 

(7.83) 


(7.84) 


(7.85) 


(7.86) 


Download free eBooks at bookboon.com 


261 














MODERN INTRODUCTORY MECHANICS 


LAGRANGIAN AND HAMILTONIAN MECHANICS 


Thus, the conditions under which H = T + U are quite different from that for H = 

dx . dU 

. A sufficient set of conditions for H = T + U is seen to be — — = 0 and - 

at <5^ 

dEl 

Investigation into the Meaning of - = 0 (Optional) 

dt 

dE 

One can separately inquire as to the conditions giving - = o> where E = T 4 

ax . dt 

convenience in the following, I will assume that —— = Q- We have 


at 


= 2 a i k + U ( q i' q i' t) ' 


jk 


'Term (1)" 

"Term (1) 1 


dE 

dt 






a jk q,q k + z , + 


j, k jfk dqi 

However, the E-L equation is explicitly given by 


dU 

dt 


"Term (2)" "Term (3)" 

a (T - U) d (d (T - u)' 

aq ± dt { dq L t 

“Term (2)”: 


= 0 


dr _ y da jk 

d ( 


q,q 


j^k 


2 ?““* 

= - 2 Z ( a ikqk + a ikqk) 


- “ 2 1 44 , - 2 2 a i* 4 , 

k,( ^1 k 

Substituting in (7.89) gives, 


"Term (4)" 

au „ v ^ik • • „ v •• d l au ^ 

a,, q ‘ q ‘ - - 2 2 357 q ' q ‘ - 2 ? ■ q “ + * y = 0 




Isolating Term (4) above, multiplying by throughout and summing on I gives: 

da. 


2 2 qkTi = - 2 tt 1 g^kTi - 2 

i,k i, j,k °^li i 


(au 

d 

' au V 


dt 

U<3± JJ 


Ti 
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This is just an alternate form for Term (1) above. Substituting, (7.93) into (7.88), we now 
have 



(7.94) 


(7.95) 


resulting in our final form, 

dE _ dU d dU . \ 

dt dt dt dq ± 


(7.96) 


(If we had not assumed 


dt 


= 0; there would be additional non-canceling terms on the 


right hand side of (7.96).) So, the condition for E = const, is much less general than 


the condition for H = const. 


— = 0|A sufficient set of conditions for = 0 are 

dU ' ^ dt 

d Xai _ -— = 0 and - = 0. Not surprisingly, this represents a union of the sufficient 

at at aq ± 


conditions for H = T + U 


ax. 


= 0 , 


au 


= 0 and those for H = const. 


dh au . dx . 

— = 0 => — = 0 since —— 

V at at at 

in mechanics than E. 


at dq ± ) 

_ ■ H is undoubtedly a more useful quantity 


) 


Hamilton’s Equations 

There is another formulation of mechanics which also comes from Hamilton’s principle. 
Lagrange’s equations are second order in time derivatives. But mathematically, we know 
that a single second order differential equation can always be replaced by two first order 
equations. This is a hint that another formulation of mechanics exists which yields first 
order differential equations. 

Let’s say we have found a “proper” (unconstrained) set of generalized coordinates, qi. The 
key to the new formulation is the replacement of the q± (t) as the dependent generalized 
variables, with a set which has twice as many members: 
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Lagrangian 
L(q ± f <3i ; t) 


Hamiltonian 

H(q if P ± ; t) 


related by 


8q ± = 



independently varied 



As pointed out above, the “canonical momenta” are in general not “real” momenta. The 
reasons are: 


1. The qi, are not generally rectangular coordinates. 

2. The potential, U, may be velocity dependent. Then, even if we use a 
rectangular coordinate, x, and T = — mx 2 , 


Px 


dL . dU 
— = mx - — . 
dx dx 


(7.97) 


Keeping these facts in mind, the Hamiltonian is 


H = 


2 PA " L ' 


where q k is considered a function of and Pk. Solving for L in terms of H, the action is 


S 



J PA 


H (p^ h) . 


(7.98) 


We will look for the stationary path, this time by varying Pk r Tk simultaneously and 
independently. 
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Variation of the action gives 

SS- fdt 2 

Now 


ti k 


. s . dti . dti _ 

6p k q k + p k Sq k - — 8 Pk - — 8 q k 

<5Pk <5qk 


(7.99) 


8 <3k = — 8 Tk r 
dt 


as before (it is qk and Pk which vary independently, not qk and q k ), so 

^2 i 2 r 2 

/ dt p k 6q k = p k 6q k | - f dt p k 6q k , 

ti T 


8S - / dt£ 


ti k 


fiPk fq k - - aq* I Pk + 


dH 

3q k 


We require 


(7.100) 


( 7 . 101 ) 


6S = 0 


( 7 . 102 ) 


which results in k = 1, . . . , s 

OH ... 

q k = - < (gives p k , q connection: q = q (p.,q.)) (7.103) 

dp k k k k i i 

p k =- r (dynamical part) (7.104) 

oq k 


These 2 s set of first order differential equations are called Hamilton’s equations. They 
replace the Lagrange equations, which may be written as 


Pk 


<9L 

oq k 


(7.105) 


Hamilton’s formulation is not as useful, practically speaking, as the E-L formulation. 
However, they are just as important, if not more so, theoretically speaking. Hamiltonian 
methods permeate all active research areas in theoretical physics and help provide a formalism 
bridge from classical mechanics to quantum mechanics. 
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First of all, let’s make sure that Hamilton’s equations imply the E-L equation, (7.105), 
above. Notice that the first equation, (7.103), reduces to an identity 


= s [ 


- L(q if q ± ,t) 

dp k \ 

i 

q + ^ 

P . ^ 

^ <9L dq i 

4k ^ ^ 

i 

1 dp k 

i dPk ' 


Pi 


cance 


1 71 






(7.106) 


The second equation, (7.104), now gives 
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’‘ = -4(? pA - L(qi 'H 

v p + — + y — 

i 1 <5q k <5q k i dcji 3q k ' 


Pi 


cancel 7S 

p k = — • 

d q k 




(7.107) 


dH 


In this form, it is particularly convenient to calculate — directly. We find 

dt 

dH ( 3H . OH . ^ OH 

* ■? ^ qi + s: Pi r* 


(7.108) 


Using Hamilton’s equations for c^and p ± now shows that the two terms in the first sum 
cancel, and we have 


dH OH 

dt at 

Also notice, using the definition of H and the E-L equations, 
dH V i • • •• \ 

- - 2 (m* + pU 

K K 

cancel cancel 


(7.109) 



v 

( aL . 

aL .. 'l 

aL 



-2 

i 

Ur q± 

+ ^T qi J 

at * 

(7.110) 

dH 

aL 





dt 

at ’ 




(7.111) 


Thus, if H ^ H(t) (or equivalently) L ^ L(t) ) then H = const. 

Here is a cookbook summary of Hamilton’s method: 

Step 1: Find T(q if q if t) and U(q ± , q ± , t) and form L(q if q if t) = T - U in 

the usual way. 

Step 2: Solve for the generalized momenta from 

0L 

p. = — ; —. (terminology: pj is 


aq 


'canonically conjugate" to 
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Step 3: Form H from 

h = J Pi^i - L (qi' t), 

i 

remembering that cp = q i |q jf p jf tj in general. 

Step 4: Write down Hamilton’s equations, 

. _ dH 

- ~ ' (definition) 

fj Pk 

p k =-. (dynamics) 

3q k 

There is a shortcut possible in this procedure. If one can establish that H = T + U (see 
(7.86) above), then we need only eliminate the q ± variables in terms of the Pi and q± 
when constructing H . Notice that if a particular coordinate, qk, does not appear explicitly 
in H, we have 

dti 

P K = (7.112) 

which implies the momentum conjugate to qk is conserved. (The same conclusion also holds 
in Lagrangian mechanics when L * L (q k ) Such coordinates are called “cyclic.” 

Let’s look at two examples of this process. 

Example 7 

Formulate the undamped pendulum in Hamiltonian mechanics. 

Solution: 

We will use the angle, 0 , as the generalized coordinate as before. I have written down the 
transformation equations between x, y and 0 in Example 2 above. These equations do 
not explicitly involve the time. In addition, the gravitational potential is velocity independent, 
which means from (7.86) that 

H = T + U, 


H = — mf 0 + mgf(l - cos 0) . 


Download free eBooks at bookboon.com 


268 





MODERN INTRODUCTORY MECHANICS 


LAGRANGIAN AND HAMILTONIAN MECHANICS 


0 is not an appropriate variable to have in H. We eliminate it by the definition, 

p 0 s ^ = m f 2 0, 

9 ae 


H = 


2m t 


+ mgf(l - cos 0) . 


Hamilton’s equations are now 

3H 


0 = 


Pb 


3p 0 mi 


,2 ' 


(definition) 


and 


® e ■ a 

p 0 =-= - mgr sin 0 . 

30 

Of course we can recover the correct second order equation, 


0 = 


Pb _ 9 


mr 


= - — sin 0 , 


by combining the two. You can see the Lagrangian technique is more direct in this case. 


Example 8 

Given the (nonrelativistic) Lagrangian for a particle (charge q, mass m) in an external 
electromagnetic field (using Gaussian units), 


L = T - U 

1 a 3 

T = - m X i i' u = t)-X k i A i (Sf t), 

2 i c 

( ())(x, t): “scalar potential”, A(x, t): “vector potential”). Derive the Hamiltonian and the 
Hamiltonian equations of motion. IsH = T + U? Is H a constant in time? 

Solution: 

Notice that the potential is velocity-dependent, so we can not write H = T + U as in 
the last example. The canonical momenta are given by 

dL . q 

P, - T- - - A, 


Download free eBooks at bookboon.com 


269 





MODERN INTRODUCTORY MECHANICS 


LAGRANGIAN AND HAMILTONIAN MECHANICS 


We now construct H as 


H = 



- L , 





+ q4> - 




After the algebra settles down, we find 


H 



+ q<t> / 


which confirms that H * T + U • Hamilton’s equations are 


x k = 


= — (p k - — A k I , (definition) 
ap k m\ c 


and 


<9H 


P k = - — = - 


q V („ q ^ \ <9A i _ ckp 


2 k- 


— A. 


<3x k me " c 1 / <9x k “ dx k 


- q 


These equations are a disguised form of the familiar Lorentz force equation for charged 
particles. To see this, form 


mx k 


Pk 



r 


mx ± 
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Using the definition of the curl from Ch.l, 


S x (V x a)]. = 2 (V X A) k 

= J *j V , K ' 


jrk 

l^.m 


k sum 


(6 ,5. - 6. 6J 

\ if jm lm jr. J 


■ 2 - v jAl ), 


mx k = q E k + — (x x b) , (Lorentz force law) 

c ' /k 

1 dA t 


where E k s — V k cJ) -- , (electric field) 

c dt 

B k s (v x a) • (magnetic field) 
dH — 

Notice that - ^ constant if 4> and/or A have explicit time dependence. 

dt 

Holonomic Constraints in Hamiltonian Formalism 

Remember how we handled holonomic constraints in the Lagrangian formalism k= 1, . . . , s 


dL* 


d <9L* 


dq k dt l <9q k 


= 0, 


L* = T - U* , 

m 

U* = U - J Mi r 

i = l 

where the constraints were 


(7.113) 


(7.114) 


Ji(qjr t) = o . (i = 1,... , m) 


(7.115) 


The incorporation of the constraints as a modification of the original potential gives us a 
hint as to how to proceed in the Hamiltonian method. A natural hypothesis is that the 
constraint Hamiltonian equations are just ( k = 1, • . . , s ) 

dH * 

apT' (7-iis) 


<j. - 


Download free eBooks at bookboon.com 


271 









MODERN INTRODUCTORY MECHANICS 


LAGRANGIAN AND HAMILTONIAN MECHANICS 


Pk 


dH* 



where 


(7.117) 


H* - H + 2 k.g. , 

i = l 

gj(qift) = o. (7.ii8) 


(7.116), (7.117) and (7.118) represent 2s + m equations in 2s + m unknowns. We 
can easily verify that this formulation is correct by reducing this set to the equivalent 
Euler-Lagrange equations. We need only repeat the steps leading to (7.107) above. There 
are essentially no new aspects of this reduction beyond what was seen above. I’ll leave it to 
you to show that these reduce to the modified E-L equations 


Pk 



(7.119) 


where 


<9L* _ <9L 

d <3k 3q k 


(7.120) 
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CHAPTER 7 PROBLEMS 

1. Consider a system of two masses in one dimension, one ( m i) attached to a wall 
with a spring (spring constant ki), the other mass ( 1 ^ 2 ) attached to the first 
mass by another spring (spring constant k 2 ). 

a) Formulate the Lagrangian of this system. 

b) Write down Lagrange’s equations and show that one obtains the earlier Eqs. 
(3.62) of the notes. 

2. a) In a given inertial system (“the laboratory”) the Lagrangian of an 

unconstrained particle is given by 

1 

L = 2 m x 2 - U ( x ) 

Write down the Euler-Lagrange equations of motion. 

b) The laboratory in (a) is inside a train, moving at a constant rate with respect 
to a fixed point on the Earth, 

X(t) = vt, v=const. 



What is the Lagrangian of the particle, L in the Earth’s frame of reference? Show that 
although L ' 7 ^ L, one obtains the same Euler-Lagrange equations as in (a). 
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3. A small particle of mass m is constrained to move on the surface of a cone with 
opening angle ot as shown uner the influence of gravity in the z-direction. 



Find the Lagrangian and Lagrange’s equations in some set of unconstrained variables. 
Show that these equations may be solved. 

4. Show, using the transformation from (Xi,yi) and (X 2 ,y 2 ) to 0! and 02 (see 
figure), that the kinetic and potential energies of the double pendulum can be 
written as 


T = m£ 2 0^ + jt 2 Q 2 2 + cos(0i-0 2 ), 

U = -mgf ( 2cos 0i + cos 02). 
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5. System 1: 



System 2: 



(Note: all the masses above are point masses.) For the two systems depicted above: 

a) Explicitly write down all the constraints in the form 
g j (Xi, t) = 0 for both systems. 

b) Find a set of generalized coordinates to describe each system. Explain your choice. 

6. A point mass m is constrained to move on a thin wire while being rotated 
at a constant angular velocity. It is attached to the center of rotation by a 
massless spring with spring constant k and unstretched length f 0 - There is no 
gravitational potential energy in this problem. 
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a) Find the equilibrium distance of the mass from the center of rotation. 

b) Write down the Euler-Lagrange equations for this system. Show that the 
motion is simple harmonic. What is the frequency of oscillations? 

7. (Harkening back to prob.17, Ch.2) A particle of mass m on a frictionless wire 
follows the trajectory y = Kx 2 (K is a constant) in the x y plane under the 
influence of gravity. 



a) Write the Lagrangian of the system, find the E-L equation and find the 
period of motion for small oscillations. [Hint: Drop terms that are small in 
the E-L equation and get it in harmonic oscillator form.] 

b) By introducing a Lagrange multiplier, find the force of the mass on the wire 
in the y-direction. (The answer is: 


-m 


( 2Kx 2 + g \ 
[l + 2K 2 x 2 j 


•) 


8. A particle of mass mi on a frictionless table is attached to another of mass 
m 2 by a rope. The mass m 2 is hanging off the table as shown (mi ^ m 2 in 
general). 



frictionless 

pulley 


Formulate this as a Lagrangian constraint problem. 

a) Lind the kinetic, potential energies and the Lagrangian equations of motion. 

b) Lind the acceleration of mass 2 and the tension in the rope. 
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9. Consider: 



The pulley is free to rotate either clockwise or counterclockwise. It’s moment of inertia is I. 

a) Set up this problem with an undetermined Lagrange multiplier to determine 
x the acceleration of mass 1. 

b) The Lagrange mutiplier will turn out to be related to the tension in the rope; 
it’s exact meaning depends on how you formulate the problem. Actually, since 
the pulley is massive, the tension in the rope on the mi and m 2 sides are 
different. Find the tension on both sides of the rope, TT and T 2 

10. Consider a mass m constrained to move on a wire that is being rotated at a 
constant angular velocity as shown (no gravity in this problem!): 



Find the force of the wire on the mass. [Flint: Introduce a Lagrange multiplier, L associated 
with the constraint.] 
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11. A particle of mass m moves under the influence of gravity along the spiral z=k0 
r=const. where k is a constant and z is the vertical axis. 


z-axis 



a) Obtain the Hamiltonian equations of motion (use cylindrical coordinates 

r, 0, z . ) • 

b) Finish this equation: z = ? ? 

12. Often the KE is a quadratic function of the generalized coordinate, 

T = f(q)q 2 , 
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where f (q) is some function of the generalized coordinate, <3. Often, too, the PE does 
not depend on q: 

U = U(q,t) . 

Construct the Hamiltonian, H. Under these circumstances is H equal to the total energy? 
Is it a constant of the motion? Explain. 

13. Given the general form of the kinetic and potential energies of a system of 
particles ("a" is the particle label, "i" is the axis label), 

T = ^ ' U = U ( x «i ' * ai ) ' 

a, i 

(Notice the potential energy is velocity dependent) construct the form of the Hamiltonian. 
Is this H conserved? Is it equal to T+U in general? 

14. Consider a mass m constrained to move on a wire which is being rotated at a 
constant angular velocity as shown (no gravity in this problem!): 



a) Find L = L(r, r, t) and construct H = H (p r , r, t ) . 

b) Is H a constant of the motion? Is it equal to T+U ? 
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15. Consider a point mass, m, constrained to moving on the surface of a hollow 
cylinder of radius R under the influence of gravity (no friction). The radius 
R = R (t), is a given function of the time. 



a) Find the Lagranian and Lagrange’s equations of motion in cylindrical 
coordinates. 

b) Construct the Hamiltonian. Is H a constant in time? Is H = T+U? 

16. A particle of mass m subject to gravity is suspended from a simple pendulum of 
fixed length E as shown. The whole system is placed in a box and accelerated 
with acceleration "a" along the x direction. 



a) Choosing 0 as the generalized coordinate (be careful!), write down the kinetic, 
potential energies of the system. Show that they are explicit functions of the time. 

b) Construct the system’s Hamiltonian. (You do not have to find the equations of 
motion.) Is H = T+U? Is H a constant of the motion? Explain your answers. 
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17. A point mass m acted on by gravity is suspended from above by a wire whose 
length is a given function of time: f (t). 



a) Find the Hamiltonian for the system and Hamilton’s equations. Eliminate Pe 
and show that the equation of motion for q is: 


^ 0+ —— sin0 = 0. 


0 + 2 

£(t) £(t) 

b) Is the Hamiltonian a constant of the motion? Is it equal to the total energy, 
T + U? Explain your answers. 


18. Remember the bead of mass m attached to a rotating hoop (in the x y plane) 
of radius ' a ' (Example 4) from Ch.7. 



Treat this as a Hamiltonian problem now. 

a) Find the Hamiltonian in the appropriate variables. 

b) Is H = const.? Is H = T + U ? Explain your answers. 
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19. Let’s think back to one of the points made in Chapter 4. There I stated that 
systems which possess a Hamiltonian have phase spaces for which the volume 
neither expands nor contracts. This was called Liouville’s theorem. Show that for 
Hamiltonian systems, 

V • v = 0 , 

_ dx _ 

where v is the generalized “velocity” ( - ~) and V is a generalized gradient, thus 
proving Liouville’s theorem for these systems. [Hint: Set = qi for i=l, . . . ,n 
and x i = Pifor i=n+l, . . . , 2n and use Hamilton’s equations for x and p.] 

20. (A challenging problem.) A particle of mass m slides on the top of a hollowed-out 
half-cylinder of radius R and mass M. The half-cylinder moves frictionlessly on a 
horizontal table, as shown. (The total momentum of the system is given to be zero.) 



a) By introducing appropriate generalized coordinates, find the Lagrangian of the 
system and Lagrange’s equations. 

b) Find the frequency of small oscillations (0 «l) of m. 
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ENDNOTES 


1 We are assuming that A, 1 exists, which we know requires det( X) ^ 0 ; we will establish this momentarily. 

2 Here it is necessary to use the Leibnitz rule for differentiation. 

3 There is even another way to state this. Using the curl concept from Ch.l, V x F = 0 is another way of 

specifying a conservative force field. 

4 Agrees with the result in Goldstein, p.525, when variables are changed appropriately. 

5 For more on parametric representations, see “Calculus of Variations” by Weinstock, ps. 34-36. 

6 Even this is not the most general type of variation! One could also alter the meaning of the 
dependent variable, X, such that x —> x + Sx (x) along the path. In the application in the 
next Chapter, X will be the time variable and there will be nothing gained by considering such 
scale change variations. 
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